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FOREWORD 


This report presents a part of the results of research studies entitled 
"Analysis of Fire Spread for the Improvement of Fire Test Methods", which 
is being conducted as a continual project of BRI for 1979-1983. 


The development of engineering design method of firesafety in buildings 
based on the quantitative prediction of fire phenomena is one of the signi- 
ficant aims of modern fire research. This approach has been successful in 
the fields of fire protection design of structural members and smoke control 
in high rise buildings, and it is now being applied to the development of 
firesafety design for preflashover fires. However, preflashover fires invol- 
ve so many physical/chemical processes that a number of processes consisting 


fire growth are still left unformulated. 


Fire plume, the process which supplies smoke and heat from fire source to 
ceiling layer and sometimes accelerates fire spread along wall if fire source 
is lacated near wall, is counted among the processes that have the most sig- 
nificant influences on firesafety in early stages of compartment fires. 
Though considerable studies have been made on fire plumes, they are almost 
limited to those far from fire source in semi-infinite space. This report has 
focused on the behavior of plume near to fire source, which has significant 
role in the transfer of energy and mass from sire source and lower air layer 
to the ceiling smoke layer. In this report, a mathematical model of the plume 
near to fire source was derived on the basis of recent experimental works, 
and some diagrams to calculate practical properties are prepared. The results 
induced herein should contribute to improve the engineering design method of 


firesafety for preflashover fires. 


March, 1982 
K. Kamimura, Director General 


Building Research Institute 
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CHARACTERIZATION OF THE INTERMITTENT FLAMING REGION OF THE UPWARD CURRENT 
ABOVE DIFFUSION FLAMES 
by 
Yuji Hasemi* 
Building Research Institute, Ministry of Construction, Tsukuba, Ibaraki. 


(BRI Research Paper, No. ) 


ABSTRACT 
The study herein aims at characterizing the intermittent flaming region of 
upward current above diffusion flames, and consists of the similarity ana- 
lysis of the region, derivation of the analytical solution of the basic 
equations for axisymmetric upward current due to buoyancy with respect to 
the experimental facts on the region and the determination of experimental 
constants which are involved in the solution by dimensional analysis and 
experiments. 
As the growth of fire and spread of smoke depend significantly on the be- 
havior of the upward current above fire source in early stages of compart- 
ment fire, the characterization of the upward current above fire source is 
a problem of increasing importance. However, though considerable informa- 
tions are available on axisymmetric upward current, most of them are almost 
limited to the region far from fire source. This report analyzes the: inter- 
mittent flaming region among the near regions to fire source, and derived 
a mathematical model of this region which enables it to predict the proper- 
ties of the region at the comparable accuracy to that for the region far 
from fire source. The report also presents some diagrams for the practical 
calculation of the useful properties such as mass flow rate in the current 


and entrainment coefficient at every height in the intermittent region. 
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I. INTRODUCTION 


In the early stages of a compartment fire, the growth of fire and spread 
of smoke depend significantly on the behaviour of the upward current above 
fire source. However, though considerable information is available to the up- 
ward current, the systematic studies have been almost concentrated to the re- 
gion far from fire orca. Recognizing that the region far from fire 
source is usually involved in the hot gas layer accumulated underneath ceil- 
ing, it is presumed quite naturally that the most meaningful regions of upward 
current in enclosure are those close to the fire source. 

The first systematic measurement of the deterministic properties, such as 
mean vertical velocity and mean excess temperature, in the regions near fire 
source was made by Terai and Niera 2? They have used methanol pool fires as 
heat source, and have found that the upward current in the lower regions is 
almost columnar shaped and the centerline velocity is almost constant with 
height. The measurement of the deterministic properties in the regions close 
to the flame or within it was followed by ecaeerey.” and Cox and Chitty? 
They both have used the diffusion flame by natural gas on a porous refractory 
burner, and have clarified quantitatively how the deterministic properties 
depend on height and heat release rate. According to their works, the upward 
current above diffusion flames appears to divide into three regions, which 
- are distinguished by the normalized height, z/Q7°. The lowest region, z/Q2h 
<0.08, corresponds to the continuous flame, where the centerline excess tem- 
perature is constant with height and the current is accelerated. The middle 
region, 0.08<z/Q%<0.20, corresponds to the intermittent flaming region, where 
the centerline excess temperature is in proportional to the inverse of height 
and the centerline velocity is constant with height. The centerline properties 
in the highest region, z/Q25>0.20, satisfies the dimensional relationships of 
the conventional models of thermal plume. McCaffrey and Cox-Chitty also found 
the dimensional relationship between the centerline properties, height and the 
heat release rate at fire source as described below. 

a A =A(z/¥5)”, 0. =B(2/Q@5) “7+ ey 
where n=1/2 for the flaming region, n=0 for the intermittent flaming region 
and n=-1/3 for the plume. The features of the highest region coinside well 
with the classical model of the plume by Yokoi?? etc, and those of the lower 
regions can explain the results of Terai-Nitta qualitatively. 


Independently from the works of McCaffrey and Cox-Chitty, the author 


= 


sa 


has analyzed Terai-Nitta's results, and has estimated that the eddy coeffi- 
cient in their upward current be tn focn He has also derived a solution of 
the equations of the axially symmetric plume with the uniform-eddy-coefficient 
hypothesis and a few dimensional assumptions derived from Terai-Nitta's data. 
Though this solution explains the most features of the data of McCaffrey and 
Cox-Chitty for the intermittent flaming region, it has not been clarified how 
the eddy coefficient is formulated with the known parameters. Recently, the 
author has compared this solution with the experimental relationships of the 
upward current in this region, and has formulated an equation to predict the 
eddy Coeeeieren ta 

This paper reviews the experimental relationsips on the upward current 
above diffusion flames, and then summarizes the derivation process of the so- 
lution of the plume equations for the intermittent flaming region. The deriva- 
tion process of the solution is modified from that shown in Ref.7), because 
the dimensional relationships derived by McCaffrey and Cox-Chitty are here 


considered. 


2. EXPERIMENTAL RELATIONSHIPS 


In order to find a clue to modeling the upward current in the intermi- 
ttent flaming region, we will here review the empirical relationships from 
previous experimental works. Tab.1 outlines the works here studied, two of 
which use natural gas diffusion flame as heat source. Terai-Nitta and we both 
have used the methanol pool fire; the former settled metal screens around the 
fire source for cutting of the draft, whilst the latter has conducted test in 


a 30x20x15m calm space without obstacles against air movement. 
2-1. Flame Height 


According to the experiments of McCaffrey and Cox-Chitty on natural gas 
fires, the flame height is proportional to Q75, and the centerline properties 
are normalized with z/Q2%. This suggests that the behaviour of the upward 
current above diffusion flames is essentially dependent on the flame height. 
Though the flame height depends on fuel types, we can demonstrate that the 
diffusion flame by natural gas and that by methanol pool fires are alike in 
the relationship between the flame height and the heat release rate. 


McCaffrey has correlated the flame height with heat release rate as 


L =0.08Q/5 (continuous flame) 


L =0.20@5 (including intermittent flaming region) Be 


ey 


As the intensity of the fire source in the early stages of fire is the order 
of Q10°v10*kW, the above formula means that the upper limit of the intermi- 
ttent flaming region may reach 1.0vl.5m in the preflashover stages of fire. 
This height is comparable with that of the neutral zone of a fire compartment 
installed with popular openings. 

The height at the middle of the intermittent flaming region above metha- 


a ; : ES. -f 3 
nol pool fire is formulated as a function of pan diameter by Terai-Nitta as 


Ly =3.5d4/5 ( 3) 


As it is expected quite naturally that Q is proportional to the square of the 


radius of the tray of burning methanol, the above relationship leads to 
Ls =0.10Q2/ ( 4) 


In connection with this, the mean visible flame height in our case is correla- 


ted with the heat release rate as 


L. =0.13Q2/5 (5) 


Comparing these data with those on natural gas fires, we find that the metha- 
nol pool fires and the natural gas fires are alike in the relationship between 


the flame height and the energy input. 
2-2. Centerline Properties 


As mentioned in the Introduction, McCaffrey and Cox-Chitty have correlated 
the centerline properties with the height from fire source and the heat release 


rate as 
n 2n-1 
w/QUS =A(2/V5)", 9, =B(z/QI5)™ ( 6) 
where A and B are constant in the respective regions. McCaffrey and Cox-Chitty 
have further pointed out that we /gB6,z(=A*/g6B), which is clearly constant in 
each region, takes the almost constant value irrespective of the regions. The 
experimental relationships concerning to the centerline properties are summa- 


rized in Tab.2. 


2-3. Transverse Characteristics 


’ 


The function characterizing the horizontal distribution of excess tempe- 
rature and vertical velocity have been approximated conventionally by Gaussian 
distribution. The horizontal spread of the upward current has been usually 


characterized by the radial distance zt or Tr at which velocity or excess 


Gham 


temperature takes the value of 1/e of that at the centerline. Here, we inves- 
tigate into the relationships between the transverse characteristics and the 
centerline properties referring to the other experiments on thermal plumes. 
Tab.3 summarizes the ratio of r3/r* for the previous experimental works 
on both of the intermittent flaming region and the conventional plume region. 
According to these data, it is clear that the width for velocity is always 
wider than that for excess temperature in the intermittent flaming region. In 
contradistinction to this, the relationship between the two widths is rather 
disordered in the plume; that is, the temperature profile is wider than that 


2) 


of velocity in the experiments of Rouse-Yih-Humphreys and Yokoiten whilst 


the velocity profile is wider than that of excess temperature in the works of 


4 and MeOa firey However, a certain order should hold 


George-Tamanini-Alpert 
between the centerline properties and the plume width ratio , ri/r*, if the 
Similarity condition is satisfied in each region. 

Approximating the horizontal distribution of velocity and excess tempera- 
ture by Gaussian distribution as 


W -az"-exp(-gr°/z "), 9 =b-2°™! -exp(-o. 49/2" 9) . {ail 


then, the integrated form of the momentum equation for the axisymmetric plume 


becomes ~ 
i 2 Si “ i-s 
dpmyar = (enya’z™ _ (-n)Wa? " _ geBz"? _ 80207"? ( 8) 
dz 0 ao, 46, ah 26 


From the above relationship, we obtain the relationship between ra/re and 


we /gB0 2, which characterizes the axial characteristics, as follows. 
a 20 


a aay pmad ese ye i «2 
g80,2  (-n)o, ~ I-n*9/ Tw ( 9) 


For the intermittent flaming region, we substitute n=0 into ( 9), and obtain 
w-/gg9 i= 2(r*/r*)* (10) 
a a 9’ Ww 
and for the plume region, substituting n=-1/3 to ( 9), we obtain 
ye 9.z=1 5 (r*/r*)? 11 
q/ 880, hak Ce (ri) 


Consequently, if the similarity condition and Gaussian distribution are both 


assumed wait. LS expectedit nat rh z re holds for we /g80 2 22 in the intermittent 


flaming region and rh mts holds for w2/gB0,z z1.5: in the plume region. Fig 2 


compares these relationships and experimental results, where is seen that 


McCaffrey's data on the plume region does not satisfy the above relationship. 


Inferring from the above considerations, it is expected that the temperature 
profile will become wider than that of velocity in the region higher than the 
area McCaffrey has observed, if we / 280.2 >1.5 holds always. Fig.2 also demon- 
strates that we /g89 2 is larger than 2(rg/r*) in the many of the experiments 
on the intermittent flaming region. This suggests that Gaussian distribution 
may not be the most likely function to characterize the horizontal distribution 
of the properties in this region. 

In connection with the scaling of the transverse characteristics, Cox- 
Chitty have shown that the width of the upward current can be also normalized 
with Qs, That is, the observed values of the width at every height can be 


represented by the functions in the following forms. 


et (2 /02>) (12) 


r4/Q 25 =£, (2/Q) (13) 


Here, we will show that (12) and (13) can be derived from the conservation of 
heat flux, the similarity requirement and the hypothesis that the axial proper- 
ties are represented in the form of (1). First, we assume that velocity and 
1+3n_1-3n 

Z ) and 


) from the second and the third conditions 


excess temperature are represented as w=A- (2/Q75)".QUS .w(r2/x 
0=B- (a SY aa a 


above mentioned. Then, the conservation condition of heat flux is given by 


> S(i-znjyo wltsn re, .,. ti el My et RED 2 + one lone dr 
Q =2C 0Q “Xx Jia B:W(r°/x z )°O(r°/x Zz )Ts5n, 1-30 
3(1-2n)/5 _1+3n,. 
RTH aoe F°A-B-W(n) “0(n)an (14) 


For this equation to hold, the power of Q in the right hand side should be 
equal to that in the left. Recognizing that the underlined part is constant, 
we can clarify how the parameter,x, depends on Q. Now, let x be proportional 


to a Then, the condition for (14) to hold is 
3(1-2n)/5 + £(1+3n) = 1 (15) 


The value of 2 satisfying (15) for n=0 is 2/5. On the other hand, & is not 
necessary for n=-1/3. Substituting x«Q2/S and ( 6) to the functions represent- 


ing the horizontal distribution of excess temperature and velocity, we obtain 


ino A és Wee /xt tong} STG ys) o Wigs /e” ee nba en (16) 


2 


o(r eee ak, 3 Binte coke eae) 


(17) 


i 


where x'=x/Q*9, 2'=z/Q25, r't=r/Q25. (12) and (13) are immediately obtained 
from (16) and (17). 


3. ESTIMATION OF THE EDDY COEFFICIENT IN THE INTERMITTENT FLAMING REGION 


Here, we will derive a mathematical model of the eddy coefficient in the 
intermittent flaming region. This procedure will employ not the statistical 
properties but only the deterministic ones, because it is difficult to measure 
the statistical properties in a flow with extreme gradient of temperature at 
present. The model of the eddy coefficient developed here involves one unknown 
parameter in common with the parameter, c, for the plume formulated with the 
concept of mixing length hypothesis by Woou This parameter will be quanti- 
fied in the following chapters by comparing the solution of the basic equations 


of the upward current with the experimental data. 
3-1. Estimation Method of the Eddy Coefficient in an Axisymmetric Flow 


In order to quantify the eddy coefficient in an axisymmetric flow using 
the measured deterministic properties, we will first derive an explicit cal- 


culation method of the eddy coefficient, K_, for one property — , whose con- 


oe 
servation is described by (18). 
og Uae aaa, 
i Fem hk ain, 2 an (KEQy? ‘ oF (18) 


We consider the balance of — in the region as seen in Fig.3. The regiong’ vy" 
is a space lying between the axis and one stream line. It is evident that the 
transfer of € through the lines AD and BC is caused by convection and the 
transfer through CD is caused by turbulence. Now, taking the limit of the 
balance of € at Az +0, we obtain the balance equation of — in the thin plate 
lying between the axis and the stream line for the stream function w* at each 
height as follows. 

r(p*) x (pr) 


d f 0& 

== rwédr =K_-r(p*)— + i! G_rdr (19) 
dz’ 9 E A eats) 0 6 

Then, the explicit equation to calculate K will be 

r(p*) NE EM 

—! rwédr - f[ Guar 

Kp ep SoO Se (20) 
6 r(e*) ob? ar 


Applying (20) to the equations of momentum and energy for an axisymmetric 


plume, we obtain 


ee 


PLY G 
oe rw dr - J, gBerdr 


dz 
K (r,2) = eee (ye) -ay/ir oe (21) 


r(p*) 


dz’ 0 


Ky(r52) = “T(yF)-00/or ie 


where r(~*) in (20),(21) and (22) means that r should be replaced with yp after 


2 


rw* and rwO are integrated on r. 


3-2. Quantification of the Eddy Coefficient in the Intermittent Flaming Region 


To estimate the distribution of the eddy coefficient by the approach above 


described, it is necessary to assume the distribution of w, € and G. from the 
experimental data. Considering the dimensional relationships discussed in the 


previous chapter, ( 7) can be modified as 


w =A-QWexp(-2,2°/9%2) (23) 


8 =B(z/QU>Jtexp(-2,r°/Q%z) (24) 
The stream function is then given by 
os Bio k Qiaz 2 
p =f rwdr tii exp(-¢ rT /Q#°z)rdr = en ey eae /Q2hz) } (25) 
W 


(1) Eddy Viscosity, K 
Eddy viscosity in the intermittent flaming region can be quantified by 
substituting (23) and (24) to (21). The first term in the numerator of (21) 


can be quantified as follows. 
r Yr 2 u/s 
f zwar ={ A QPexp(-26,2°/Q2) rar = Ae 1-exp(-2¢,,27/072) } (26) 
W 


Replacing r with € using (25), 


d ee eae 1 = A-QYy(1- 2 
yes r a0, ( ee } Qe ( AQdea eae) 


Then, differenciating (27) by z, 


r(w) ce : 
ua rwear Eee (28) 
dz 9 Q% 72 


The other terms in (21) will be quantified by using (23) and (24) as 


Z i 75y1 2 geBQuhe 2 
J gbordr =] B8B(z/Q ) exp(-¢,r /Q%z)rdr = “5c t1-exp(-t4r /Q26z)} (29) 
Iv) 


== 


-2AS sé 2 

erga CXR (ir (04a) (30) 
Q'\z 

The eddy viscosity is then 


r-ow/or = 


; g6BQY{1-exp(-¢,r°/Q%z) }/2z, = 5 W122 
7 


Z wh “ay 
2Aq exp(-¢ r /Q%z)/Q’z 


(2) Eddy Diffusivity of Temperature, Ky 
In the common way as the derivation of the eddy viscosity, the eddy coe- 

fficient for the energy conservation equation is quantified by substituting 

(23) and (24) to (22). The first term in the numerator of (22) is transformed 


as follows. 


a 2 eB OR? -B-Q3h 
A-B 2 vat Be OQ Oe 
f rwedr 2S ore eMac /QVPz}dr = PA CREAY banana io /Q2823] 
(Coie 
A+B-Q3 20,0 w 70 W 
= —— +—{1-(1-— =) } (32) 
Racer te) A Qe 2 
r() 27s 
i werdr = - aA texp (-c gr /Qh2) (33) 
The denominator is quantified as 
r°90/or = -2Br or°- exp (-car°/Q5z)/QS2" (34) 


Then substituting (33) and (34) to (22), we obtain 


ee = = DAE {1-exp(-5 r/Q25z) } (35) 
Do w2C att AtwSer Mg 


3-3. Dimensional Structure and Simplification of the Eddy Coefficient Formulae 


The central problem in modeling the eddy coefficient is to clarify how it 
depends on the heat release rate at fire source and to find a formula describ- 
ing its spatial distribution. To study this problem, we here separate the for- 
mulae of the eddy viscosity and the eddy coefficient for temperature, (31) and 
(35), into two parts; the one representing its dependence on the heat release 
rate and the other representing its spatial distribution. Normalizing r and z 
in (31) and (35) by Q%, then we obtain 

. 2 2 TO ae 
gBB{l-exp(-c,r' /z')}/2c,-A {l-exp(-c.r! /z')} /4c 
| aN eed: Fath, lei 2 ae Se (36) 
2A rex (- rt hey ee. 
ow Sat ee 


-— ALE! {1-exp(-¢, 81/2") (37) 


ACW Sar 


K. = Q¥s 


The above equations means that K and Ky are proportional to Q*%. The under- 
lined part, representing the spatial distribution of K, and Ky, may be too 
complicated. Recognizing here that the functions characterizing the horizon- 
tal distribution of w and © are no more than practical approximations, an 
approximation with a simpler expression is desirable for the more essential 
modeling of the eddy coefficients for solving the basic equations of an axi- 
symmetric upward current. According to Fig.4, representing the calculated re- 
sults of K /Q% and K,/Q*, K, is almost uniform in the intermittent flaming 
region. K is not so uniform, but it is not so different from Ky in the same 
region. Thus, the eddy coefficient in the intermittent flaming region will be 


approximated again by 


K «2 K, ~ k+Q3/ (38) 


where K is about 0.0030.004 [m2/s: kw3°] in Terai-Nitta's. data. The value. of 
k2k/Q35 will be quantified again in Chapter 5 by comparing the solution of 
the equations for the upward current, which will be derived with the condi- 


tion (38), and the experimental data shown in Chapter 2. 


4. SOLUTION OF THE EQUATIONS FOR AXISYMMETRIC UPWARD CURRENT WITH THE UNIFORM 
EDDY COEFFICIENT HYPOTHESIS 


This chapter derives a solution of the basic equations of the axisymme- 
tric upward current with respect to the experimental relationships shown in 
Chapter 2 and 3. The principal problem here is to derive the formulae repre- 
senting the horizontal distribution of w and 8, because the dependence of the 
centerline properties on the height and heat release rate at fire source has 
been already clarified. The basic approach is to derive the ordinary differen- 
tial equations of W(n) and O(n), where near2/xz, from the basic equations of 
the upward current and then to solve them. 

First, we rewrite the momentum equation using stream function instead of 
u and w, and decrease the number of variables. 


2 3 2 
BW nob _ ae 20, ap _ 30°) Cs ze m 39 
“are Sex Lees dr oroz K& Se Pal a2? eek (39) 


As the centerline velocity in the intermittent flaming region is constant with 


height, the form of ~ becomes 


wig if} i} nN ; 
v Fe nie =J raW(n)dn sazJ xW(n)dn/2-=az9(n), a(n)sf xW(n)dn/2 (40) 


Using $(n) above defined, (39) will become 


Lot 


sa 


"3 3 326 
a® 38 92 Ae 299d eos 4 30° dP BUS 
a°(o y 25 ee “9 z=) ‘ a or or i taa7) musts ar3” ror a or2 
+rtgB6 (41) 


Furthermore, applying the following transforms for obtaining the ordinary 


differential equations of $(n), 


007 dO" dn 4 4 ei des, Sona. 


or dn or XZeaH r dn 
370 "ana $49" Nqey2ty5Rdae de, 


or. =e Re ee ee hn 
a3o _ mel nee 0° d°n dn , dd 99n _ 8n9d3o | 12n? d?o 
or 3 =! ees or2 or dn or? x dno r3 dn2 
er dt ie 
oz dn oz he ee z dn 
376 _ 4% an an g7ae 32n iy tntyd “hr 2n do (42) 
aroz. dne. oraz dn oroz rz ddansaft TZ an 


then, we have 


2 
4a2n2z ; +8ak n? al a i$ + rtgpe = (43) 


For solving (43), it is necessary to derive the ordinary differential equation 


of O(n). The energy conservation equation is rewritten with stream function as 


a, de 
7 "an 


K 
90 . 2anz do 00 = 2 ae (44) 


-O) 57 r2 dn oz or 


Considering that the axial temperature is propotional to the inverse of z, 
we can derive the following ordinary differential equation of 9(n) from (44). 


dé 


qe ay gaussian . & dol 
ors Tan -9) 5° ‘ Aiden AS eral s re (33) 


where is applied the following transforms. 


or Ze Ahh ate ya alg 

JD b dO on b do ted. 

OZ co ne ion) oe eas dnzete cia 

2.98, _ 4bn do | d*6, _ 4bn d_, do 

or’ or rz ‘dr. Sane fore dn "dn (46) 


(45) can be transformed further into 


d do 
ance Kean + aod) =0 (47) 


(47) means that the value in ( ) is constant at every height. Considering that 


do/dn=0 and @=0 at n=0, the integration of (47) on n gives 
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2KgndO/dn + a0 =0 (48) 


Then, the relationship between 0 and © will be 


_ 2KQ" ao 
aQ dn (49) 


Substituting (49) into (43), then we obtain 


Heed? | oBbr* . «sgBbx2 
Kane an=  San-z2° 8a 30) 


where K =K = Ko. As the right hand side of (50) is constant, we find the ano- 


ther relationship between © and 4 as follows. 


we 8aK .d2o 
~~ gBbx2 dn2 ; (51) 


Substituting it into @9), we can eliminate © from (49). 


d36 d26 _ : 
Nand * HOgpz =O, usa/2K meg 


Recognizing that the n-th derivative of f(x)=6x/(x+o) is given by 


n 
et ee sree niséo 
dx” ( 1) (x+o) "41 (535) 


we find the solution of (52) in the form of 


3n 


"= atta) Sou 


Then, the stream function is quantified by substituting this into (40) as 
follows. 


3n 6Kzn 
= az@ = 2K: ypz———~. = 
v Punta) nto (S5) 


@(n) is also formulated explicitly by substituting (54) into (51) as 


Bakad’o.) . 8aK 60 _ 960K? (56) 


_ = 
gBbx* dn? —s gBbx* ~— (nto)? gBbx* (n+) 3 


From the above results, w and 9 are formulated as 


poi: jye960K* ne an (57) 
z gBx2 (nto) 3z 

wee Lv _ _120k , (58) 
Beet as (ato be 


In order to quantify the parameter x, we apply w and 9 above obtained to the 


conservation condition of the vertical energy flux. 


I~ 


115202nC_pK3 dn 
ere OP 


-o 115202K3r 


anC pf rwodr = 21C pf ———— — dr = —-— 
pl p°? 0g6x3 (#0) 52 Bx? 0 (nto) ° 
288nC_pK3 
ce Pp = 
F EO) Q 
gBo~x 
Then, 
2nCypK3 
c.. depres 
Oo gBQ 
Substituting this into (57) and (58), 
w= __ §8Q — of 
~ “2a Clak 
(n+o) 
3} 
@ = et 


3nCpp Kz (nto) 3 
all 


Recognizing the conditions,n=r*/xz and x«o+, it is evident that w and 6 
calculated results of (61) and (62) are independent from the value of o. 


Applying o=1 for the simplicity, then the main results will become 


= ape c0Onpane | 
aie oe (n+1)2 


i) made eR 
3C,,e Kz (n+1) ° 
= 6Kzn/(n+1) 
X<= Lay2nC eK 3/280 


(59) 


(60) 


(61) 


(62) 


as 


(63) 


(64) 


(65) 
(66) 


Further substituting K=k-Q¥>, we have the following equations, only one un- 


known parameter ''k'' will be quantified in the following chapter. 


gb Qs 
aut bk (n+1)2 


Q7% 
3™CyP kz (n+1) ° 


v = 6kQ¥%z-n/(n+1) 
xX = 129% ¥2nC, ok3/g8 


5. CORRELATION BETWEEN THE MODEL AND THE EXPERIMENTAL RELATIONSHIPS 


(67) 


(68) 


(69) 
(70) 


Though it is evident that the model represented as (67)%(70) satisfies 


some qualitative features of the intermittent flaming region of the upward 


«tts 


current above diffusion flames, we have to study,at least,the following two 
aspects for the further verification of the model. One problem is to clarify 
how the model agrees with the quantitative features of the experimental up- 
ward current, and the other one is to quantify the model of the eddy coeffi- 
cient. In this chapter, we compare the model formulated as (67) (70) with 
the experimental results on some centerline and transverse characteristics, 


and then we estimate K/Q*%° by several approaches. 
5-1. Centerline and Transverse Characteristics 


Here, we derive some centerline and transverse properties which do not 
involve the eddy coefficient from the model, and compare them with those es- 
timated from the experimental data. 

(1) Rate of the Width of Excess Temperature Profile and that of Velocity Pro- 
file 

Defining the width of excess temperature profile or velocity profile 
by the radial distance r* where 6/0, or w/w, takes one certain value ¢e 
(0 <e <1, e=l/e in the analyses of Terai-Nitta, McCaffrey and Cox-Chitty), 


the curve giving the width becomes 


ré = ¥x(e¥2-1)z C7} 
oS Vx (e¥3-1)z (2) 


The rate of rh /ry is then 
redre = ¥(e%8-1)/(é4-1) (73) 


From (73), it is evident that this rate is not dependent on the heat release 
rate nor the eddy coefficient. r5/re calculated for the various values of e, 
seen in Fig.5, demonstrates that it is almost constant for e<0.2 and rp/ry 
is about 0.78 for c=l/e. This value agrees well with the experimental values 
shown in Tab.3. 
(2) Transverse Velocity and Excess Temperature Distribution 

The experiments show that the transverse velocity and excess temperature 
distribution can be approximated practically by Gaussian distribution. We 
here compare the functions characterizing the transverse distribution in the 
model (67)v(70) with Gaussian distribution. 

According to (63) and (64), the transverse velocity and excess tempera- 


ture distribution is given by 


1/ (n+1) (74) 
1/(n+1)° (75) 


W(n) 
O(n) 


ate 


sa 


The comparison of (74),(75) with Gaussian distribution is shown in Fig.6. 
. e ° * * mo * * 22 

The rate of gto a which is correlated with uxit® by o,/ og (rg /r*) is 
about 9/16 in this figure. This expresses well the experimental rate of pro- 
file widths. 

2 
(3) we/gp6 2 

According to the model, one parameter characterizing the force of buo- 
yancy, we / gB6 2 is 1.5. Comparing this with those listed in Tab.2, we find 
that this is a little larger than the results of Terai-Nitta, Cox-Chitty and 


the present experiments, but it is a little smaller than McCaffrey's data. 
5-2. Determination of K/Q3 


As the dimensional structure of K for the intermittent flaming region 
has been determined as K=k -Q3 in Chapter 3, the principal problem in formu- 
lating K is to determine the parameter k. Before discussing for formulating 
K, we have to note something on the difference between Q in the model and 
that in experimental data. 

The energy flux appearing in the model (67) (70) arises from the conser- 
vation condition of energy flux and means the energy flux dissipating by con- 
vection, whilst the energy generated by combustion of methanol or natural gas 
dissipates not only by convection but also by radiation and thermal transfer 
through the methanol pan or the burner wall. Considering this, the dimensio- 


nal relationship on K above mentioned should become 
K(Etk:03°5= eine) 9? = k/e4>- 0 (76) 


where on is the heat release rate by combustion and Q. is the energy dissipa- 
ted by convection. The parameter, @, is then defined as c=Q./Q- Of course, 
this parameter depends on the type of fuel. Considering only the flow chara- 
cteristics, it should be more essential to estimate k'=k/e** than to estima- 
te kK itself. However, as it is difficult to measure Q. and the other heat 
losses independently, the heat release rate estimated from the consumption 
rate of fuel has been conventionally used as the characteristic heat release 
rate or heat flux in many analyses of plume experiments. Also we estimate 
only Kk in this section. 

The value of k can be estimated with the measured excess temperature and 
velocity by transforming the equations on the profile or centerline properties 
to the explicit form on k and then substituting the measured properties into 
it. To estimate k from the centerline velocity or excess temperature, for exa- 


mple, the explicit formula to get k is derived from (67) or (68) as follows. 


apie. 


ame oo am eanta shea A cto (77) 


2) 
Q@.Z. =l 
E Qs i a ( a ) im Ti (78) 
k = —s——_ = =——— = ————_ 
8 3mC 3m7C_oB 
3mCe az T p° gi p° 


The formulae to estimate kK are also derived from the formulae of the width 


of the profile of velocity and excess temperature as follows. 


Ce H/ VP) "28 3 
~ 288mCpp (e¥?-1) 2 (2/Q7)? (79) 
S. (r4/Q*)*g8 1B 


288nC,p (e!8-1)2 (2/Q2B)2 (80) 


The values of k=K/Q% estimated by the above four methods from the data of 
Terai-Nitta, McCaffrey, Cox-Chitty and our experiments are summarized in Tab. 
4. The table shows that the difference of kK due to the estimation methods or 
fuel species is very small except for the case of Terai-Nitta. Taking k= 
0.0025 for the simplicity, the model (67) (70) can be rewritten as 


4 gB Q¥ys 
Ne nC. Heats as 81 
0.0057C 0 (n+1)? (81) 


@ = a"8 
0.00751C,,p2(n+1)* ae 


n= r2/xz, x & 0.002129 ¥nC o/g8 (83) 


Fig.7 shows the correlation between the widths estimated from the model and 
those estimated from experiments. The agreement between the model and the ex- 
periments is good for rx, whilst the experimental data of Fs are somewhat 
scattering around the theoretical curve. Fig.8 shows the comparison between 
the calculated isotherms and those. generated from McCaffrey's experiments. 


The calculated isotherms are more swelled than those by the experiments in 


the high temperature region; the reason for this is thought as a result of the 


ignorance of the compressibility in the derivation of the model. 


6. MASS FLUX, ENERGY FLUX AND ENTRAINMENT IN THE UPWARD CURRENT 


The practical quantities other than the centerline properties are such 
integral quantities as mass and energy flux. In this chapter, the model (67)~ 
(70) are employed to calculate them with respect to the compressibility of 


the-flow. As the consideration of compressibility makes it impossible to cal- 


=35- 


culate some integral properties, we will prepare the calculation diagrams of 


the integral properties as a practical result of this reporty 
6-1. Mass and Bulk Flux 


The mass flux between the axis and the radial distance,r, at the height, 


z, 1s defined as 
M(z,r) = 20f pwrdr (84) 
0 


Employing the model (67)%(70), and considering also the compressibility due 


to the rise of temperature, M(z,r) will be formulated as 


M(zexr) = 20pg{ nemedr = 0. 03ap gO 8 eee (85) 
ol oT +59 ; JOG 61336Q7F/nCpp_2 (nel) Dine 


The right hand side of the above equation may not be obtained analytically. 

To derive a practical calculation method of M(z,r), we divide the both sides 
of (85) by Q. Then, M(z,r)/Q is written as a function of the normalized height 
»z'=z/Q2/s as follows. 


M(z,r) _ nage dn | cs 
Q combi f {1+1338/™CpP_2' (+1) °}(n+1)? a 


The bulk flux in the upward current is also defined and formulated as 


zs "2rw 3/5 dn 
F(z5r)e= anJ wrdr = mxz] eso = 0.037Q z | ane 
: 0.031Q%r2z (87) 


r7+0.00212Q7zVnCpPm/ BB 


Dividing the both sides of (87) again by Q, we obtain the relationship between 
F(z,7))/ Oeands 2! sand). 
Cae 0.03mr'2z! 


Q r'*+0.00212z! TCpPn/ BB 


The entrainment coefficient,a, which is conventionally defined as dF(z,~)/dz= 


(88) 


2aF (z ,°)/r* hale is then quantified as 
r*-dF(z,~)/dz TC,0, Yt 
ty IW 9S ee pea eT aL 
oL SF (z,*) OOM brlay aa) (89) 


While the entrainment coefficient has been assumed as constant in the conven- 
tional analyses of turbulent plumes, (89) means that it decreases with height 


in the intermittent flaming region. 


2f6> 


6-2. Energy Flux 


The ennergy flux between the axis and the radial distance,r, is formu- 


lated as (90), if the compressibility due to temperature change is considered. 


r 
H(z,r) = 27C f owOrdr = pag pie 
4° 1Q dn 
= +08 ———————————E eS (90) 
Hof, {1+1338Q*/mCpo_2 (n+1)?}(n+1) ° 
Dividing the both sides by Q, we obtain 
H(2;r) dn 
I (91) 
Q Les f {1+133B/mCpo_2z' (n+1) 2} (n+1)? 


6-3. Calculation Diagrams of Mass and Bulk Flux, Entrainment Coefficient and 


Energy iE bux 


As M(z,r)/Q, F(z,r)/Q,a and H(z,r)/Q are decisively determined for one 
set of the physical properties of air, we can here prepare the diagrams for 
the practical calculation of these properties. Recognizing that r=r* and r+ 
are usually taken as the characteristic radial distance of axisymmetric up- 
ward current, we show the diagrams of M(z,r*)/Q, M(z,*)/Q, F(z,r*)/Q, F(z,~)/Q, 
H(z,r*)/Q, H(z,~”)/Q and a. 

As it seems impossible to calculate M/Q and H/Q by an analytical approach, 
we have to apply a numerical approach to get them. For the practical calcula- 
tion with reasonable precision, we divide the integration interval into two 
parts; the one is the part where excess temperature is considerably high and 
the other is the one where excess temperature is not so high that the follow- 


ing approximation is reasonably assumed. 
p = p,/(1+80) = py (1-88) (92) 
Then, we can obtain M/Q and H/Q by integrating the first part numerically and 


integrating the latter part analytically. M(z,~)/Q, for example, can be cal- 


culated as follows. 


is} 
M(Z5°) - g.03mp2'Lf © at 
Q Se Uf {1+1338/Cpem2" (n+1) 2} (n+1)* 
+ Ae Pe a dn : : a) 
Na {1+133B/1CpP_2' (n+1) }(n+1) 
M dn = 1, __ 1538 ;_an 


= 0.03 92'[f_ ————»—_._.,*_ 1 : 
E Y {1+1338/mCpepz! (nt1)3}(n+1)* ne Tp? 2! (ntl) 
dn ee 1338 
0 {141338/mCppn72' (n41)3}(nt1)2 Natl AnCpog2! (ngtl) 


ae em (93) 


= 0.031992" nie 


=i72 


(nety2) 


ya 


The part where we employ the numerical calculation is underlined in (93). The 
integral properties calculated with the assumption ,P_=P0/(1+89, ) and 
Pn mae eptnea? # are shown in Fig.9 Fig.11. M/Q is almost proportional to 2/Q? 


and it is practically estimated by the following simple equation. 
M(z,°)/Q = 1.17(z/Q2")-0.022 | (94) 


The agreement between the estimated entrainment coefficient and Terai-Nitta's 


data is not good, but our result coinsides well with the estimate of McCaffrey. 


7. CONCLUSIONS 


The behavior of the intermittent flaming region of the upward current 
above diffusion flames can be predicted practically by the model represented 
by (67)%(70), which is based on one solution of the equations of axisymmetric 
buoyant flow. Practical integral properties, such as mass and energy flux or 
entrainment coefficient can be also obtained diagrammatically from Fig.9 
Fig.11.Though it is an interesting problem why the eddy coefficient in this 
region depends on qi’ it is not necessary to investigate further into the 
structure of the turbulence in this region as far as we are concerned in the 
engineering deterministic properties. The results of this report hold for 
methanol pool fire and diffusion flame above natural gas burner, but it may 
not be applicable to the fuels which are much different in the shape or height 


of the flame generating from them. 
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TERMINOLOGY 


Cp :specific heat[kW s/kgK] 

F :bulk flux in upward current [m3/s] 

H heat flux in upward current[kW] 

K :eddy coefficient [m2/s] 

I, :flame height [m] 

M :mass flux in upward current[kg/s] 

Q :heat release rate at fire source[kW] 

W(n) :characterizing function of the transverse distribution of vertica 
velocity[-] 

d :diameter of methanol fire source[m] 
:acceleration of gravity[m/s2] 
:radial coordinate[m] 

r'  :radial coordinate normalized with Q25[m/kw2] 

re :width of upward current[m] 

u :radial velocity[m/s] 

a’ :turbulent fluctuation of radial velocity[m/s] 

W :vertical velocity[m/s] 

w! turbulent fluctuation of vertical velocity[m/s] 

x :characteristic scale length of plume width[m] 

Zz :vertical coordinate[m] 

ay :vertical coordinate normalized with Q25 [m/kw2 ] 

O(n) :characterizing function of the transverse distribution of excess 
temperature [-] 

6(n) :characterizing function of the transverse distribution of stream 
function|[- ] 

a :entrainment coefficient|[-]| 

B :expansion coefficient [KJ 

0 :densidy [kg/m] 

n :nondimensional parameter on the horizontal distribution of physical 
properties 

6 :excess temperature 

subscripts 

a :axis 

mm :mean 3 


:vertical velocity 


sexcess temperature 
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Fig.8(b) Theoretical and*Experimental Isotherms(Q=4.9kW) 
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Fig.9 Normalized Mass and Bulk Flux vs Normalized Height from Fire Source 
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Fig.10 Entrainment Coefficient vs Normalized Height from Fire Source 
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Fig.11 Normalized Heat Flux vs Normalized Height from Fire Source 
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Table.1 Experiments on the Upward Current above Diffusion Flames 


Heat Source, Intensity Measurement 
velocity excess temperature 


1) C-A thermocouple 


Present pool fire of methanol, 11.8kW §& 26.6kW 
pool fire of methanol 


natural gas, 21.7kW & 44.9kW 


Bidirectional 


tani? 


McCaffrey”) 
6) 


Correlation Method C-A thermocouple 
ee C-A thermocouple 
natural gas, 17.9kW & 47kW Correlation Method thermocouple 


Bidirectional 
Cox-Chitty 


Table.2 Summary of Axial Characteristics of Three Regions in Upward Current 


wo /Ue=a-2 , @=B:2!07 5 I (on 
j 2 
Flame Intermittent Plume nal 
Flaming ————————— 
g 
0.03°<z'r< 0:08 | 0.08 <-zt < 0.20 0.20 < z! 2289.2 


1/2 -1/3 
B 
Present - = 1.90 70 = - | 0.72 
Terai-Nitta he 700vg00 | 1.5 56 2 “ 0.56 
McCaffrey 6.84 797 1.93 63 Tor? 21.6 0.86 
Cox-Chitty 6.83 880 1.85 70 1.08 23.6 0.69 
ease 2 ~ : : : 1.42 29.7 | 1.00 
Yokoi!) : i | : : tay 24.6 0.81 
George-Tamanini”) E : Bok : 1.24 24.6 0.640.69 


Alpert 
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Table.3 Plume-width Ratio for Upper Regions of Upward Current 


* * 
ee 
Intermittent | 
; Plume 
Flaming 
Terai-Nitta - 
McCaffrey 0.83 
Cox-Chitty - 
Rouse-Yih- 
Hees 
Humphreys 
Yokoi Lia 
George-Tamanini- 0.92 


Alpert 


Table.4 Estimation of K/Q¥ from Previous Experimental Works 


McCaffrey Cox-Chitty 


Present 


a eP 0.0028 0.0026 
2mCypA 
1 
3mCppB 0.0025 0.0025 
{ Pye si 0.0027 
288nCpp (e¥2-1)? (z/QP) ~.0..0033 
(rg /Q7°)*g8 ys 
0.0025 0.0024 


*Zeencpp (@ 8-1) 22/Q) 
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FORWARD 


This report presents the result of the reseach on the 
Statistical analysis of fire spread phenomenon, which was 
conducted as the part of continual studies of BRI for 1976 - 
1980, entitled "Design System for Fire Safety in buildings". 

Concerning fire spread phenomenon, many model fire tests 
have been carried out so far, besides real fire incidents have 
also been investigated in detail for many years. Such tests 
results and investigation reports are now being accumulated as 
are possible to be statistically analysed by an electric 
computer. 

memchiasestuday, by using such data,»the effects«of various 
factors on fire spread are statistically analysed and the 
phenomenon is tryed to be predicted for the given factors. 

The result of such statistical studies on fire spread will 
make the basis for establishing the reasonable design system for 
fire safety in buildings, together with the phisical or the 
chemical studies on this field. 

I would like to acknowledge with appreciation those who 


supported on this study. 


May, 1982 
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Tsukuba, Ibaraki, Japan 


(Research Paper No96 May 1982) 


ABSTRACT 


In this paper, statistical techniques which can be applied 
to the estimation of fire spread are introduced, together with 
their application examples to Occurrance of Flashover in Fire 
Tests and Development of Real Fire in Dwellings. 

These techniques are called "Theory of Quantification", 
since numerical values are given to the factors expressed in 
nominal scale so as to estimate fire spread most precisely, based 
on the assumption of a linear sum model of the factors, where 
coefficients are given to the factors expressed in interval scale 
aseweias iComparitson-of ‘the ‘effectsVofHthe “factors can “be 
achieved by comparing the ranges or the partial correlations 
calculated by using the given values to the factors. 

In the application example of these techniques to Fire Tests 
Data, the domain where flashover occurres are shown on the two 
dimensional Euclidian space using linear sum of the facters. 
Concerning Real Fire Data, the discrimination if a fire develops 
up to the phase that a room of fire origin is wholly involved is 
conducted. 

The development 'of isuchanalysi's ‘of: Fire Test) Data or Real 
Fire Data which are accumulated so far will give more valuable 


and direct informations to the-estimation of fire spread. 
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I. INTRODUCTION 


Fire spread is very complex and manifold phenomenon dominated 
by many factors including uncertain ones such as natural 
conditions, fire extinguishment activities, human behabiour,etc. 
Consequently, in order to analyse and estimate fire spread 
enoughly, not only determonistic phisical approches but 
probabalisti¢eorestatisticalapproaches.are needed. 

Many successful Fire Tests and inquiries of Real Fire Scenes 
have already conducted. Tokyo Fire Department is now conducting 
very detailed inquiries of fire scenes including interior linings, 
structures, ignition sources etc...which are available to analyse 
fire spread. In Japan, fire reports by Fire Defence Agency, Tokyo 
Fire Department and some other fire departments are now being 
accumulated as the data base which can be analysed by electric 
computers. 

When analysing such data concerning fire spread, avairable 
statastical.technigues should be examined well. Since fire, spread 
_has qualitative aspects specified by descrete phase development or 
occurrence of flashover besides quantitative aspects such as 
temperature rise, burning rate, flashover time. Similar does the 
dominate factor, that is, such factors as types of interior 
linings, structures or kinds of ignition sources are qualitative 
factors, while, water content of materials, temperature of outside 
air etc. are quantitative ones. "Theory of Quantification" 
developed by statician C.HAYASHI [1] is very useful for analysing 
the phenomena where qualitative variables and quantitative ones 
are thus mixed. 

After statistical multi-valiate techniques including "Theory 
of quantification" are introduced at the following paragraph, the 
results of the application to Occurrence of Flashover in Fire 


Tests and Development of Real:Fire in Dwellings are reported. 


II. STATISTICAL :-TECHNIQUES 


In statistical terminology, the subject to be analysed or 
estimated, that is fire spread in this study, is called dependent 
variable, while, the factors are called independent variables. 
These variables are either qualitative or quantitative. 
Qualitative variables have levels given with nominal scale and 
quantitative variables have values given with interval scale. 

What kind of multi-variate analysis techniques is available for 
some subject depends upon the types of dependent variable and 
independent variables. Table 1 shows the available techniques 
corresponding to the types of variables. All of these techniques 
premise linear sum model of factors, and the effect of each factor 
is estimated by way of maximizing whether the correlation 
coefficient ("a -"e') or the ‘correlation®ratio: ("d=—9@ 44 
Technique (c) covers.(a) and (b), and technique (£) covers*(aj@ane 
(e), so only techniques (c) and (f) are explained below. The 
other techniques can be led easily as the special cases of these 
techniques. 

Either of the techniques premises the following linear sum 


model of factors. 


Table 1 Statistical Multi-variate Analysis Techniques 


Dependent 
Var. 


Nominal 
Scale 


Interval 
scale 


Independent 
Varn 


(a) Multiple Correlation | (d) Multiple Correlation 
: . (b) Theory of Quantifica-— | (e) Theory of Quantifica- 
phelinen eR eras tion type I tion type IL 


Interval Scale & (c) Theory of Quantifica- | (f) Theory of Quantifica- 
Nominal Scale tion type I' tion type II' 
are mixed 


(b),(c),(e),(f) were developed by C.Hayashi 


Ci) a. = 


Rk 
z 
te = 


| 14 
z 6 Gk) Ae erty | ee cee ste 'C 
j=l k=1 d 
where 

Component Variable on sample zt ( sample score ) 


Number of factors expressed with nominal scale 


-(G7k) =({1 : if sample t is”classified into category k 


ae 
L 
R 
a : Number of categories of factor j 
3 
by factor 7 


0 : otherwise 


Sgt : Category Value which is given to category k of 
factor g. This value is estimated by the method 
explained later. 


r Number of factors expressed in interval scale 
Wy Weight given to factor Z ( interval scale ) 
Zo3 Vatue of factor lL of "Sample tr 

C Constant 


The method to estimate Xone and W, and C depends upon the 
type of the dependent variable. 


(1) The case that dependent variable is interval scale. (c) 
Xi and W, are estimated so that Peason's product moment 

correlation coefficient between sample score q and dependent 

variable y may be maximum. The correlation coefficient eva is 


expressed as follows. 


620) p FS 
yo Oo Cs 8} 


where, n is the total number of samples. Y,@ are the averages 


qi-y, a and Oy, Oq are the standard deviations -of*y ,/a’ 


respectively. 

The “condition that Pya is maximum leads to the condition that 
<p 0 0 a s 
ony yea (eet leo sta faasys Verelets cere) » —e- = 0 C2 = 1,2, 
0 Xyyv ) 0 Wz 
*vi51,r -) “are satisfied simultaneously. a 


These expressions lead to the simultanious linear equation 
concerning Xy,, and Wz. This simultanious linear equation is an 


arbitrary one as it is, by the fact that whether some sample is 


sa 


classified into arbitrary-one of the ctegories of each factorsgon 
not depends up on the state of classification into other 
categories. Consequently, this equation is solved generally by 
adding the following conditions, the number of which is R. 
kj 
Cay one Koy ney = 0 Gules Loe; pine ne? 


where, NK is the number of samples which are classified into 
categoryske? by? factonu ae 

Further, the correlation coefficient is constant even 1f @is 
transformed linearly, so constant value C and the scale of the 
solution can be determined arbitrarily. If they are determined as 


Y= @ = C,,and ¢. => , the result of the analysis can be 


Rye y 
interpreted easier since the solution thus obtained is eaqual to 
that by least square method, namely, the solution solved under the 
condition that the square sum of error term, which is given as the 
following expression, is minimum. : 
Cl) S = ; 2 =: : ini dotnet 8 
ve? ae ai 


2 
1 cae i 


L 


(ii) The case that dependent variable is nominal scale : (f) 

This method is a kind of discriminant analysis using the 
correlation ratio as the criterion of preciseness of 
discrimination. X ik and Wz are estimated so that the correlation 
ratio may be maximum. The correlation ratio n*? is the ratio of 
'between class varience' of to the 'total variance’ o*, which are 


expressed as follows. 


2 
Ss 
is nye s 
fe) 
n 
where, oe 1: isi Lome aes) 
Qk t 
T= 
tm 
n 
2 tr, = — \2 
63% 2D —l a aa) 
b pos ry t 


n : Total number of samples 


n, : The number of the samples which is classified 
into category t by dependent variable. 


ix: The average of sample score a of the samples 
belonging to category f. 


ce : The number of categories of dependent variable 


The condition that n* is maximum leads to the condition that 


on n* 
= 0 (Sete eoee nee = leon we ares joo eee 
dX > 3 3 3 Uu 3 gWs 0 
Ct = 1,2.,¢¢¢6e5 mod} are isatisf ted \simulitaniously. 
Substituting the expression (5), 
oO < 
6) 3b | es Te) Ct = 1 6e5crer ,R s VE=silgiek-er sk ) 
OX OX u 
uv uv 
2 
o°b _ 2 307 a BIS 
Z L 
These simultanious equations lead to the eigen value problem 
after all, where n* is the eigen value and ae Laas 
" ; Ret eee ees ik; 


Xo4 aes hn Kix eee XRkp Wy suis cas Wz aula as Wy ] 


is the eigen vector. 


In order to solve this problem uniquely, besides the 
condition (3) is necessary to be considered, Gd S5sandpG.arehalso 
necessary to be determined appropriately since the correlation 
ratio is constant even if @ is transformed linearly. Here, the 
conditions that 9° = n* and C = 0 are added. 

At this eigen value problem, the spectral factorization of 
the matrix concerned give the solutions, the number of which is 
eagul to the rank of the matrix, where the eigen vectors are 
independent with each other. When discriminating the categories 
of dependent variable, namely, when predicting the phenomenon, 
better discrimination can be achieved generally by taking up not 
only the solution which gives the maximum correlation ratio ( 
eigen value ) but some other solutions together. 

After Xx, wz and C were determined, the effect of the 
factors can be compared with the values of partial correlations 
between dependent variable and each factor ( When caluculating the 
partial correlation, G, is substituted for the value of dependent 
Variable, if .it is.nominaliscale?). -=Thefcomparisonzofntheseffects 
only among the factoers of nominal scale, can be achieved also by 


the range of category value X >, within each factor . 


By the way, Theory of Quantification type I' (c) is 
mathematically eaqual to the multiple correlation analysis using 
1-0 dummy variable for the independent variable of nominal scale. 
However, it would be a remarkable point that the comparison of the 
effects of qualitative factors with the partial correlations or 
the ranges is generated by the conception to quantify the 
qualitative variables, not merely to calculate the partial 


regression coefficient of 1-0 dummy variable. 


TII. Analysis of Plashover sby ‘Fane .TesttData 


Occurrence of Flashover can be analysed by applying Theory of 
Quantification type I' and type II' introduced at the previous 
paragraph. 

The Fire Test Data analysed here are those done by Dr. 
F.SAITO, the main purpose of which is to examine the effect of the 
thickness of plywood on smoke generation. Flashover time is also 
observed in this test. Such results of the tests are reported at 
lieterature [3]. 

The model box used was 95cm width x 46cm high x 43cm depth in 
size, and plywood was set on the cieling and the back and side 
walls. Cribs were arranged in the corner of the box. Small piece 
of soft fiber board ( lcm x lcm ) dipped with alcohol was used for 
ignition. The figure of the model box is shown at appendixI.. 

The tests were designed to examine the effects of the 


following factors. 


1s RatetofS0pening -AreastoSFrontoArea S45 =so27hevels 
(LJILZ 42 je172 
2..- Size of Fuel. =s-3S—9e7levets 


(1) 8 cribs, the length of which are 6 cm ( 4 layers ) 
(2) 16 Cribs, the length of which are 14 cm ( 4 layers ) 
(3) 25 Cribs, the length of*which are*17sem: (°5elayerss) 
Each Crib has the section 2cm x 2cm. 

3.1 ThacknéssnotsPiywoodest- = .oaievels 
(1 )23% Ommeod 2)¢.54,5mmy4.89n79 -btamictt)... 9a Ommet5j463 220mm 


Rate of water content of Crib and Plywood were also measured, 


and they were also taken up as the factors of flashover together 
with the above three factors. 


(i) Analysis of Flashover Time 

Theory?) of Quantification,»type jI' can be applied in this case 
Since flashover time, which is dependent variable, is given with 
interval scale, and the factors, which are independent variables, 
are given with both nominal scale and interval scale. 

The result of the analysis of flashover time is shown at 
table 2. 

Comparing—the-valuestof_partial) correlation, the effects of 
Suze of Fuel is (hargest among 5’ factors, and the effect*%of 
Thickness: om PlyYwoods icy nextfito at. |Onjthe othershand, Rate of 
Opening Area and Rate of Water Content have little effect. of 
course such results are based on the levels of factors taken up. 

The category values of Size of Fuel show that the larger is 
the Size of Fuel, the smaller is the category value. This means 
that the larger is the fuel, the shorter is flashover time. 

Seeing the category value of thickness of Plywood, the value 
fore. (mm ao par: icularlly small, therefore you can see’ that 
fieashover timesas particularly short in this case. The relation 
between Thickness of Plywood and category value is not monotonous, 
the value is maximum in the case of 7.5mm. 

Figure 1 shows the relation between obserbed flashover time 
and estimated one. The correlation coefficient which incdicates 


the preciseness of the analysis is 0.785. 


(ii) Analysis including the case that flashover did not occur 
Theory of Quantification type II' can be applied in this case 
by arranging the dependent variables as follows. Namely, 
flashover time is classified into the following five categories. 
4 ijyeloo sec. De low=( 299200 — 299 sec. (3) 800,- 399 sec. (4)=400 
= 499. sech apt5) 600 sec. above. Further, the category (6)!no 
flashover is added. : 
As the dependent variable thus arranged is nominal scale, the 
analysis is done so that above categories may be well 


descriminated. Table 3 shows the two solutions which give the 


Table 2 Result of the analysis of flashoyertime by Theory of 
Quantification type I' 

* AV = 9.4%, SD = 1.7% ** AV = 7.9%, SD = 1.6% 

**xk Weight is calculated using standerdised value of Rate of Water Content 


Loi = (X,, ~ AV, )/SD, 
Category Partial 


Openings 1/4 30 0.089 
1/8 4. seen) Ome eo 0.063 

Crib 6 cm 

Length 14 0.761 
17 

Thickness 3.0 mm . 

of 5.5 
Plywood yee 2.788 0.573 
9.0 

12.0 


Rate of Water 
Content of 
Pl ywood** 


-0.035 -0.024 
( Weight ) 


36 
21 
17 
16 
11 ‘ 
1] : 
9 
Rate of Water 
Content of OU s 0.052 
Crib* ( Weight ) | 
4 
4 
4 


Fig. 1 


Relation between 
Flashover Time and 
Estimated Value 
(Correlation Coef. 
= 0.785 ) 


Obserbed Flashovertime (min. ) 


0 5 10 
Estimated Flashovertime ( min: ) 


largest two correlation ratios ( eigen values ). These two 
solutions, namely axis I and axis II are independent with eath 
other since they are obtained by spectral factorization of the 
matrix concerned. The correlation ratios, 0.770 of axis I and 
0.357 of axis II are fairly good. 

Either the partial correlations and the category values of 
axis I can be interpreted almost similarly with the results of the 
analysis of flashover time. 


Figure 2 shows the distribution of sample score a;. The 


distribution on the direction of axis I shows that the 


Tab.3 Result of Multiple Correlation Ratio Analysis(Quality & Quantity Mixed) 
*AV = 9.4%, SD = 1.7% WAY 217.82, SB = 1 16% 
**xk Weight is caluculated using standerdized value of Rate of Water Content 


3 2), = (X, p> AV,)/SD, 


AXIS I] 


Correlation 0.770 0.357 
Ratio 

Nominal No. of | Category Partial | Category Partial 
Openings 

WxHx1/4 gar -0.022 -0.138 

WxHx1/8 37 0.020 0.042} 0.043 0.123 0276) 10.157 


Z2e390) 0.803 0.852 | 0.454 


i.0364 0.4573 0.949 } 0.381 
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-0.076 -0.136 
Crib* (weignt) | 
Water Con- 
tent of 70 0.058 0.102 -0.216 -0.231 
Pl ywood** (weight) (weight) 


Water Con- 


Figure 2 Distribution of Sample Score a 


Note: Ellipses show 90% confidence bounds of two dimensional 

distribution of sample scores based on the assumption of 

normal distribution. 
shorter is the flashover time, the larger is the sample score, and 
the category (6) no flashover separates remarkably. Axis II shows 
that the longer is the flashover time, the larger is the sample 
score except two categories of (5) 500 sec. above and (6) no 
flashover. 

The ellipses in figure 2 show 90% confidence bounds on two 
dimensional distribution of sample scores based on the assumption 
of normal distribution. As being expected from this figure, the 
prediction .of flashover time or presence of flashoverg¢eanobe 
achieved by determining the domains appropriately on the 


multi-dimensional Euclidian space where sample scores distribute. 
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IV. ANALYSIS OF REAL FIRE SPREAD 


Real fire spread is analysed by Theory of Quantification type 
II in this paragraph. 

Real Fire Data used here are residential fire reports 
inquired by Tokyo Fire Department from 1974 to 1977, the ignition 
place of which are whether living rooms or kitchens. 

The subject to be descriminated is if a fire develops up to 
the phase that a room of fire origin is wholly involved ( Whole 
Room Fire Phase ), the statistical index of which is if the burned 
area of concerned fire exceeds 80% of the floor area of a room of 
PIce Origin. 

Seven kinds of factors were taken up as independent 
variables, namely, (1) Ignition Pattern (2) Rate of Combustible 
Presa onpvertical Faces 3) §6ombusri biti eyeot Cerling. (4) Rate of 
Combustible Area on Horisintal Face (5) Openings ( YAVH ) (6) 
Structure (7) Extinguishment in the Early Stage of Fire Spread (8) 
Avért” Time" to Fire Station.” The~categories of these factors are 
explained invdetalil ~at“appendixmLi. | Mhough four -factors out-of 
these factors are basically expressed in interval ‘scales, they are 
transformed into qualitative factors by deviding values 
appropriately. Such transformation is effectual especially if the 
phenomenon is not linear for the concerned factor. 

The result of the analysisimby Theory of Omantification=type 
II is shown at Table 4. Comparing the range and the partial 
correlation, you can see that, for the occurrence of Whole Room 
Fire, the effect of ‘Alert Time to Fire Station' is the highest 
among the factors, and the effects of Extinguishment in the Early 
Stage of Fire Spread', ‘Ignition Pattern', 'Rate of Combustible 
Area on Vertical Face', and ' Combustibility of Ceiling" are 
comparatively high, however, 'Rate of Combustible Area on 
Horisontal Face' show low effect. 
| Using the! category vlues given) as’ the result of the analysis 
by Theory of Quantification, the discrimination if: a fire develops 
up to the phase that a room of fire origin is wholly involved can 
be achieved. Fig.3 shows the distribution of sample scores. From 


this figure, you can see that Whole Room Fire occurs with the 


higher probability as the sample score increases, further it may 
be also understood that there would be a discriminant point 
whether Whole Room Fire occurs or not. Based on this distribution 
of sample scores, the discriminant point ao= 0.23 is given under 
the condition of minimizing the error of discrimination. Using 
this discriminant point, the rate of success of discrimination is 


Tek 


Table 4 The result of the analysis of Whole Room Fire Occurrence 


“ 1 RANGE*OILES 
5 2 IMFLAMABLE M. ,OIL 
& 3 APP.+FIBER, WASTES ) (OL. & 
= IGNITION 4 TOBACCO-FIB. , WAST. (> 0.316 0.187 
=| =ATIERN 5 MATCH+FIBER,WAST ; 
=a ’ ° 
ke 
6 
g 7 
D— 0.285 0.147 
0.210 0.099 


0 - 75 (3%) 
135)=1 100 


D 
46 | -0.027 ; 


0 - 2.5 (m%) 
2.5 - 5.0 

5.0 - 10.0 
10.0 - 


FIRE SPREAD FACTOR 


0.038 
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( 0.104 
=0.033 : 


2g 
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Correlation Ratio = 0.30 
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-_DISCRIMINANT POINT 
Qo = 0.23 


FREQUENCY 
on 
oO 


SERRE RRR 


| 


-1.0 -0.5 0 0.5 1.0 
SAMPLE SCORE 


Whole Room Fire did not occur 
SaaS Whole Room Fire did occur 
Figure 3 Frequency distribution of sample scores 


For example, if ‘Ignition Pattern" is 'l; range and the likes 
fowor! or fats', 'Rate of Combustible Area on Vertical Face' is 
reo, ‘Combustibiiity of Ceiling  is-*3; combustible 
material', 'Rate of Combustible Area on Horisontal Face' is '4; 
Sage Opening’ is” "3; 6m 4% ', **Structure’ is’ '3; wooden 
erructure’, ‘Extinguishment in Early Stage of Fire Spread" is ‘1; 
being done'. and ,'Alert Time to Station',is '2; 4 minutes',. the 
sample score comes to be 0.36, and the rate.of success .of 


discrimination is 79% .when Whole Room Fire is, considered to 


occur. : 


V. CONCLUDING REMARKS 


It was made out that the multi-variate analysis techniques 
introduced in this paper is useful for the factor analysis and the 
estimation of fire spread, the results of which can be easily 
interpreted in relation to the phenomenon by showing them with a 
diagram. 

The development of such analysis of Fire Test Data and Real 
Fire Data which were accumulated so far will give more valuable 


and direct informations to the estimation of fire spread. 
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ee 
FACTOR EXPLANATION OF CATEGORY 


1 Ignition pattern from small kitchen range,range,or water heater 
to oils and fats (typical ignition pattern at kitchen). 


5 2 Ignition pattern from certain ignition source to imflamable material, 
# oils and fats or explosive substance excluding the pattern stated 
< as 1 
™ | IGNITION PATTERN 3 Ignition pattern from appratus to fiber goods or wastes. 
a seer 4 Ignition pattern from tobacco to fiber goods or wastes. 
3 5 Ignition pattern from match, lighter,candle or stick of incense 
4 to fiber goods or wastes. 
aS 6 Ignition pattern from certain ignition source to furnishings or 
H furnitures. 
7 Ignition pattern except above 1 - 6 or the ignition pattern is 
unknown. 
RATE OF COMBUSTIBLE AREA 1-0) = 25° (Re) 2 -254- 50 
ON VERTICAL FACE * 8) S075 4 75 - 100 ( Fire Outbreak Room ) 
ee : 
) . : 
& 1 Non-combustible material 
a oe eee 2 Fire retardant material 
fay EIL 3 Combustible material ( Fire Outbreak Room ) 
fa 
<< | RATE OF COMBUSTIBLE AREA {10 - 75 ( 8) 
: ON HORIZONTAL FACE * 2°75 = 100 . - ( Fire Outbreak Room ) 
n ; 
10- 2.5 (m%) 2 2.5 - 5.0 
fx] : 
fe OPENINGS ( SAVH ) 3°5-0-- 10.0 4 10.0 - ( Fire Outbreak Room ) 
fa 
-1 Fire proof structure 
STRUCTURE 2 Fire protective structure 
3 Wooden structure ; 
% f 
ae EXTINGUISHMENT IN THE EAR-| 1 Being done 
=e LY STAGE OF FIRE SPREAD*** | 2 Not being done 
AB . 
ae ALART TIME TO STATION **** epi seat vewvne 
* In the case that there are combustibles such as combustible furnitures,they 
are projected to the vertical face ( wall ) and the horizontal face ( floor ), 


and the projected part is regarded as combustible.Next,the projected area and 
the area of combustible interior finish materials are summed up on the vertical 
face and the horizontal face,and the rate of the each summed up combustible 
area to the area of the each face is calculated. 

cad The testing method of combustibilities of interior finish materials is decided 
by the law.The order 1 to 3 corresponds to the order of combustibilities. 

*x*kk Pire extinguishment activity in the early stage of fire spread ( generally 
before flashover ) with fire extinguishment equipments such as fire extingui- 
shers. 

*x*k*xx Elapsed time from the ignition time ( estimated ) to the alert time in fire 
station. 
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FOREWORD 


This paper reports the result of the research on the 
dynamical analysis of deterioration processes of polymeric 
materials, which was conducted as part of continual studies 
of BRI for 1977-1982, entitled, "Effective Use of Inorganic 
Polymers. for Building Materials". | 

Recently the importance of durability of building materials 
and components has been widely appreciated, and the rational 
evaluation criteria for the progress of deterioration have been 
pursued in various ways. A lot of phenomenological data have 
been accumulated. However, there still remain major difficulties 
to be solved for the analytical methods of deterioration processes. 

In this paper, the deterioration process of polymers by 
the simultaneous action of UV-light and diffusive oxygen is 
explained theoretically based on an unsteady-state dynamics, 
taking some experimental data into consideration, and the 
parabolic law concerning exposure time is derived for a typical 
path of the progress of deterioration. 

The result may be incomplete in some respects, as compared 
with the whole problems of deterioration, but this paper will be 
useful as a basis for the evaluation of durability of building 
materials hereafter. 

I would like to express my hope for the development of 
research on this field and acknowledge with appreciation those 


who have supported on this study. 


December, 1982 


_K. Kamimura, Director General 


Building Research Institute 
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DETERIORATION MODEL OF POLYMERIC MATERIALS 


WITH SPECIAL ATTENTION TO AGE AND DETERIORATION DEPTH* 
by 
Toshio Fukushima** 
BUILDING RESEARCH INSTITUTE 
MINISTRY OF CONSTRUCTION 
Tsukuba, Ibaraki, Japan 
(Research Paper No.97 December 1982) 


ABSTRACT 


The deterioration process of polymeric materials by the 
Simultaneous action of UV-light and diffusive oxygen is 
analyzed theoretically based on an unsteady-state dynamics, 
taking account of the progress of deterioration inward from 
surface. It is supposed that the deterioration proceeds 
photochemically by depolymerization and photooxidation reactions 
of polymer molecules photoactivated by absorbing UV-light 
near the absorption band characteristic of polymeric materials. 
As a result it is shown under some assumptions that the depth 
of the deteriorated layer increases approximately in proportion 
to the square root of exposure time. 

This /t law (parabolic law) has previously been observed 
in the cases of neutralization of concrete and oxidation of 
metals, but this theory suggests that the parabolic law should 
also apply to the case of the deterioration of polymeric 
Materials. 

Experimentally many polymeric materials show the power law 
of exposure time (t”, n=0.5~1.0). The difference between theory 
and experiment is considered to be due to the complicated 
deterioration mechanisms which involve depolymerization and 
photooxidation. 

The reduction in flexural strength over time of polymeric 
materials is explained by the reduced contribution of the 


deteriorated surface layer to the strength. 
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I. INTRODUCTION 


Polymeric materials are widely applied in buildings 
such as paints, adhesives, waterproofings, finishing materials, 
thermal insulating materials, and solar collector components. 

Though they have many attractive properties, they are 
subject to gradual deterioration when they are used under 
natural weathering, and they deteriorate rapidly under the 
Simultaneous action of UV-light and diffusive oxygen(say, ozone) 
in accelerated aging tests. 

On the other hand, societal needs for the conservation of 
natural resources and energy has given increased importance 
to the concept of durability as performance over time in many 
fields, especially in building materials, components and/or 
elements[1-4]. At the same time, in connection with the prediction 
of lifetimes of buildings, there has been an increasing demand 
for criteria for evaluating the progress of deterioration and 
for identifying deterioration mechanisms of building materials. 
In order to evaluate the performance over time, it would be 
necessary to consider the change of material properties and 
the effect of the deteriorated layer on the performance, 
assuming no change of performance requirements. 

The deterioration processes of polymeric materials are 
very complicated. They show various types of deterioration 
behavior, depending on the conditions under which they are used, 
because they are influenced by various extrinsic factors(UV-light, 
oxygen, water, pollutant gases and freezing-thawing action,etc.) 
and various intrinsic factors(temperature and humidity conditions, 
types and compositions of polymers, treatment and maintenance 
conditions, etc.). Much phenomenological data on the 
deterioration processes of polymers in buildings has been 
accumulated. It would be reasonable for us to establish rational 
selection criteria for the progress of deterioration 
based on the data. Taking the complexity of the deterioration 
phenomena into account, however, it may be helpful to try to 
understand the deterioration mechanisms of polymers ynder natural 
weathering theoretically by developing an appropriate model 


based on fundamental principles. 

This paper deals with the complex deterioration process of 
polymeric materials by the simultaneous action of UV-light and 
diffusive oxygen based on an unsteady-state dynamical analysis, 
paying special attention to the progress of deterioration 
inward from surface. The aim is to derive the law of the 
increase of the depth of the deteriorated layer with time, and to 
consider the reduction in flexural strength as the deterioration 
advances. It is shown that the degree of deterioration varies 
continously with time and depth and that the increase in the 
depth of the deteriorated layer with time is approximately in 
accord with the parabolic law [5-6], as is often observed in 
the cases of neutralization of concrete [7], and oxidation of 
metals [8]. 


II. OUTLINE OF THE MODEL OF THE DETERIORATION PROCESS 


The fundamental assumption of the model for dynamical 
analysis are as follows: 1) Only polymer molecules which are 
activated into excited states by absorbing the UV-light (photon 
flux, no) near the absorption band characteristic of polymeric 
materials can take part in photochemical reactions (reaction 
efficiency, n; molar concentration, C*). 2) Photochemical reactions 
are i) depolymerization of activated polymer molecules (lst order 
reaction, R\=kK,C*),, 11) photooxidation by the collision sof 
diffusive oxygens: (molar concentration, .Ca;,diffusion flux, a. 
diffusion coefficient, Da), and activated polymer molecules (ane 
order reaction, R2=k,C*C,) 3) Deteriorated polymer molecules 
(molar concentration, Cp) do not diffuse and remain in their 
Original positions. 4) The influence of temperatures on the 
photochemical reactions is included in materials property 
constants suth as rate constants (k,, k2)and diffusion 
coefficient. (Da) according to the Arrhenius law. 

5) Deterioration proceeds inward from the surface and the degree 
of deterioration varies with exposure time (t) and depth below the 
surface (x). As the deterioration advances, deteriorated 


polymer molecules having IR-active functional groups such as -CO, 


-OH are created, and the degree of deterioration is reflected in 
the increase of absorption (a) or the decrease of reflection (R) 
in IR spectra. 6) The deteriorated layers result in a reduction 
in flexural strength of polymeric materials with the progress of 
deteriration. Based on the fundamental assumption decribed 
above, the following three basic processes are considered. 

a) Absorption of UV-light: 

Considering the material balance within the infinitesimal 
plate, x~x + dx, shown in Fig.2, on the assumption that 
UV-light (wavelength, 1; photon flux, nq) falls on the 
polymeric material surface as decribed in Fig.1, the following 
differential equations can be obtained (see Appendix A 
for the derivation of these equations) : 

oe" 7 t= "=n /Ng) On/dxe—.k C* mks C*C, (1) 

on/ox= ~Cof,n (2) 
(No, Avogadro’s constant; Co, initial concentration of 
undeteriorated polymer molecules; Eye molar absorption 
coefficient). Equations (1) and (2) lead to equation (3): 


dC*/dt= ACq€, exp (-Co€,x) mkyC® —koC*Ca (3) 


(A = nno/No) 


b) Diffusion of Oxygen: 

From the material balance of oxygen within the plate 
show in Fig.2, and the following equation is obtained: 

9Ca/dt= Dpad*Ca/dx* —k, C*Ca (4) 
c) Process of Photochemical Reactions: 

Overall reaction rate (R=Ri+R2) is the rate of creation of 
deteriorated polymer molecules, and it is also equal to the rate 
of increase of the absorption coefficient in the absorption band 


of IR-active functional group (say, -CO group, wavelength, i°). 


Consequently, 
R= kKIC* + k2C¥Cy = 0Cg/9t = Va, >/dt (5) 
Here, 
ki= Aiexp(-Ei/kT), k2= Az,exo(-E2 /kT) (6) 
Da= Da, ©XP (-E3/kT) (7) 


a 


(T, absolute temperature; k, Boltzman constant; E,, E2, Es, 


activation energies) 


III. BASIC DIFFERENTIAL EQUATIONS FOR DYNAMICAL ANALYSIS 

For the purpose of simplyfing the analysis, two assumptions 
are made, however, the results never fail to have generality. 
The assumptions are: i) compared with the thickness of the plate 
of polymeric materials (L), the deteriorated depth (6) can be 
considered to be small enough for the thickness of the plate to 
be regarded as virtually infinitely large; ii) in photooxidation 
reaction (R,= k;C*Ca), we can consider the reaction to be quasi- 
lst-order (R,= k, Ca). Taking these considerations into account 
together with the initial conditions and boundary conditions, 
basic differential equations for dynamical analysis are summarized 


as follows: 


9Ca/It= Dad*Cy/dx*?-k,“Ca (8) 
IC*/Ot= AC, exexp(=Cetxx ly mak Cami. 6 Ce (9) 
dC,/dt= k,C* + k. “Cy (10) 


[Initial Condi erons |: t=0; C*=Ca=Cp=0 

[Boundary Conditions]: [I] t>0, x=0; Ca=Ca, 
[II] t>0, x+ +”; Cy=0 
[III] t>0, x=5; 93Cp/ax=0 


IV. RESULTS OF DYNAMICAL ANALYSIS 
Solving the simultaneous differential equations under the 
given initial conditions and boundary conditions by the Laplace 
transformation method, analytical solutions (11)-(16) are 
obtained in the form of dimensionless concentrations (see 
Appendix B for the derivation of these equations) : 
dx=Ca(x, t)/Cay= (1/2) {exp (-xvk, “7Da) erfc(x/2y/Dqat -vkK, “t) 
+ exp(xvk, ~/Dq) erfc(x/2/Dat + vk, 7t) } 
(11) 
greC* (x, £)/C,¢= {1-exp(-k, t) }- {exp (-Cy€, x) — (Ca, /C*) kK, “bat 
| (12) 
Op=Cplx, t)/C *= [k WPAN cle rempdé st) at | 
-{exp(-C,€)x)/k,} + {(k, “-k,)CA,/C,* } 
ee dnat 4 1 -exp (kB) sber On e/ Croton, 
(13) 


V=dp/dt= {l-exp (-k,t) }-exp(-C,€ x) 
— (Cao /Co*) ky “exp (-k, t) ba (14) 
Here, 
Cy *=C, Ae ,=C nny€,/N, 
erfc x= 1-(2//m) S* exp(-E2) dé (15) 
Profiles of functions based on calculations by a computer are 


shown in Fig.3~8. 


V. DISCUSSIONS 
The degree of deterioration varies with time and depth, 
and it is influenced by the photochemical reaction constants 
(k,;, k,~) and diffusion coefficient (Da). In Fig.4, 6, 7, and 8 
are shown the changes of the space-distribution of each 
concentration function with time, and the influence of 
photooxidation and oxygen diffusion on the distribution of 
oxygen are shown in Fig.3 and 5, respectively. Deriving the time 
dependence of surface deterioration, we obtain equation (16): 
oB(0, t)= [k,t -{l-exp(-k,t) }]/k, + (Ca,/Cy*) 
eT ki ska) tots l-exp(-k,t)] 
k,/t (Cay/Co*) 
= k,Ca,t (16) 


This shows that surface deterioration advances approximately 
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in proportion to exposure time, and that the progress of 
deterioration depends upon the surface oxygen concentration 
and the rate constant for photooxidation, though it depends, 
of course, linearly upon the number of photons absorbed, and 
consequently, upon the intensity of incident UV-light. 

That the deterioration proceeds continuously inward from 
surface with time can be observed in IR absorption spectra and 
in the change of materials properties. For example, Watanabe et 
al. [10] have measured IR sbsorption spectra for sliced samples 
of outdoor-exposed organic polymer plate of polystyrene (PS) 
and polyoxymethylene (POM) and other polymers. These experimen- 
tal data are compared with the theoretical results in Fig.9.a) 
and 9.b). The fact that the progress of the deterioration is 


reflected in the change of the IR absorption spectra for the 


deteriorated surface layer has been confirmed by the author 
for inorganic polymers as well [11-12]. Fig.10.a) shows the 
change of the IR absorption for the hardening of a film of 
sodium silicate water-glass (Na,0-3.25i0,) for the surface 
treatment of asbesto-cement board, and Fig. 10.b) for the 
deterioration of films of sodium silicate wter-glass for the 
surface treatment of steal plate during immersion in 50°C 
hot water (drying at 120°C). The IR spectra were measured by 
an infrared spectrophotometer equipped with a micro-computer 
(HITACHI model 260-50) by the surface reflection-absorption 
method using a T1C1-T1Br total reflection prism. 

In this theory, thermal oxidation is not taken into account 
directly as an elementary process. But, from the above 
discussions, we can explain the progress of deterioration of 
polymeric materials by the chemical attack of penetrative 
reagents (mainly oxygen) on the activated polymer molecules 
(mainly by the absorption of UV-light, but possibly also by 
thermal activation). Eurin [13] has given a valuable review of 
the degradation of organic polymeric materials, in which he has 
pointed out the importance of photooxidation and thermal 
oxidation reactions. Though he ingeniously explained the 
progress of photochemical degradation, he did not refer to the 
progress of deterioration inward from surface. The author would 
like to point out that the deterioration progresses continuously 
from the surface, and that the evaluation of the degree of 
deterioration is very important in comparing the durability of 
building materials. 

On the other hand, using the boundary condition [III] i.e. 
Iop/dt| _.=0, we can obtain approximately the time dependence 
of the depth of deteriorated layer (see Appendix C for details 
of the derivation of this equation): 

6= /Bt B=6 (Da tek pF KD) (17) 
From this result, it can be shown that the depth of the 
deteriorated layer increases approximately in proportion to the 
Square root of exposure time. This /t law (parabolic law) was 


obtained for the first time for the progress of deterioration 


of polymeric materials as the natural derivation based on the 
unsteady-state dynamical analysis, assuming the simultaneous 
action of UV-light and diffusive oxygen. Fig.1ll.a) and b) show 
the comparison of experimental results with the theory by the 
least-squares plots. The experimental data are from Kubota et 
al. [14]. It can be seen that many polymeric materilas show 

the deterioration by the power law of exposure time (t™: n= 0.5 

~1.0). The difference between theory and experiment is 
considered tq be due to deterioration mechanisms which involve 
depolymerization and photooxidation. If the deterioration is 
determined by depolymerization, it is expected to proceed 
linearly with time, and if it is controlled by photooxidation, 
the progress of deterioration obeys the power law. Usually, the 
deterioration of polymers occurs by a combination of the two 
reactions, so that the deterioration of many polymers obeys the 
power law. To fully understand the deterioration mechaniss, it 
is necessarymfor us to look at them in more detail. Recently, 
the author has observed the same parabolic law experimentally in 
the deterioration by water of water-glass films [11-12]. 

The parameter 8 which determines the progress of 

deterioration is a function of materials property constants ( 
Da, k, 


evaluated approximately by the parabolic law, if extrinsic 


kaw ), and the progress of deterioration could be 


factors and these material property constants are clearly 
quantified. Further, since these constants are functions of 
temperature as shown in equations (6), (7), the effects of 
temperature on the progress of deterioration can be considered 
to be included implicitly in the parabolic law. 

Such unsteady-state dynamical analysis should be also 
applicable to other types of deterioration processes. The author 
has just applied the same analytical method to the carbonation 
and neutralization processes of concrete [15-16]. 

Defining a deterioration progress factor (X) as the ratio 
of the depth of the deteriorated layer to the plate thickness, 
and considering that deteriorated layers do not contribute to 


the flexural strength of the plate, the following law for the 


reduction in flexural strength with the advance of deterioration 


can be obtained: 


X(t)= 6(t)/L; (02x51) (18) 
F(t)= Ales (t) |? ALA [1-045 (de) vegan 
= F,{1-X(t) }? (19) 


where A is a constant, and F is the flexural strength of the 
plates in the initial undeteriorated state. It is assumed, as is 
often observed experimentally, that the flexural strength is 
proportional to the second power of the thickness of the plates. 
According to Watanabe et al. [10], deteriorated layers near the 
surface of polystyren plates exposed to natural weathering for 
five years, far from not contributing to the flexural strength, 
influence the inner undeteriorated layers, and reduce their 
specific flexural strength. This experimental fact suggests that 
the deterioration progresses continuously deeper into the 
polymer plates than is expected. This point must be studied more 
in detail, and the influence of deterioration on various kinds 
of mechanical strength such as. compressive, tensile and flexural 


strength, must be considered. 


VI. CONCLUSIONS 

The deterioration of polymeric materials by the 
Simultaneous action of UV-light and diffusive oxygen was 
theoretically analyzed based on an unsteady-state dynamics. The > 
theoretical results are compared with experimental data. 

The deterioration of polymeric materilas progresses inward 
from surface, and the degree of deterioration varies with time 
and depth. The increase in the depth of the deteriorated layer 
is approximately in accord with the parabolic (vt ) law. 

The deteriorated layer influences the mechanical strength 
(say, flexural sterngth), resulting in the reduction in the 


total strength of polymeric materials. 
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APPENDIX A: Absorption process of UV-light 

Let us consider a polymer plate of thickness L and the 
cross section S, as shown in Fig.l, and take the x-axis 
perpendicular to the plate surface. 

When UV-light (wavelength, A; intensity, I,,) falls on the 
surface of a plate at an angle 6, a part of the light is 
reflected from the surface at the rate of R(6, \X), and the rest is 
transmitted into the plate. Consequently the energy flux at the 
surface for the UV-light transmitted into the plate is given as 
follows: 

6,,= {1-R(8, A) }1I,,/S (Al) 
Representing the energy flux at the depth x below the surface as 
6,(x) and the linear absorption coefficient as kK) the energy 
absorbed per unit time and unit area within the infinitesimal 


plate (thickness, dx) shown in Fig.2 is: 


- d®,= |, - Oy |ypgy= 7 (db)/dx) ax 
= Kpopixidx . (A2) 
Consequently, the following equation can be obtained. 
d,/dx= -K 7%) (x) (A3) 


kK, is usually constant in an undeteriorated transmitting media, 
but if the deterioration advances, it is considered generally 

to depend upon the depth from surface (x) and time (t), i. e., 
Ky 3 Kae eer Or convenience, however, it is assumed to be 
constant throughout this theoretical model. 

The linear absorption coefficient (k,) can be correlated 
with a molar absorption coefficient (€,) as follows: 

Oo Coe, (A4) 
where Cyhistthegini coral concentration of undeteriorated ploymer 
molecules. 

From equations (Al), (A3), and (A4), the energy flux ata 
depth x from the surface (6)(x)) is: 

@,(x) = 0) )exp(-Cy€,x) (A5) 
Consequently, the photon flux at a depth of x from the surface 
n(x) is obtained, using the Plank constant (h), the velocity of 
light (c), and the wavelength of the incident UV-light (\), is: 

n(x) = 6)(x) /h(c/) 


= FD 


(A), /hc) exp (-C, €)) 


Ho expl=C ec) x) 

on/ox = —C €,n n,= A® s/he (A6) 
where n, is the photon flux at the surface of the plate. 

According to the law of photochemical equivalence, absorbed 
photons generate as many activated polymer molecules, some of 
which can participate in photochemical reactions (reaction 
efficiency, n), while the rest of which are likely to be 
thermally or radiately deactivated into the ground state. 
Consequently, the number of active polymer molecules generated 
within the infinitesimal plate per unit time is n-(-dn/dx) -Sdx, 
the mole number being (nN,)-(-dn/dx)-Sdx (N : Avogadro's 
constant). 

The accumulation rate of active molecules within the 
infinitesimal plate is represented as (dC*/dt)-Sdx (C*, molar 
concentration of active polymer molecules), and the anihilation 
Tate is R-SdwwRy annihilation wate “of active polymer molecules 
and creation rate of deteriorated polymer molecules). 

Considering the material balance within the infinitesimal 


plate, the following equation is derived: 


(9C*/3t)Sdx = -(n/N,) + (dn/dx)Sdx —- R°*Sdx 
so that: 
(3C0%/at). = -(n/N,)\-dn/ox — R hay) 


APPENDIX B: Derivation of dimensionless functions 

Fundamental differential equations are summarized as 
follows, together with the initial and boundary conditions: 
[DIFFERENTIAL EQUATIONS] 


BCa Jats aDnoGCa/oX=—aak Cp (B1) 
ac*/at = AC, c)exp(-Cye yee kee - k, “Ca (B2) 
sCp/st = k,C* + k, Cx (B3) 
[Initial Conditions] ; t20,.x>0: Ca=CB=C*=0 (B4) 
[Boundary Conditions]; [I] t>0, x=0, Ca=Ca, (B5) 
| [II] t>0, xe: Ca=0 (B6) 
[ITII} t>0, x=6: dCp/dx=0 “ (By) 


Converting equation (Bl) into the form of Laplace 


= iS 


transformation (Gay Laplace transformed function of Ca; 
Caos Sie ememchden, e)athe : 

sCq = Dad?Ca/dx? - k,“Ca 

d?Ca/dx* = (stk, “)/Dp (B8) 
Solving this second order differential equation, 

Cy (s) = pe’* + ge “* (w = V(stk, )/D, . 

Pp, gq: integration constants) 

From the boundary conditions (B5), (B6): 

p+q=Cy, /s, p= 0 
Consequently. 

Ca(s) = (Ca,/s)exp(-x/(s +k, )/Dag ) (B9) 

Doing the reversed Laplace transformation, the following 
analytical solution for the concentration function of oxygen 
Ca(x, t) is obtained: 

{Ca (x, te) Cay) = oa(x, t) 
(1/2) {exp (-xvk,7/Dy ) 

erfc(x/2vDat - vk, t) 

+ exp(xvk, 7/Dq ) 


- .erfc (x/2/Dat +i vky? tha (B10) 
and the flux function of oxygen Na(x, t) is: 


NA(x, t) = DadCa/dx 
= (Cao VR, 77DA) On + (Cyy/2) VDayTE 
[exp (-xVk, ~/Da) -exp{-(x/2VDqt —- vk,7t)*} + 
exp (xV/K,~/Dq) -exp{-(x/2/Dnt + VK,7t)7}] (B11) 
From equation (B2), the following Laplace transformed form 
is derived: 
: sc* = AC, €, expiO<C feyx i Asti k,c* - k, “Ca (B12) 
(C*, Laplace transformes function of C*; 
Ce = Pay eGas, Eydte) 
Inserting equation (B9) into (B12) 
C*(s) = AC ,€,exp(-C,€,x) /s(stk, *) (B13) 
- k, “Cy ,exp{-xV(stk, 7) /Dy /s (stk, ~) 
By reversed Laplace transformation 
{Ca*(x, t)/C)} @abl-exp(-k ft) F7K} 
“{Aevexp HC Saya yt k, ~(Ca,/C,) da lx, sd pal 
Inserting the equation (B9), (B13) into (B3), 


~ 24 ~ 


s7otiCa(k, t) /Cpe=——i1- exp (-k,t) }Ae, exp (-C,€,x) 
= Keer beng /Colon(x, tt) > exp (-kit) 
By integrating concerning times (t), 
(Ca lx, tt) /C)} = [tat le-expi-k oe) ]/k , "Ae ,exp(-C,€, x) 
eK, kg tGn ./ie bie bndt 
+ {I-expt-kt)t (Ca, /G.) bn (B14) 
Summarizing the analytical results as dimensionless 
concentration functions, the equations (11), (12), (13), and (14) 


described in section IV are obtained. 


APPENDIX C: Derivation of the parabolic law 
Let us derive the deterioration depth as a function of 
exposure time from the boundary condition [III] shown in 
equation (B7): 
9CB/9x| -_.= ielenasiy 4 eb fete - {1-exp(-k, t) }lexp(-C,c, x) 
Pk ok,) np Coe rgexiaeda tc, tidat }| 


+ l-exp(-k,t)} (Cao/C,*)-d¢a/dx|,_5=0 = (C1) 
We ‘can neglect the first term at comparatively short exposure 
pemommas yet expl— Kk, tC)... SOpkhart : 
2 us 
Oa Kad efoxl fbn (xy st) dt Fi) 


+ (fi-stexp (=k, t) }odn/ox Lm C2) 
2oa/2x| ,_¢ = vk,“/Dy [exp( 6vk,7/Da)erfc(6/2vVDat + vk, “t) 
- exp(-6vk | 7/Da) erfc(6/2vB,t - vk, 7t)] 
+ (1//mTD,t) [exp (6VKy */Da) exp{-(6/2VDat + vk,~t) *} 
+ exp(-é6vk,7/Daq) exp{-(6/2V/Dat - k,“t)?*}] 
(C3) 


Making use of the following approximation: 
erfc(6/2/Dat = Vk, 7t)= 1- (2/Vm) (6/2/Dat = Vk, 7t) 
exp{-(6/2/Dat = vk, 7E)= 1-(1/2) (8/2/Dat. * ve, 7) 
Let oxp( Oks tyeeekre 


doa/ox|,_5 * VKy 7 Dy {(6//Dat sinh (6yK,~7Da) 
4°2Vk)-t cosh ( §YK | ~7Dy) 
Hd 1/2 SDReMed2cosh'Syvk ;“7Da) 


(67/4Dat) cosh(é/K, 7Da) ; 


a 75 = 


- (6¥K,7/Da) sinh (6/k, ~/7Dq) 

k,t cosh(évk, 7/Da) 

~ Wy /mDg “(6 7Da) 7k, & + 27K, “t) 
+ (1/2/nDqgt ) {2-(6*/4Dat) - (k,~/D,) 6? - ky} 
= (4k,“t - 1) (6*/8DatYmDat)+ (3k t + 2) /4D 


+ (Sk Mer + 2)i/2¥ tae (C4) 
S52a/dx| pg = (4k, “t + 1) (6?/4D,/mDAE 
ey Ki, /mD_t Tt at) (Coo 
Using the relation: 
r te 
a/dxi Sq oal(x, t)dt}| 3 = I, (86, /ax| ,og) dt 


and inserting equation (C4) and (C5) into (C2), we can obtain 
the following equations: 

(kya ka 26Ft (4k Cees CEeaeIBDA Ce, tle Ly 

+e k ital GOR (I met 1) ot 4D ELS Ee en 0 (C6) 


6*" S"4Datikiy (ky te2ky = 2k) 2 2k 1a 
(ka * = 92k) cP aaeRE 2k 

= 4Dat(A + B/{4k,“(k, - 2k,“)t + 3k, -. 2k,7}] 
(C7) 
A =a(k, + 2k) 7) /(k + Seehkame (C8) 
B = 8k,*/(k, - k,*) (C9) 

AS, a, result: 

6 = VBt (C10) 
B =-Da (5k, 2k) hice = 2k) (hn ks “4S er 
(cl 


- 16 - 


FIGURE CAPTIONS 
Fig.l. Incidence of UV-light upon the surface of a 
polymeric material. 
Fig.2. Material fluxes within an infinitesimal layer of 
a polymeric plate. 
Fig.3. Influence of photooxidation on the distribution of 
oxygen. 
Fig.4. Time dependence of the distribution of oxygen. 
Fig.5. Influence of diffusion on the Oiestriupution of 
oxygen. 
Fig.6. Time dependence of the distribution of the activated 
polymer molecules. 
Fig.7. Profile of the rate of production of deteriorated 
polymer molecules. 
Fig.8. Time dependence of each concentration function at 
constant distances from the surface of a polymeric plate. 
Fig.9.a) Change of IR absorption with distance below the 
surface of polystyrene (PS) polymer exposed outdoors (at 
Choshi) for three years (from Watanabe et al. [10]). 
Fig.9.b) Change of IR absorption with distance below the 
surface of polyoxymethylene (POM) polymer exposed outdoors 
(at Choshi) for five years (from Watanabe et al. LLod) 
Fig.10.a) Hardening of a film of sodium silicate water-glass 
(Na,O-3.2Si0,) for the surface treatment of asbesto-cement 
board. 
Fig.10.b) Deterioration of films of sodium silicate water- 
glass for the surface treatment of a steal plate during 
immersion in 50°C hot water (drying at 120°C). 
Fig.l1l.a) Least-squares plot using the power law for the 
depth of the deteriorated layer of polyvnyl chloride (PVC) 
after various periods of outdoor exposure. 
Experimental data are from Kubota et al. [14]. 
Fig.ll.b) Least-squares plot using the power law for the 
depth of the deteriorated layer of various polymers after 
various periods of outdoor exposure (at Choshi). 


Experimental data are from Kubota et al. [14] 3: 
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FOREWORD 


This paper presents the experimental results of tests on two half- 
scale 3-story shear walls with connecting beams, which were conducted as 
support tests of a full-scale 7-story reinforced concrete building tested 
in 1980 under the "U.S.-Japan Cooperative Research Program Utilizing 


Large-Scale Testing Facilities". 


The load was applied to the shear walls by operating one horizontal - 
and two vertical jacks during the tests, in order to produce similar 
stresses in the shear wall as those of the 7-story structure subjected to 


waters! load. 


Detailed measurements were made of horizontal and vertical displace- 


ments and of strains of reinforcing bars. 


The elastic and inelastic behavior of the shear walls, which includes 
the effect of connecting beams, is discussed in detail based on test and 


analytical results. 


The shear walls yielded in flexure preceeding shear failure. Thus, 
these tests are the first in Japan for a flexural type multistory shear 


wall with large scale, 
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PLANER TESTS ON REINFORCED CONCRETE SHEAR 
WALL ASSEMBLIES 


— U.S.-Japan Cooperative Research Program — 


by 
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M. Yoshimura *2 
H. Lsoishi *3 
S. Nakata x4 


Building Research Institute, Ministry of Construction. 


Tsukuba, Ibaraki (Research Paper No.98 January 1983) 


ABSTRACT 


This study was done as support tests for the pseudo-dynamic test 
of a full-scale seven story reinforced concrete building scheduled 
for 1980 under the "U.S. -Japan Cooperative Research Program Utilizing 
Large Scale Testing Facilities". 

Tests were made on the shear wall-beam assemblies of the lower three 
stories of the internal structural frame in the ridge direction of the full- 
scale seven story building. Test specimens were of two kinds. One was a 
half-scale test specimen (Wl) with reinforcements corresponding to the 
full-scale design building. The other (W2) was of the same cross-section 
and shape as Wl but with increased reinforcements and increased strength 
for the sake of comparison, since the full-scale seven story building has 
a low base shear at ultimate strength. The boundary beams at each story 
of the test specimens were connected to each other by tie bars with both 
pinned ends to prevent vertical displacement, so that bending moments of 
the boundary beam would decrease wall moment. 

The moment distribution applied to the specimen was similar to that 
of the wall in the seven story buildings under the external forces of 
inverted triangular distribution. The controlled moment together with 
the axial force was applied to the test specimen at the top of the third 
story, by adjusting the load of two jacks installed vertically on each 
side of the top of the test specimen... Load hysteresis consisted mainly 
of repetitions to displacement values where considerable plasticization 
was expected from a preliminary seismic response analysis of the ene 


scale seven story building. 


Detailed measurement was made of the horizontal and vertical displace- 


ments at the left and right ends of the respective stories, the stretching 


forces of the boundary beams in the hinge region, strains in the main 
and shear reinforcements of the beams and columns, and strains in the 
vertical and horizontal wall reinforcement. 

The test results were as follows: 
i) In both Wl and wW2, the load-deflection hysteresis displayed bending- 
type behavior with a decreasing loss of strength from cycle to cycle at 
the same displacement level. Final failure was a collapse of the column 
concrete on the compression side accompanied by a sharp decrease in strength. 
The rotation angle of the shear wall at the time of rupture was 1/50 for 
W2 and 1/40 for Wl, but in the case of Wl, a stable load-deflection hyster- 
esis was maintained only to 1/50. 
ii) The respective stories produced a nearly identical rotation angle. 
Particularly in Wl, the deformation was mainly, determined by the rotation 
due to the elongation and shrinkage at the base of the shear wall. 
iii) In both Wl and w2, all the vertical wall reinforcement at the bottom 
of the wall and the main reinforcements of the tensile side columns yielded 
before the rotation angle reached to 1/200. The main reinforcement of the 
boundary beam and the adjacent slab reinforcement also yielded before the 
rotation angle reached 1/200. 
iv). The test specimens had no beams in the wall, so the shear cracking in 
the upper and lower stories penetrated the slab and produced continuous 
cracks throughout the upper and lower stories. However, the absence of 
beams had little influence on the rigidity and strength. 
v) The shear strength of the test specimens was considerably higher than 
the maximum flexural strength, so in the pseudo-dynamic test of the full- 


scale seven story building, the shear wall is not expected to fail in shear. 
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1. PURPOSE 


This study was carried out as support tests for the pseudo-dynamic test 


of the full-scale seven story reinforced concrete building, scheduled for 
1980 under the "U.S.-Japan Cooperative Research Program Utilizing Large 
Scale Testing Facilities". 

This full-scale seven story building has three spans in the ridge 
direction and two spans in the transverse direction, and has a shear wall 
through out all stories at the center of the internal structural plane in 
the ridge direction (See Figs. 1.1 and 1.2). In such a building, it is 
generally predicted that the properties of the shear wall will greatly 
influence the aseismic performance of the building. In fact, in the 
preliminary analysis of this structure, it was pointed out that the shear 
wall would carry about 85% of the total horizontal force acting on the first 
story in the elastic region and bear 60% of the maximum flexural strength 
of the building. The details of this building were decided basing on 
discussions between the U.S. and Japan from the structural planning stage. 
Resultantly, it incorporates various points can not be seen in our country, 
such as omission of the beams in the shear wall. This full-scale seven 
story building has a base shear coefficient of about 0.24 (except for the 
live load as weight) at the maximum flexural strength and is thus designed 
to resist strong ground motions with ductile behavior. Consequently, the 
shear wall must have a considerable plastic deformation capacity. 

The most important focus of thses tests is on the shear wall and the 
beams of the lower three stories of the full-scale seven story test building. 
For this objective, two specimens were made in half- scale. Specimen Wl has 
reinforcement corresponding to that of the lower three stories of the full- 
scale seven story building. Specimen W2 has the same cross-section and shape 
as Wl, but it has increased reinforcement and thus a greater strength. 


The main objectives of this study are listed below. 


i) To obtain load-deflection curves for preliminary analyses of the 
full-scale seven story building and to analyze the test results, or, 
more particularly, to confirm the deformation characteristics and the 
types of failure. 

ii) To make clear the vartical displacement and rotation of the shear 
wall and the effects of the boundary beam and to examine the effects of the 
transverse beam in the full-scale seven story building. 

ii) To examine what effect the absence of the beams in the shear wall 
is likely to have, and note any problems with the reinforcing method, 
and more particularly, in the anchoring method of main reinforcing bars 
in the boundary beams. 

It should be noted that as part of the U.S.-Japan Cooperative 
Research Program, similar test were carried out on the beam—column 


assemblages, of the full-scale seven story test building. 


2. TEST "SPECIMENS 


2.1 Description of Test Specimens 

The focal point of these tests is on the shear wall-beam assemblies 
of the lower three stories of the sull-scale seven story test building. 

There were two specimens. One was a one-half scale specimen (Wl) with 
reinforcement corresponding to that of the full-scale building. The other 
specimen (W2) was of the same cross-section and dimensions as Wl but with 
increased reinforcement to increase the strength, because the full-scale 
seven story building has a low base shear at reserve strength. 

The length of the boundary beams connected with the shear wall, was 
one-half of the span, taking the center of the beam of the full-scale build- 
ing as the inflection point. Further, considering that the behavior of the 
beams and shear wall, and hence the test results, would be greatly affected 
by the presence or absence of floor slabs, particularly a shear wall without 
beams as in the case of the test specimens, floor slabs were provided at 
each floor with a width as specified in the Building Code Requirements for 
Reinforced Concrete, Architectural Institute of Japan. This width was ef- 
fectively equal to that specified in the U.S.A. For accurate evaluation 
of behavior on the beam-column joint, a cross beam was provided of a length 
equal to the slab width. The beam provided at the top for applying the force, 
and the foundation beam for taking the reaction were made as rigid as possible. 
The both top and bottom main reinforcements of the boundary beams were bent 
down at the inside of the column. The wall reinforcement and main column 
reinforcement were lapped at the floors of the second and third stories, 
respectively. The concrete covering for reinforcement was one-half of the 
specified value because the test specimens were of one-half scale. Concrete 
placing was done at the floor surfaces of the foundation, first story, second 


story and third story. 


Drawings of the reinforcement arrangements and a table of the rein- 
forcement for Wl and W2 are shown in Figs. 2.1.1, 2.1.2, respectively and 
also in Table 2.1.1. The difference in the reinforcement arrangements bet- 
ween Wl and W2 was mainly the amount of main reinforcement for the column, 
the amount in W2 being 2.5 times that of Wl. As for the other reinforce- 
ments, i.e. the wall reinforcement, and main and the shear reinforcement 
of the boundary beam, the arrangement of W2 was from 1.5 to 2.0 times that 
of Wl. 

The schedule of work and the testing of the specimens is shown in 
Tab ke? al. 

2.2 Mechanical Properties of the Materials 

The concrete was standard, using river sand and river gravel with a 
design strength of Fo = 270kg/cm?. Three test cylinders were prepared for 
both the compression and splitting tests from the concrete for the founda- 
tion, first story, second story and third story. They were cured indoors 
until the day of testing under the same conditions as those of the test 
specimens. Cylinder sizes were 10cm$ x 20cm for the compression tests and 
15cm x 30cm for the splitting tests. As for the reinforcement, deformed 
bars (SD35) De, D1O and D13 were used for the main portions of the test 
specimens and D22 for the main reinforcement of the foundation and loading 
beam. The test results for the reinforcement and concrete are shown in 
Tables 2.2.1 and 2.2.2. In Table 2.2.1, the values for the reinforcement 
are the mean values of three tests of the respective sizes. Stress-strain 


curves for the reinforcements are shown in Fig. 2.2.1. 


3. EXPERIMENTAL PROCEDURE 


Bal, Application of Force 

An inverted triangular distribution was assumed for the distribution 
of external force on the shear wall of the full-scale seven story building. 
According to the results of preliminary seismic response analysis of the 
full-scale seven story building (using as an input seismic wave the NS com- 
ponent at the first floor of the Construction Engineering Department Build- 
ing at Tohoku University during the Miyagi-Ken Oki Earthquake in 1978, with 
the maximum acceleration set at 350cm/sec?, the ratio between the moment at 
the top of the third story to that at the base of the wall varies greatly 
through the response time history. This ratio was approximately within the 
range from 1/1.5 to 1/4 when yielding occurred at the base of the shear wall. 
On the other hand, assuming an inverted triangular distribution of external 
force and a yield mechanism of the boundary and transverse direction beams, 
the ratio of the moment in the shear wall at the top of the third story to 
thateat the base of the wall is 1/2.21. 

The bending moment diagram when the inverse triangularly distributed 
external force acts on the one-half scale seven story shear wall assemblies 
is shown in Fig. 3.1.1. The dashed line represents the bending moment when 
the external force is assumed to be of a uniform distribution. 

The loading apparatus is shown in Figs. 3.1.2 and 3.1.3. The horizontal 
load was applied as a concentrated load at the top. The boundary beams were 
connected to each other by tie bars with both pinned ends to constrain their 
vertical displacement. Two center hole jacks were used to apply the vertical 
load on the left and right sides of the force transfer beam at the top of the 


third story. This system allowed the left-and right-vertical loads (Nz, and 


Nr) to vary along with the horizontal load, according to the formulae (1) - (4) 


shown below so that: 

a. the axial force in the first story wall would correspond to the axial 
force of the full-scale seven story shear wall assembly (d09 = 264,000/ 
(20 x 500) = 26.4kg/cm), and 

b. the moment distribution would correspond to that occurring when an 
inverted triangularly distributed external force was applied to the full- 
scale seven story shear wall assembly. 

In calculating this moment, it was assumed that the boundary and trans- 
verse beams would yield when the base of the shear wall yielded. Before 
yielding, their loads were assumed to increase linearly with the horizontal 
load and the effects of their bending moments were taken into account. 

Vertical loads of Test Specimen W1 : 

Ny w=" 24, 6%) -"05356P (€) , (1) 
Nee =e2o.6* te SOP (t) (zy 

Vertical loads of Test Specimen W2 : 

Nr = 24.6* - 0.364P(t) (3) 

NR = 25.6* + 0.364P(t) (4) 

(* The difference of 1 ton between Ny, and Np is aué’to the 
weight of the loading apparatus) 

The bending moment diagrams for the test specimens, obtained according 
to the foregoing loading method, are shown in Fig. 3.1.4. The bending moment 
diagrams at maximum flexural strength are shown in Fig. 3.1.5 together with 
those of the half-scale seven story shear walls. Since they agree with each 
other, this loading method seems well suited to represent the moment distri- 
bution of the lower three stories of the seven story shear wall assemblies. 
Because of limitations in the loading apparatus, for specimen W2, the control 
of the vertical load was possible only up to 58 tons, which was just short 


of the yield load. 


3.2 Data Measurement 

The measurement points for strain and displacement, are shown in Figs. 
3.2.1 and 3.2.2. The items measured are given below. Figures in parentheses 
show the gauge number in Figs. 3.2.1 and 3.2.2. 

i) Strains in Concrete, steel reinforcement and tie bars were measured with 
electric resistance wire strain gauges. 
a) Strains in the main reinforcing bars of the column, beam, and longi- 
tudinal wall reinforcement (1 - 48). 
b) Strains in the column hoops, beam stirrups and horizontal wall rein- 
forcement (49-68). 
c) Strains in the slab reinforcement (69-76). 
dad) Strains on the concrete surface (77-109). 
e) Strains in the tie bars connecting boundary beams (140-157). 
ii). Horizontal and vertical displacements were measured by inductance type 
displacement transducers (110-139). 
iii) Jack loads were measured by load cells (158, 148 and 149). 

The marks A and |m represent the displacement transducers provided for 
measuring rotation angles in the hinge regions at the beam and column ends, 
respectively. They were positioned at a distance equivalent to their depth, 
25cm away from the beam and column ends and the rotation over this distance 
was measured directly. The measuring jigs for the horizontal and vertical 
displacements were attached to the transverse beams. The measurement point 


for vertical displacement was located 10cm away from the column surface. 


4. TEST RESULTS 


4.1 Loading History 

Figs. 4.1.1 and 4.1.2 show the loading histories recorded for Wl and 
W2. Tables 4.1.1 and 4.1.2 show the maximum loads and the values of the 
rotation angle due to horizontal displacement (mean valué of the horizontal 
displacements at the loading side and the opposite side of loading of the 
shear wall) of the respective stories for each load cycle. Load histeresis 
consisted mainly of repetitions to displacement values where considerable 
plasticization was expected from a preliminary seismic response analysis of 
the full-scale seven story building. Neither Wl nor W2 had any appreciable 
decrease in the strength during repetitions of the same displacement amplitude 
before their collapse. 
4.2. Crack Patterns and Modes of Failure 

Flexural cracks were observed at the base of the column for specimen 
Wl in the first story at a horizontal load of 0 = 12 tons (T = O/tk = 4.8kg/cm?) 
and 14 tons (T = 5.6kg/cm2) in the positive and negative loading respectively. 
W2 had a hair crack at the base of column at a positive load of 10 tons (T = 
Akg/cm*) . The crack became clearly observable at the load of 20 tons (T = 
8kg/cm7) . At negative loading, the crack appeared at the load of 20 tons. 
Thereafter, flexural cracks corresponding to the hoop spacing was observed 
gradually in the upper part of the column. In Wl, a flexural crack produced 
at a position of about 1.5D from the base of the column developed into a 
shear crack in the wall, under both positive and negative loads of 24 tons 
(t = 9.6kg/cm2) . The angles of these cracks were nearly horizontal, probably 
due to the influence of bending stress, so the cracks should correctly be 
called flexural shear cracks. Wl also had a slight shear crack at the bottom 
of the column on the tensile side under a negative load of 20 tons. In W2, 


a similar shear crack in the wall was produced under a positive load of 16 


tons (T = 6.4kg/cm?) and a negative load of 20 tons(T = 8kg/cm7) . Increased 
load was accompanied by shear cracking in the wall, while flexural cracks 
developed upward. A flexural crack occurred in the column at the second 
story level under a positive load of 28 tons and a negative load of 32 tons 
for Wl, and under a positive and negative load of 32 tons for W2. On the 
second and third stories, shear cracks in the wall were produced almost si- 
multaneously under a load of 32 to 36 toms (t = 12.8%14.4kg/cm?) . There- 
after, a flexural crack developed in the column at the third story level. 
The shear cracks in the walls of the second and third stories were of an 
angle of about 30°, smaller than the 45° angle of cracks produced by pure 
shear, but greater than the angle of the shear cracks in the first story. 
The foregoing was the phenomena observed before reaching a rotation 


angle of 1/400 of the shear wall, as defined by the following formula: 


R = 6p/h = 6p/4750mm 

where R = rotation angle of the shear wall; 
Sp = horizgntal displacement of the top of the third story, and 
h = height from the top of the third story to column base of the 


first story. 

In the following, the term "rotation angle" means the rotation angle 
of the shear wall unless otherwise specified. 

Up to this Bo irttts all of the boundary beams had flexural cracks. After 
a rotation angle of 1/400, these cracks increased gradually. At a rotation 
angle of 1/100, clear shear cracks were observed in the first story column 
of both Wl and W2 and crushing of the concrete was observed at the base of 
the column and at the compressive bottom,sides of the boundary beams. In 
aa the concrete of the first story column began to break up eae 
due to shear compression at a rotation angle of 1/50, and with three repeti- 
tions under 1/50, fracture of the longitudinal wall reinforcement took place. 


Thereafter, at a rotation angle of 1/40, the majority of the vertical wall 


reinforcement and one main reinforcement of the column fractured, and during 
unloading, the concrete at the base of the column failed in compression. In 
W2, on the second repetition at a rotation angle of 1/50, the concrete at the 
base of column failed in compression as in Wl. The fact that the rotation 
angle of W2 at the time of the collapse was smaller than that of Wl may Be 
explained by the higher axial compressive stress of the column on the com- 
pression side because of the larger maximum flexural strength. 

The test specimens had no beams in the wall so that shear cracks produced 
in the upper and lower stories penetrated through the slabs to form continuous 
cracks throughout the upper and lower stories. However, these cracks were of 
small width and are not considered to have had appreciable effect on the be- 
havior of the shear wall. 

Tables 4.2.1 and 4.2.2 list the crack patterns, and Figs. 4.2.1 to 4.2.9 
show the cracks after reaching specific rotation angles. 

4.3 Strains in the Reinforcement and Concrete 

Figures 4.3.1 and 4.3.2 show the yielding patterns of the reinforcement 
in Wl and W2. The main reinforcements at the base of the column and most of 
the vertical wall reinforcement yielded before a rotation angle of 1/400 and 
completely yielded by 1/200. Judging from the gauges aftached at an interval 
of 15cm to the anchored portion of the main column reinforcement, the rein- 
forcement yielded at the 15cm from the surface of foundation for Wl and the 
30cm for W2, and it seems that a considerable slippage of the reinforcement 
occurred at ultimate stage. As for the boundary beam, almost all the rein- 
forcement, including the slab reinforcement, yielded before a rotation angle 
of about 1/4200. 

As for the shear reinforcement, those at the base of the column yielded 
after a rotation angle of 1/200. Those at the upper part of the first story 
column and the bottom of the second story column did not yield and had a 


strain value of less then 1000u. The shear reinforcement at the beam-column 
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joint on the second floor also did not yield and had a strain value of 
more than 1000”. Concerning horizontal wall reinforcement, it was observed 
that in Wl only the horizontal reinforcement at the lower part of the wall 
yielded at a rotation angle of 1/200. In W2, however, horizontal reinforce- 
ment on the first story began to yield after a rotation angle of 1/200 and, 
before a rotation angle of 1/100 was reached, almost all of the horizontal 
reinforcements on the first story had yielded. While the structure had no 
beam in the wall, yielding occurred in some of the vertical and horizontal 
reinforcement near the corner of the beam-column joint. Summing up the 
yielding patterns of the reinforcement in both Wl and W2, the main reinforce- 
ment of the base of the column and the vertical reinforcement at the base 
of the wall yielded before a rotation angle of 1/400, the main reinforcement 
of the boundary beam and the slab reinforcements yielded before a rotation 
angle of 1/200, and finally, the horizontal wall reinforcements and the hoop 
SE the column only yielded at or beyond the maximum flexural strength. 

Figs. 4.3.3 through 4.3.8 show the shearing stress of the wall plate 
as a function of the horizontal load obtained by wire strain gauges on the 
concrete surface shown in Fig. 3.2.2. Distribution of stress in the same 
horizontal plane are shown in Figs. 4.3.9 and 4.3.10 for the respective loads. 
The dashed line in the drawings shows the stress distribution obtained accord- 
ing to the I-beam primary elemental theory of bending. At the central plane 
of the wall plate on the first story, the experimental values are in agreement 
with the values calculated for both central and end parts. For the shearing 
stress beneath the slabs on the second and third stories, the experimental 
value at the central part of the wall plate is also in good agreement with 
the value calculated. However, at both ends, the experimental values are 
considerably smaller than those calculated. Further, the shearing stress at 


the loading side varies little with an increasing horizontal load. The value 


of the shearing force obtained by integrating the shearing stress distributed 

on the horizontal plane of the wall agrees relatively well with the value cal- 
culated for the central plane of the first story, but for the plane beneath 
the slab, the value is considerably smaller than the value calculated. Com- 

paring the shearing force acting on the column beneath the slab (obtainable 

by deducting the shearing force acting on the wall plate from the horizontal 
shearing force) with the shearing force acting on the column at the central 

plane of the first story, it may be said that a greater shearing force acts 

on the column beneath the slab. 

4.4 Elastoplastic Behavior 

Figs. 4.4.1 through 4.4.6 show the relation of horizontal load Q to 
absolute horizontal displacement § of the respective stories of Wl and W2. 
Figs. 4.4.7 through 4.4.10 show the relation of the horizontal load Q to the 
story displacement 6p of the second and third stories. Both Wl and W2 give 
a stable spindle-shaped loop and have no notable decrease in the strength 
for load repetitions within a rotation angle of less than 1/100. Further, 

Wl and W2 have substantially the same rotation angles at the first, second 
and third stories. 

Figs. 4.4.11 and 4.4.12 show the relation of top vertical displacements 
toy of the column of the third and first stories to horizontal load Q at the 
loading side for Wl, and the similar relation of vertical displacement ROV 
to horizontal load Q at the opposite side of loading for W2. At the top on 
the tensile side of the third story, both Wl and W2 gave values of a stretch 
of Scm or more for a rotation angle of less than 1/50. However, on the com- 
pression side, the longitudinal shrinkage was less than lcm. For both Wl 
and W2, the vertical displacement of the second floor is not very different 
to that of the top of the third story, so it may be said that the greater 


part of the vertical displacement occurs on the first story. 


Skeleton-curves for the horizontal and vertical displacements of the 
respective stories up to a rotation angle of 1/100 are shown together with 
the analytical results in Figs. 4.4.13 and 4.4.14. For the stress-strain 
relationship of concrete, the e-function proposed by Dr. Umemura was used 
and a bilinear model with the second gradient taken as 1/200 of the elastic 
rigidity was used for the stress-strain relationship of the reinforcement. 
Slippage of the reinforcement at the foundation was not taken into consider- 
ation in analysis. The test results for both horizontal and vertical dis- 
placements agree well with the analytical results. It may thus be said 
that Wl produced only flexural deformation and that shearing deformation 
was negligible. However, in the case of W2, shearing behavior was noted 
in the first story. 

Figs. 4.4.15 through 4.4.20 show the chord rotation 6; obtained from 
the vertical displacements of the top, third and second floors. 

81 = (19vi-1 7 RSvi -1)/2y 
where 8; = angle of rotation of the floor of the i-th story; 

L°vi - 1 = vertical displacement at the top of the left side column of 

ene i" 2) = SCOry: 

ROvi-1 = vertical displacement at the top of the right side column 

ofythe* (i - Ly*story; and 
&£y = distance between the measurements of vertical displacements. 

In the drawing, the dashed line indicates that proper measurements 
were not obtained because of an insufficient stroke of the displacement 
transducer. Both Wl and W2 had angles of rotation of substantially the 
same values and rotation angles due to the horizontal displacements were 
previously shown in Figs. 4.4.1 through 4.4.6. Figs. 4.4.21 through 4.4.24 
show the values of the net horizontal displacement Spr defined by the follow- 


ing formula, up to a rotation angle of 1/100 for the third and second stories. 


SpintieSnaeqe faa 
where Sp; = net horizontal displacement of i-th story; 

Spy = story displacement of i-th story; 

6; = angle of rotation of the i-th floor; and 

hj = story height of i-th story. 
The net displacement was small enough, and it seems that the second and 
third stories displayed a rigid behavior. 

From the foregoing, it can be seen that deformation of the shear wall 
of Wl showed bending yielding type behavior governed mainly by the rotation 
due to elongation and shrinkage of the base of the first. stony.....,.1dn, the, case 
of W2, it may be said that it showed similar behavior to Wl, although a 
kind of shearing deformation was produced in the first story. In the first 
story portion of W2, the so called "anisotropic plate characteristic" was 
seen, where a vertical displacement occurs as a shearing force acts on a 
wall plate with shear cracks. 

4.5 Elastoplastic Behavior of Boundary Beams 

Figs. 4.5.1 through 4.5.8 show the relation of shearing force to hori- 
zontal load of the boundary beams, calculated in accordance with the follow- 
ing formulae from strains of the wire strain gauges attached to the tie bars 


connecting boundary beams. 


p2] = AsEg(€z - €2). = 7-89 x 210ton x (eq E92) 
= L657(€4) TaE2) ton 
and BQ2 = AsEseg = 7-89 x 210ton X Ep = 1657e5 ton 
where B21, B22 = shearing forces of the beams of the first and second 


stories, respectively; 
Ag = cross-sectional area of the tie bar connecting the 
boundary beams; 
Eg = Young's modulus of,,the tie bar; and 
Ej, €g = strains of tie bars at the first and second story, 


respectively. 
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The relation of assumed shearing force to horizontal load of the boundary 
beam for control of the vertical loading jacks provided at the top of the 
respective test specimens is shown by the broken line in the above mentioned 
figures. The shearing force produced in the beam while the top side of the 
beam is under tension is small until the wall begins to yield, but increases 
sharply after yielding of the wall. The forces measured when tension existed 
on the bottom side of the beam and during unloading when tension existed on 
the top side of the beam agrees well with the assumed values. Thus it may 
be said that the test specimens have developed stresses approximately equal 
to those assumed (See Fig. 3.1.4). 
Figs. 4.5.9 through 4.5.16 show the relations between the shearing 

force of the boundary beam stated above and the rotation angle of the boundary 
beam calculated from the vertical displacements of the left and right columns 
of the shear wall and angles of rotation 8 as described in paragraph 4.4. 

RBOs J? RBOMn Fe Rods #992, Sp yi. +.237584. 

LBOi = LSvi ~ 912B = pSyi ~ 137504 
Here, 

RBOir LBOi = displacements of the left and right boundary beams 

on the i-th floor; 
Rov, Lov = left and right vertical displacements on the columns 
of i-th floor; 


6; = angle of rotation o£ the i-th, floor: and 


Lp clear length of the beam. 

Except for the deviation of the loading axis because of an initial stress 
produced in the beam while attaching the tie bars, the four left and right 
beams display nearly similar loops for both Wl and W2. The beams are sub- 
ject to yield at a rotation angle of 1/200 and to eae nei on there- 


after. During tension exist on the bottom side of the beam, rotation angles 


approximately equal to those of the shear wall are produced, and during 
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tension exist on the top side of the beam, a rotation angle about two times 
that of the shear wall is observed. This characteristic presents a remark- 
able pinching phenomenon as can be seen in the results of the member test 
on the beam-column assemblages carried out separately to this study. 

The test specimens included no transverse beams. However, considering 
the amount of elongation and shrinkage of the columns on the tensile and 
compressive sides, it seems sufficient to use the shearing force corresponding 
to the maximum flexural strength of the beam for calculating the restaining 
force by the transverse beam located at the tensile column side at maximum 
flexural strength of the shear wall, and to use the smaller value of the 
maximum flexural strength of the beam, for bending moment calculations for 
the transverse beam on the conpressive side of the shear wall. 

4.6 Strengths 

The values for cracking strength and ultimate strength of the test 
specimens are shown in Table 4.6.1. The analytical results in the table 
are based on the following formulae. 


i) Flexural crack moment 


Maa Fee ae Nae AS 
QBc = Mpc/ (M/Q) 


ii) Shearing force at flexural shear cracking 

Opse = 0-265 - t - WYFG + Mpc/(M/Q - 0.5D) 
iii) Shearing force at shear cracking 

gee Balke” 1) (Bf. eae aR ore 


" 1 + 1.05t/B 
CBee y £1177 78 


cok'= 2) are = 0, 00S46r. eee 


iv) Shearing force at maximum flexural strength 


0.9atsoy D + 0.4 wsoy D + 0.5ND{1 - N/(BDF,) } 


Mem 


Ory = Mpy/ (M/Q) 


v) Shearing strength 


O23 
Qsm = AREER Pe rae + 2.7V sowyPwe + 0.1loo} bej 
M/(Qq) + 0.12 


With respect to the flexural cracking strength, the test results agree 
well with the analytical results. According to the analysis, both Wl and 
W2 reach yield strength without producing shear cracks, but according to 
the test results, shear cracks are produced in all stories. The reason 
for this difference the test specimens were subject to the influence of 
bending stress as can be seen from the fact that the angle of the crack is 
considerably smaller than 45°. According to the tests, these cracks are 
to be called flexural shear cracks. With respect to the maximum flexural 
strength, the values agree well with those given by the e-function method. 
The values obtained by the approximate formula were somewhat lower. 

The analytical shearing strength is far beyond the maximum flexural 
strength in the case of Wl. Further, judging from the crack pattern gener- 
ated in Wl, the shear strength seems to be enoughly stronger than compaired 
with the flexural strength. It may thus be considered that in the pseudo- 
dynamic tests of the full-scale seven story building, the shear wall is 


unlikely to be subject to damage by shear. 


5... CONCLUSIONS 


Tests on the shear wall portion of the lower three stories of a seven 
story building which contains a shear wall were conducted, and the following 
results were obtained. 

i) In both test specimens, Wl and W2, the load-deflection hysteresis of the 
shear walls showed a behavior of the flexural type with a small loss of 
strength from cycle to cycle at the same displacement. However, a rapid 
collapse resulted because the final failure was caused by compressive crush- 
ing of the concrete at the base of the column surrounding to the wall. The 
rotation angle of the shear wall at collapse was 1/40 in the case of Wl and 
1/50 in the case of W2. However, in Wl, a stable load-deflection hysteresis 
was maintained only until a rotation angle of 1/50. 

ii) The respective stories produced approximately aes) rotation angles. 

In Wl, particularly, the deformation was governed largely by. the elongation 
and shrinkage at the base of the shear wall. Also, in Wl, the values of 

the horizontal and vertical displacements agree well with the analytical 
values by the e-function method by Dr. Umemura, so that it is considered 
that the shearing deformation has little effect. 

iii) Vertical reinforcements at the base of the wall yielded entirely at 

a rotation angle of 1/200 in both Wl and W2. As for the hoops of the column 
and the horizontal reinforcements of the wall, at a rotation angle of 1/200, 
the hoops of the column footing and horizontal reinforcement at the wall 
footing began to yield for Wl. For W2, some of the hoops of the column 
footing and horizontal reinforcements at the first story wall began to yield 
at the rotation angle of 1/200, and the majority of these reinforcements 
yielded after a rotation angle of 1/100. 

iv) The boundary beams produced a greater rotation angle than that of the 


shear wall and yielded generally at the shear wall rotation angle of 1/200. 


Slab reinforcements on the second floor yielded almost simultaneously 

with the main reinforcements of the beam. The beams on the second and 

third floors produced substantially similar deformations. 

v) The measured maximum flexural strength agreed well with the analytical 
result by the e-function method. The calculated shear strength of the res- 
pective test specimens is enoughly stronger than flexural strength. Thus, 
it is not expected that the shear wall will be subject to shear failure in 
the pseudo-dynamic tests of the full-scale seven story building. 

vi) As the shear wall had no internal beams, the shear cracks in the upper 
and lower stories penetrated through the slabs to cause continuous cracks 
throughout the upper and lower stories. However, absence of the beams had 
little effect on the rigidity and strength of the wall. 

vii) There were no transverse beams in the test specimens. However, from 
the amount of elongation and shrinkage of the columns at the tensile and 
compression sides, it seems sufficient for calculating the restraining force 
by the transverse beam at the maximum flexural strength of the shear wall, 
to use the shear force at yielding of the beam for the tensile column side, 
and a value of the shear force somewhat discounted for the compression column 


side. 
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APPENDIX A: 

Results of Measurement of Reinforcement and Concrete Strains, Displace- 
ments and Loads of Wl. 

The following. load-strain curves and load-displacement curves of the 
reinforcement and concrete of Wl correspond to Figs. 4.3.1 and 4.3.2 in the 


text. 


APPENDIX B: 

Results of Measurement of Reinforcement and Concrete Strains, Displace- 
ments and Loads of W2. 

The load-strain curves and load-displacement curves of the reinforce- 
| ments and concrete of W2 shown in Figs. 4.3.1 and 4.3.2 in the paper are 


shown. 
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Fig. 1.2. (Continued) Full-Scale Structure 
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‘Fig.2.12 Definitions and Reinforcement of Specimen W2 


Fig.2.11 Definitions and Reinforcement of Specimen W1 
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Table 2.1.2 Schedule of Construction Model Work and Tests 


aetens 790 
Concrete placing ifor the foundation. 


Feb.15 | Molding work for the first story. 
Tensile tests of reinforcement to be used. 


Feb. 23 Molding work for the second story. 
Concrete placing for the second story. 


Feb.29 Reinforcement arrangements for the loading part and 
third story. 
Molding work for the loading part, third story. 


Apr.21-28 Wl loading tests. 
Compression and splitting tests of test cylinders. 


Reinforcement arrangements for the foundation. 


Molding Work for the foundation. 


Table 2.2.1 Mechanical Behavior of Reinforcing Bar 
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Pig.2.2.1: Stress Strain Relationship of Reinforcing Bars 
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Fig.3.1.1 Bending Moment(t.m) at Yielding of Half-scale Seven-story Shear Walls 
under Inverted Triangle Lateral Force Distribution, Where Dushed Line 
Shows the Bending Moment When the Lateral Force is Assumed to Be 
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Fig.3,1,3 Detall of Test Set-up 
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Fig.3.1.4 Bending Moment of the Specimens 
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Fig.3.1.5 Bending Moment at Yielding of the Specimens (Unit:t.m) 


* (Dusted lines express these of the half-scale seven-story shear walls.) 
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Fig.4.1.1 Loading history of Wl 
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Fig.4.3.3 Laceral load versus shearing stress relationship of wall of Wl (No.1) 
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Fig.4.3.5 Lateral load versus shearing stress relationship of wall of Wl (No.3) 
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Fig.4.3.6 Lateral load versus shearing stress relationship of wal] of W2 (No.1) 
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Pig.4.3.8 Laceral load versus shearing stress relationship of wall of W2 (No.3) 
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Fig. 4.3.9 Shearing stress distribution of central and upper 
part of the lst wall and upper part of the 2nd wall of Wl 


note : x=distance from center of wall 
( a=horizontal half clear length of wall 
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Fig. 4.3.10 Shearing stress distribution of central and upper part 
of the lst wall and upper part of the 2nd wall of W2 


note : x=distance from center of wall 
a=horizontal half clear length of wall 


Fig.4.4.1 Lateral load versus top floor deflection relationship of Wl 
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Fig. 4.4.2 Lateral load versus 3rd floor deflection 
relationship of Wl 
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Fig. 4.4.3 Lateral load versus 2nd floor deflection 
relationship on Wl 
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Fig.4.4.4 Lateral load versus top floor deflection relationship of W2 
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Fug.4.4.5 Lateral load versus 3nd floor deflection relationship 
of W2 
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Fig.4.4.7 Lateral lqad versus relati 
deflection relationship of 
3rd story of Wl 


Fig.4.4.6 Lateral load versus 2rd floor deflection 
relationship of W2 
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Fig.4.4.8 Lateral load versus relative deflection Fig.4.4.9 Lateral load versus relative 
relationship of 2nd story of Wl deflection relationship of 
3rd story of W2 


Fig.4.4.10 Lateral load versus relative deflection relationship of 2nd story 
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Fig.4.4.13,Comparison of analytical displacements with experimenta] results of Wl 
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Fig.4.4.15 Lateral load versus chord rotation relationship of 
2nd floor of Wl obtained by vertical displacements 
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Fig.4.4.16 Lateral load versus chord rotation relationship of 3rd floor 
of Wl obtained by vertical displacements ,;5,, and ,6,, ‘* 
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Fig.4.4.17 Laterl load versus chord rotation relationship of top floor 
of Wl obtained by vertical displacements ;5,, and pdy3 
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Fig.4.4.18 Laterl load versus chord rotation relationship of 
2nd floor of W2 obtained by vertical displacements 
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Fig.4.4.19 Lateral load versus chord rotation relationship of 3rd 
floor of W2 obtained by vertical displacemsts 
Love and Rov2 
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Fig.4.4.20 Lateral load versus chord rotation relationship of 
top floor of W2 obtained by vertical displacements 
toys and bys 
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Fig.4.4.21 Lateral load versus tangential 
drift relationship of 2nd story 
of Wl 
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Fig.4.4.23 Lateral load versus tangential 
drift relationship of 2nd story 
of W2 
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Fig.4.4.22 Laterl load versus rotation 


angle relatiohship of tho floor 
of W2 obtained by vertical 
displacements ;6,, and péy3 
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Fig.4.5.1 Lateral load versus shear force relationship VilgudeS.2)) Gaterall loadiversus (shear force relationship 
of boundary beam of 2nd floor at loading side of Wl of boundary beam of 2rd floor at opposite 


side of loading of Wl 
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Fig.4.5.3 Lateral load versus shear force relationship 
of boundary beam of 3rd floor at loading 
side of Wl 
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Fig.4.5.4 Lateral load versus shear force relationship of boundary beam of 3rd floor at 


, opposite side of loading of Wl 
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Fig.4.5.5. Lateral load versus shear force relatiohship 
of boundary beam of 2nd floor at loading side 
of w2 


Fig.4.5.6 Lateral load versus shear force relationship 
of boundary beam of 2nd floor et opposite side 
of loading of W2 
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Fig.4.5.7 Lateral load versus shear force relationship of boundary beam 
of 3rd floor at loading side of W2 
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Fig,4.5.8 Lateral load versus shear force relationship of boundary 
beam of 3rd floor at opposite side of loading of w2 
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Pig.4.5.9 Shear force versus deflection relationship of 
boundary beam of 2nd floor at loading side of Wl 
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Pig.4.5.10 Shear force versus deflection relationship Fig.4.5.11 Shear force versus deflection relationship 
of boundary beam of 2nd floor at opposite of boundary beam of 3rd floor at loading 


side of loading of Wl side of Wl 
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Fig.4.5.12 Shear force versus deflection relationship of boundary beam 
of 3rd floor at opposite side of loading of Wl 
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Fig.4.5.13 Shear force versus deflection relationship Pig.4.5.14 Shear force versus deflection relationship 
of boundary beam of 2nd floor at loading of boundary beam of 2nd floor at opposite 
side of W2 side of loading of w2 
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Fig.4.5.15 Shear force versus deflection relationship of boundary beam 
of 3rd floor at loading side of W2 


Fig.4.5.16 Shear force versus deflection relationship of boundary beam 


of 3rd floor at opposite side of loading of W2 
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test structure, difficult choices were encountered due to 
difference in design codes and practices in the United States 
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for beam-to-column connections are a few examples. 

In order to choose better and economical detailing for the 
frame parts of the full-scale test structure, half-scale beam- 
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ABSTRACT 

One-half scale reinforced concrete test specimens were 
tested to study the behavior of beam-column subassemblages 
designed under Japanese and U.S. construction practices. 

The specimens represented second-floor interior and exterior 
beam-column assemblies and top-floor exterior beam-column 
assemblies of a seven-story full-scale test structure. 

The placement of beam and column longitudianl reinforce- 
ment was the same in all specimens. The amount and arrange- 
ment of lateral reinforcement was varied in specimens follow- 
ing the U.S. and Japanese design requirements. More web rein- 
forcement in beams and columns was required by the U.S. design 
code. 

Within a deformation range predicted by a nonlinear 
earthquake response analysis of the prototype structure, the 
behavior of the two specimens designed by the two codes was 
guite similar. Specimens with slab behaved in a manner dif- 


ferent from the corresponding specimen without slab. 


* Head, Structural Division, Building Research Institute 
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1 INTRODUCTION 

A U.S.-Japan Cooperative Research Program Utilizing 
Large Scale Testing Facilities in the United States and Japan 
was initiated in the summer of 1977. The planning group of 
the program prepared recommendations(1) to improve seismic 
safety practices through studies and.to determine the relation- 
ship among full-scale tests, component tests, and analytical 
studies. 

It was recommended(1) that a full-scale seven-story 
reinforced concrete building structure representing good 
current practice by tested, and that a series of coordinated 
experiments associated with the full-scale tests ke conducted 
in Japan and the United States on reinforced concrete joint 
assemblies, walls, and frames. The results of all associated 
tests in both countries should be fully correlated with each 
other and with the results of the tests made in Japan on the 
full-scale seven-story structure. 

In the course of designing the full-scale seven-story 
test structure, difficult choices were encountered due to 
difference in design codes and practices in the United 
States and Japan. Differences is the equation for design 
shear capacity, column tie detailing, and the reinforcement 
requirements for beam-to-column connections are a few examples. 

In order to choose better and economical detailing for 
the frame parts of the full-scale test structure, it was 
recommended that half-scale beam-to-column subassemblages 
be tested using U.S. and Japanese construction practices. 


It was not intended to test the subassemblages beyond 


a deformation range expected* during the Tull-scale teste 
Results from the tests will be used to model the stiffness 
properties of structural members in a nonlinear dynamic 
analysis, which is used to determine a probable maximum 
deformation range of the full-scale structure during an 
earthquake under consideration. 
The (ese results will also be used to study the 

correlation between tke behavior of half-scale beam-column 
subassemblies and the behavior of beam-column parts within 


the full-scale test structure. 


Weel ool oPECIMENS 

2.1 Prototype Full-Scale Test Structure 

The full-scale seven-story reinforced concrete test 
structure represents a portion of a building having dimensions 
common to earthquake resistant construction in the United 
States and Japan. The lateral load resistance of the structure 
is provided by interacting structural wall and frames; two 
moment resisting frames and one structural wall-frame as | 
BOW 1 Fig. 2,2. the structure has” threéé™ spans’ in the 
direction of loading. The structure is symmetric about the 
Cenler. 

Spans are 6m, 5m and 6m, in the direction of loading, 
and 6m on centers in the transverse direction. Inter-story 
heights are 3.75m in the first story and 3.0m from the second 
to the seventh stories. The column dimensions (0.50m x 0.50m) 
are the same throughout the structure. The beam dimensions 
(0.30m x 0.50m) are common in the longitudinal direction. 
The transverse beams have dimensions of 0.30m x 0.45m. 

The amount of longitudinal reinforcement is the same in 
each story and in each span. All columns are reinforced 
with 8-D22 (No.7) deformed bars (Fig.2.2). Gross reinforcement 
ratio P_ is 1.24%. The amount of longitudinal reinforcement 
is different at beam end and center (Fig.2.2). The beam end 
is reinforced by 3-D19 (No. 6) deformed bars (reinforcement 
Foatic, OL 0.05%) at the top, and by 2-DL9 (No. 6)” deformed 
bars (reinforcement ratio of 0.43%) at the bottom. The slab 
reinforcement is the same in each floor (Fig.2.3). The slab 


is reinforced by D190 (No. 3) deformed bars in double layers. 


Additional negative reinforcement (D13, No. 4) is provided 
perpendicular to beams. 

Nominal strength of all reinforcing bars (SD35 in Japan, 
and Grades 50 in U.S.) 1513.50 Lerten (S0ksi.,. 0r 13 40MPale 
Specified compressive strength of concrete is 270KG/ ome 
(4,000psi or 36.5MPa) - 

The reinforcement details of the full scale structure 
have not been decided, but will be finalized with the infor- 


Mation obtained from the current tests. 


2.2 Dimensions of Test Specimens 

The behavior of the entire structure will be effected 
by the amount.of lateral reinforcement ands ane method Of 
reinforcement eee. The United States and Japan require 
different amounts of lateral reinforcement and reinforcement 
detailing. Thereforce, it was desired to study the effect 
of different reinforcement practices. 

Three types of beam-column subassemblies were chosen 
as test specimens for this purpose; 
(a) interior beam-column assemblies (I-series), 
(b) exterior, beam-column assémblies (E-series), and 
(c) exterior beam-top column assemblies (T-series). 
These three series of specimens are marked in the elevation 
diagram in Fig. 2.l1.b. The inflection points were - assumed to be 
located at the mid-span of beams and at mid-height of columns. 
The specimen was taken as the portion bounded by the inflection 
points. 

The dimensions of each test specimen were one-half 


those of the full-scale structure. However, the amount of 


reinforcement was not precisely scaled due to the limitation 
in available bar sizes. Slab reinforcement in the half-scale 
Specimens was placed in a single layer. 

The column and beam sections with longitudinal rein- 
forcement are shown in Fig. 2.4. Although the amount of 
longitudinal reinforcement at the end and at the center of 
a beam was not the same in the full-scale structure, the 
test specimens had uniform longitudinal reinforcement along 
the beams. Three D10 (No. 3) bars were used at the top 
(tensile reinforcement ratio of 0.65%), and two D10 (No. 3) 
bars were placed at the bottom (tensile reinforcement ratio 
of 0.43%) of a beam section. Four D1l3 (No. 4) bars were 
used in the corners and four D10 (No. 3) bars were placed at 
the middle face of a column section. The gross reinforcement 
ratio was 1.27 percent. 

Three specimens in each series were designed using (a) 
Japanese design practice with slab, (b) U.S. practice with 
slab, and (c) Japanese practice without slab. The slab 
width was chosen to be 75cm on the basis of the effective 
width of T-beams. The effective width requirements are 
comparable in the two countries. Dimensions and reinforce- 
ment for test specimens are summarized in Table 25 bend 
their details are shown in Figs. 2.5 to 2.7. 

The major difference in U.S. and Japanese design 
requirements are the amount and arrangement of lateral 
reinforcement, and the method of anchoring beam longitudinal 
reinforcement in the beam-column connection. The en digte of 


beam lateral reinforcement in the specimens was determined 


by the minimum requirements in both countries because the 

expected shear forces in the members were small. It isa 

Japanese construction practice in the exterior beam-column 

connection to bend down the beam bottom longitudinal rein- 

forcement and to anchor the reinforcement in the lower 

column (Figs.2.6 and 2.7) because the concrete is normally 

cast from a slab face to the next slab face above at a time. 

It was decided to use deformed bars for lateral rein- 
forcement. The smallest deformed bar available was D6 (No. 2), 
and was used in the specimens. 

Architectural Institute of Japan Standard (AIJ Standard) 
for structural Calculation of Reinforced Concrete Structures 
-1979-(2) requires the amount of beam lateral reinfercement 
as follows; 

(a) the web reinforcement shall be at least 9-mm diameter 
plain bars or at least 10-mm diameter deformed bars 
(D10), 

(b) the web reinforcing bar shall be spaced at not more than 
one-half the overall beam depth and at not more than 25cm, 

(c) the minimum web reinforcement ratio shall be not less 
than 0.2 per cent. 

therefore, the spacing of beam web reinforcement in Japanese 

specimens was determined to be 12.5cm, one half the overall 

Gepth of the beam secticn. 

The spacing of beam web reinforcement in U.S. specimens 
was determined as one-fourth of the effective depth of the 


beam section as specified in ACI standard 318-77(3). 


2.3 Material Properties 

High-early strength concrete was used in the Specimens. 
All specimens were cast outdoors, from the same batch of 
concrete, in the up-right position. The specimens were 
left outdoors during the curing period in Fepruary, Dut the 
temperature never got down below 4°C (39°F). 

The average compressive strength of concrete was 340kg/cm* 
(4,900psi or 33PMa), which was 20 per cent higher than the 
specified strength of 2T0Ke/ cn (4,000psi or 26.5MPa). The 
beam-column subassemblies were tested at 17 to 55 days from 
concrete casting. The increase in concrete strength was 
observed to be small during this period. 

Reinforcing bars (D6, D10, and D13) were tested in 
accordance with Japan Industrial Standard, similar to ASTM 
in U.S. The stress-strain curves of three kinds of reinfrocing 
bars are shown in Fig. 2.8. The yield stresses of D6 (No. Zee 
Dw tNO* 54" and D113 (No. 4) were ee feo om (54.7ksi or 370MPa) 
3.75ton/em* (54.4ksi or 368MPa), and 4.07ton/cm”(59.1ksi or 
399MPa), respectively. These values were approximately by 
5 to 10 per cent higher than the nominal strength of F5ton/cm- 
(Jo.0nslt Or S43MPa) . 

D1O0 and D13 reinforcing bars showed clear. yield plateaus, 


and Dé bars yielded gradually. 


3 LOADING 

3.1 Apparatus 

Varying forced displacements were applied at the free 
ends of the beams. The loading apparatus for series I (interior), 
E (exterior) and T (top) beam-column subassembly specimens 
USesnown}an fig eels 

Constant axial load was applied to the column of series 
I and E specimens, simulating the gravity load corresponding 
to the working load conditionsin.the full-scaleystructuma, 
The level of axial. stress was 38.4kg/em~ (5S7psi or 3.7 73eee 
in series I specimens, and SER Ska ien (527psi or 3.56MPa) in 
series E specimens. No constant axial stress was applied in 
series T specimens, but the column axial stress varied with 
lateral load. 

A 100-ton capacity actuator with + 30-cm stroke was used 
to apply the constant) axial load in the columneesthesioag 
was measured by a 100-ton load cell. An electrical pump with 
automatic control was used to apply hydraulic pressure in the 
1C0-ton actuator. The column axial load fluctuated during 
the €xperiment® by: P22 atoeinomith exepeci@ved) leveia 

Two 30-ton capacity actuators (+30-cm stroke) were used 
to displace the beam free ends of series I specimen. The 
same actuator was used in a series E or T specimen. Two 
30-ton capacity load cells were attached to these actuators. 

The existing heavy steel reaction frame was used for 
the test. The reaction frame is made of 900 x 300 x 16-mm 
H sections for columns, and 600 x 300 x 12-mm H section for 


beams. Some extension pipe members were added to the actuator- 


load cell assembly to fill the space between the top and 
bottom beams of the reaction frame. 

The top end of the series I and E columns and the free 
end of the eines T beam was supported by a vertical roller. 
The other ends of all specimens were either supported through 
a mechanical hinge, or connected to a loading device through 
a mechanical hinge. The location of mechanical hinges and 
roller in the specimens corresponded to the aeeuneam deaeion 
of inflection points in the full-scale structure. 

The beam of a series T specimen was subjected to varying 
amplitude axial load. In an actual structure, the shear 
force acting in a exterior column must be resisted by the 
adjacent beam and slab in the form of axial stresses. 
Therefore, the axial load acting in the beam of a series T 


specimen was thought justifiable. 


32294 Program 

Before the tests, the response of the full-scale proto- 
type structure was computed under earthquake motions. One- 
component model (4)., the model which assumed all inelastic 
deformation to ke concentrated at member ends, was used to 
represent the deformation of reinforced comcrete members 
including the structural wall. The degrading trilinear 
hysteresis model(5) was used to model the moment-rotation 
relationship at member ends. The NS component of the 1978 
Miyagi-ken Oki earthquake recorded at the Tohoku University 
Station was used for the analysis. 

The response was analyzed:(a) for a mediumintensity, 


and (b) for a major intensity earthquake, by scaling the 


acceleration amplitude of the earthquake record. A maximum 
acceleration of Near aeten (0.2g) was chosen for a medium 
intensity earthquake, and a maximum acceleration of 450cm/sec 
(0.45g) for a major intensity earthquake. 

The maximum beam end rotation in the second flcor level 
was approximately 1/200 radian from the medium intensity 
earthquake, and approximately 1/100 radian from the major 
intensity earthquake. The existing shear wall in B-frame 
(Fig. 2.1) was effective in reducing the response deformation 
of the prototype structure. With this information at hand, 
it was decided to test the specimens causing the beam end 
rotation up to 1/50 radian. It was also decided to apply 
at least two cycles of load reversals ata given displacement 
amplitude to study the deterioration in resistance. 

The loading. program is schematically shown in Fig. 3.2. 
In a series I specimen, the displacements at both beam ends 
were moved by the same amplitude, but in the opposite directions. 
The relation between the beam deformation and story deflection 


is diseussedpin pSectéyronn4 tis 


aes eee 


4 INSTRUMENTATION 
4.1 Displacement Measurement 

Two independent reference frames were attached to a test 
specimen to measure the deformation. Figure 4.1 shows the two 
reference frames attached to a series I specimen. 

The darkened rigid aluminum reference frame A was attached 
to the sides of a transverse beam immediately adjacent to the 
column. The column end lateral displacements were measured 
with respect to the frame A. 

The beam end displacements were measured with respect 
to the shaded rigid aluminum beam B, attached to the top 
and bottom of the transverse beam immediately adjacent to 
the column. 

One reference frame moved relative to the other during 
the loading due to the distortion of the beam-column joint 
panel. The angle change between the two frames approximately 
represented the shear distortion angle of the beam-column 
joint panel if the shear distortion was assumed to distribute 
Mircormly within the, joint. 

The flexural rotation in the critical region of a beam 
was defined by measuring the longitudinal deformations at 
the top and the bottom of the beam. A displacement transducer 
measured the deformation between the column face and a bolt 
embedded at 10cm from the column face. This measurement 
included a deformation attributable to the slippase of beam 
reed nal bars within the beam-column connection. 

The displacements were measured by Berita ee type 


transducers, with + 50-mm stroke and 30 -mm resolution for 


column-end displacements, beam-end displacements and beam 
flexural rotation measurements. The relative movements of 
the two reference frames were measured by the same type 
displacement transducers, but with +e5- mm? stnokegand _— -mm 
resolution. 

The location of displacement transcucers is shown in 
Fig. 4.2 for the three specimens. The usage of the transducers 


is summarized below. , 


Specimen I, Specimen E, Specimen T 
Column Deformation Oiten Ge 01,02 O01 
Beam Deformation Ginae 03 02 
Joint Panel Deformation Os-O6 O4-Os 03-04 
Beam End Rotation KizGs 2970n gmat ay Os-O6 


The displacemeuts measured by the transducers were used 
ta define the member deformations and story deformation. 
Consider a series I specimen. 

Beam-end flexural rotation 0, was defined by the beam 
longitydinal displacements in tnedcuiticad region, for example 


Tree Gs art \ Sar fn deer: ge a Be 


10’ 
Ob = ( S19 =: Dg ) v/ hy (4.1) 
in which 


h, : distance between two transducers for displacements 


C9 and S10 7 


8p : flexural rotation over a 10-cm gauge length 
The shear distortion angle, Ws , Of a beam-column 
joint panel was defined by the relative displacements of the 


two reference frames, for example 9. and O°. ina series I 


5 6 
specimen, as shown in Fig. 4.4, 
1 Pa ie pe ee (4.2) 
in which 
es : distance between two displacement transducers 


for displacements Oc and Jes 
The measured colum-end and beam-end displacements 
included the displacement contribution from the distortion 
of the beam-column joint panel. The column-tip displacement, 
ae measured from the beam top face could be defined as 


faltows, Fig..4.5; 


10} =j {6) + ea 
a 1 ‘hes V4.3) 
where 
Oo, : column-end displacement measured from reference 
frame A; 


y, : shear distortion angle of beam-column joint 
panel; 


a : one-half of height of beam-column joint panel. 


Similarly, the beam-tip displacement, op , Measured from the 


column side face was defined as follows, Fig. 4.5, 


0 EPR as fe Ns (4.4) 


iniwhich 


o, : beam-end displacement measured from reference 
frame B; 
b : one-half of width of beam-column joint panel. 


Although the lateral movement at the column tips was 
prevented during the experiment, the story displacement of 
the prototype structure could be determined from the displace- 
ment measurements of the specimen. In the prototype structure, 
the inflection points in beams will not annie tania move in 
the vertical direction. Therefore, if the beam tips were 
assumed to have stayed horizontal, the story lateral displace- 
ment between the two column inflection points can be determined 


as follows (Fig. 4.6); 


° story Pry F 0. tT M$ ha es laa) 
hee Hee + Y Soe ) 7x (ho eee (1 oe ) 
(4.5) 
in which 
wat; Cas iret and o, : column and beam tip displacements 


shown in Fig. 4.2; 


na RA 


y., : Shear distortion angle of beam-column joint 


Pp 
panel; 
he > one-half of column clear height; 
Xy one-half of beam clear span. 


Sometimes it iS convenient to know a story deformation 
angle, R, rather than a relative story displacement. The 
story deformation angle is defined as the relative story 


displacement, a , Givided by the inter-story height h.- 


story 


=a0) e 
3 story Gh sisal 
The damage of nonstructural elements in a structure due to 
earthquake motion is often related to this angle. 
Story shear, 9, waS approximately obtained from the 


moment equilibrium of the subassemblacge as given below (Fig. 4.7): 


P =e 
ge —ti 4. xi (4.7) 
2 
where 
Pits P. ; lateral loads applied at:.beam tips; 


£2: spanwhkength: 


The P-Aj effect of column axial load on story shear was 
neglected. The story shear acting in the prototype s$/ructure 
is not exactly the shear given by Fig. 4.7 because Fig. 4.7 
does not represent the configuration of the subassemblage 


within the deformed prototype structure. 
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4.2 Strain Measurement 

Strain gauges were used to monitor the behavior of 
longitudinal reinforcement in the beams, columns and slab. 
The position of strain gauges is shown in Fig. 4.8 for 
Specimen I-l. 

The strains in beam and slab longitudinal reinforcement 
were measured at the column face. The strains in column 
reinforcement were measured at beam top and bottom face 
levels. In addition, the strains in beam longitudinal 
reinforcement were measured at the center of the beam-cloumn 
joint panel to examine the deterioration of bond between the 
reinforcement and concrete in the beam-column joint. 

The strain gauges used were high elongation YL-10 (Tokyo 
Sokki Kenkyujo, TSK) gauges. The gauge length weed 10-mm, 
and the gauge resistance was 120%. After the deformation 
on a reinforcing bar was removed by file, the location for 
a strain gauge was smoothed by sand paper and cleaned by 
acetone. A strain gauge was pasted onto the cleaned bar 
surface by CN adhesive, and was coated with heated wax in 
two or three layers for water proofing. Furthermore, the 
gauge was protected by a glass-fibre cloth with B2 adhesive. 
The two types of adhesive were provided by TSK. 

Strain gauges were connected to either of two swictch 
boxes, and to an Automatic Strain Meter TDS-256 (TSK). The 


strain readings were digitized and punched on a paper tape. 
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5.1 General Behavior 

I-Series Specimens (I-l, -2 and -3) were cross-shaped, 
representing a second-floor interior beam-column subassemblage 
of “the “prototype ‘building. Specimens I-l “and I-2 had-a-slab, 
and Specimen I-3 did not. Specimen I-l, designed according 
to Japanese design code (2), had a smaller web reinforcement 
ratio than Specimen I-2, designed according to U.S. design 
Practice (s)as Shown’ in ‘Table r2eL. ©The ‘three specimens had 
the same amount of longitudinal reinforcement. Specimen I-3 
had the same web reinforcement ratio as Specimen I-l. 
The web reinforcement ratio in beams, columns and beam-column 
connection was varied by the spacing of lateral reinforcement, 
using the same bar size (D6 or No. 2). Additional sub-ties 
were placed in Specimen I-2 column to support middle logitudinal 
reinforcement, increasing substantially the web reinforcement 
ratio’in the columns (Table 2.1). 

General behavior of Specimens I-1l and I-2 was similar, 
but different from that of Specimen I-3. 

Figures’ 5/1 to 5.3 show ‘cracking ‘patterns of these three 
specimens in the first cycle (story deformation angle 
Re] P7200y7"Gn the third*cyele” (R= 1/2100), and after test. 
The story deformation angle is defined in Section 4.1. 
Tensile yielding of beam reinforcement was observed in the 
first cycle in all specimens, 

Specimens I-l and I-2 (Figs. 5.1 and 5.2) showed almost 
identical crack patterns, indicating little effect OF the 


difference in the amount of web reinforcement on the behavior 


in this range of deformation. 

Flexural cracks were first observed in the beams or slab. 
Slab cracks developed perpendicular to the beam axis extending 
over the entire slab width. Diaggqnal shear cracks in the 
beam beneath the slab were observed at approximately R = 1/200, 
starting from the slab side only. The slab contributec to 
the flexural capacity when the beam top was in tension, hence 
the larger shear force was developed to cause shear cracks. 
However, when the beam bottom was in tension, the slab 
contribution to the flexural capacity was not large enough 
to cause shear cracks. 

When the deformation increased to approximately R = 1/50, 
the cracks developed along the beam-slab boundary line (Figs. 
S206 Andeso2eC ie Wide flexural cracks in tit candela lower 
part concentrated at the column faces. No compressive crushing 
was observed in these two specimens. Minor Filta cracks 
were observed in the column at the beam faces. 

In the case of Specimen I-3, the shear force in the beams 
was limited by a relatively low ieee ae yielding capacity. 
Therefore, the cracks were all due to bending without compressive 
crushing. The-beams did not develop diagonal shear cracks. 

No flexural crack was observed in the columns, nor in the 


beam-column joint panel. 


5.2 Hystersis Relationship 
The vertical displacements were applied at beam ends of 
the subassemblage to simulate the lateral load induced in 


the prototype structure. The magnitude of vertical displacements 


hg. 


(6, ) at both beam tips were the same, but in the opposite 


B 
directions. The inflection point of the beam in the prototype 
building were assumed at the mid-span. 

Figures 5.4, 5.5 and 5.6 show story shear interstory 
displacement relationship of Specimens I-l, I-2 and I-3, 
respectively. The points at which yielding of the beam and 
slab reinforcement was first observed iS also indicated in 
the figures. 

Specimens I-1 and I-2 showed almost identical behavior 
up to a displacement amplitude of R = 1/100 rac. The two 
specimens snowed conspicuous slip behavior in the fifth cycle 
(R = 1/50) of loading. Strength decay was observed in the 
repeated cycle at the same large amplitude in these specimens, 

Specimen I-l having a larger strength loss (Figs. 

5.4 and 5.5). This slip behavior must be closely related 

to the wide flexural cracks concentrated in the beam lower 
portion at the column faces (Fig. 5.7). The subassemblage 
could not develope resistance until either of the two cracks 
closed during loading. The smaller slip behavior in Specimen 
I-2 might be attributable to the extra lateral reinforcement 
provided in the beam-column connection, rather than to the 
larger amount of beam or column att reinforcement. 

In the earthquake response analysis, the maximum story 
deformation angle was estimated to be R = 1/100 even from 
the most severe earthquake motion (Section 3.2). Within 
fbi range of deformation, all Ea I specimens showed 
fat and stable hysteresis loops. The behavior of the three 


specimens is judged to be satisfactory. 


— 1) 


When the full-scale structure is tested to the 
capacity of the actuators available at the Building Research 
Institute, Ministry of Construction, Japan, the + 50-cm stroke 
capacity of the sixth floor actuators will limit the deformation 
of the full-scale structure during the experiment. If the 
displacement is assumed to vary linearly from the base to 
the sixth floor, the average story deformation angle becomes 
1/397 .Saeqdt theasthucturess steestedeunder olad reversals 
at this displacement amplitude, the slip behavior of the 
subassemblage may cause a_ slip-type hysteresis of the 
Structure, 

The second cycle resistance at 6 centimeter displacement 
(Point B in Fig. 5.4) in Specimen I-1l was lower than the 
previous maximum resistance (Point A). This was caused by an 
error in controlling beam ends displacements at point C. 

The left beam was moved upward by 5.5 cm and the right beam 
downward by only 2.5 cm. Therefore, the left beam was 
severely damaged at the point C causing a reduced resistance 
aipousitiap. 

Hysteresis shape of Specimen I-3 was different from 
those of Specimens I-1l and I-2. Hysteresis loaps cf Specimen 
I-3 were fat and stable, and strength decay was not observed. 
The effect of slab on the hysteresis behavior was significant. 

Loading was reversed at small displacement (R = 1/400), 
in théefilrst ha Pi landictheyieast “half of the gthird cycle 
(R = 1/100). These small loops are shown in Figs. 5.4, 5.5 
and 5.6 with broken lines. These loops were thin and showed 


slip behavior in Specimens I-l and I-2. 
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The overall UereeL Gn Ge the subassemblage eS composed 
of contributions form the beam, column and beam-column JOLOe 
deformations. The beam deformation sharéd a large pertion 
of the subassemblage deformation, particularly after yielding. 
At R = 1/50, measured contributions of the beam deformation 
to the total deformation were 90-92 % in Specimens I-1 and 
I-2, and 96-98 % in Specimens I-3. Column and joint deformations 
of Specimen I-3 were small in comparison with the beam, as 
indicated by no crack observed in column and joint. In the 
case of Specimens I-1 and I-2, column deformation was larger 


than that of Specimen I-3 due to the larger force levels. 


5.3 Strain Measurement 

Strain in reinforcement was measured,mainly at the 
critical sections of the beams, columns and slab. The story 
shear-strain relationship of the six strain gauges shown in 
Pip s<0 15 stow in Figs. 5.9 to 5.11 for Specimens I-1l, 
I-2 and I-3, respectively. In the case of Specimens I-1 and 
I-2, the tensile yielding of beam bottom reinforcement was 
observed in the first or weLt8G cycle of loading (Figs. 5.9 
and 5.10). Tensile yielding of beam top reinforcement and 
slab reinforcement was observed later, in the second or 
fourth cycle. In the case of Specimen I+3, Fig. 5.11, top 
and bottom reinforcement yielded during the first or second 
cycle of loading. 

The reinforcement slippage was examined by the strains 
in beam reinforcement at the center of the heanecorarn FOINC. 


These strains remained in the elastic range, which indicated 
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that bond failure did not propagate to the center of the 
beam-columi joint (Figs. 59 -to ss Pr) : 

Common among all I-series specimens, strains in the 
column reinforcement remained in the elastic region. 
Strains in column reinforcement in specimens I-l and I-2 
Changed their signs due to large moment reversals while in 
Specimen I-3, strains in column reinforcement remained com- 
pressive due to high axial load and lower bending moment 
reversals. 

The distribution of strains in slab reinforcement at 
the critical section is. shown, jin Figs. -5..12_ and 5.13 se 
strain distribution was relatively uniform prior to the 
yielding (R = 1/200). Once the yielding of a reinforcing 


bar started, the distribution shape became irregular. 


5.4 Measured and Calculated Strengths 
Measured and calculated yield and maximum moments at 
the critical beam section are given in Table 5.1. Measured 
yield moments were determined when the first yielding of 
tensile reinforcement (beam or slab) was observed. Measured 
maximum moment were the largest moment observed during loading. 
Yield and maximum moments were calculated in two ways; 
(1) Approximate yield moment equations suggested in 
Commentary of AI J Standard (2) 


for beam, 


Me = 0.9 PEGaa (5S 
and for column, 
M = 0.8 oO ‘ “e ° 
¥ 8) as i +O 5 ae Cu N/bDF .) So 


where, 


a, ? area of tension reinforcement including slab 
reinforcement; 

oy wetecistlc fLetdgottoss "OrUreInmrolceing = steel, 

a : beam effective depth; 

D : column overall depth; 

N?.>. ecolumn) axial. load; 

b : column overall width; 

Fo : compressive strength of concrete. 


(2) Flexural Theory 

A linear strain variation across the section was assumed 
and the stress-strain relationship of the longitudinal steel 
(including strain hardening) and concrete were considered. 
Figure 2.8 shows the stress-strain relationships for steel 
used in computation. 

Maximum moments in Table 5.1 were evaluated at the 
deflection when the maximum moment was observed in test. 
The deflection was evaluated on the assumption of linear 
curvature distribution along the beam. Maximum beam-end 
moment was observed at a relative story displacement of 
approximate 4 cm (R = 1/50), which corresponds to a beam-end 
curvature of 1.0 se Soom weds, Therefore, maximum moments 
in Table 5.1 were calculated for a curvature of 1.0 x 1077 
(l/fem). Extreme compression fiber strains at this curvature 
are listed. Table 5.1 shows that calculated values are 
almost identical to the corresponding test values. 


In Specimens I-l and I-2, measured beam top yield moments 


are lower than the calculated values; 


this is because yielding of all the slab reinforcement 
occurred slightly after the yielding of the beam reinforcement. 
In the case of beam bottom tensile yielding, the values of 
measured yield moments in Specimens I-1l and I-2 were about 


30 percent larger than that in Specimen I-3. 


2. ero Unmeal Vi 

Second-story interior beam-column subassemblages were 
tested in this. serres to study the effects of (a) the amount 
of web reinforcement in beams, columns and beam-column 
connection, (b) the sub-ties for column middle reinforcement, 
and (c) the slab. 

Within the deformation range predicted by a nonlinear 
earthquake response analysis, the amount of web reinforcement 
suggested by the Japanese Code (2) was sufficient. The crack 
patterns and stiffness characteristics of the subassemblages 
were not Significantly affected by the amount of the web 
reinforcement used in the study. Slightly better performance 
of Specimen I-2 was observed in a displacement range where 
the story deformation angle R exceeded 1/50, probably 
because of better confinement in the beam-column connection. 

In the specimens tested, the major damage occurred 
within the beams. The damage of the columns was limited to 
minor flexural cracks. Therefore, the sub-ties, which gave 
support to the column middle reinfocing bars, did not affect 
the behavior of the specimens. 

The slab contributed to the resistance of the subassemblage. 


The yield moments observed in Specimens I-l and I-2 witha slab 
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were more than twice that observed in Specimen I-3 without a 
slab, when the slab was stressed in tension. When the slab 

was in compression, the yield moments observed in the Specimens 
with a slab were more than 1.3 times that in the specimen 
without a slab. Due to the slab contribution to resistance, 
diagonal shear cracks were observed in the specimens with a 
Slab. Slip-type hysteresis behavior was also noted in the 
specimens with a slab. especially at large displacement 


reversals. 
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6 ~RESULIS, OF E-SERLEol toon 
6.1 General Behavior 
E-Series Specimens (E-l, -2 and -3) were of T-shape, 
representing a second-floor exterior beam-column subassemblage 
of the prototype full-scale building. Specimens E-1 and E-2 
were constructed with a slab, and Specimen E-3 without a slab. 


The bar size was kept the same in the three Specimens (E-l, 


2 and 3). The amount of longitudinal reinforcement in beams 
and columns was the same in the three specimens. Specimen 


E-l was designed in accordince with the Japanese design 
practice (2), and Specimen E-2 according to the U.S. design 
practice (3). A major difference in the two Specimens (E-l 
and 2) was the amount of lateral reinforcement in beam, 
columns and Posneootdan connection as shown in Table 2.1 and 
jE ye 2.6. The web reinforcement ratio was varied by using 
different spacing of lateral reinforcement. Additional 
sub-ties were provided in Specimen E-2 columns to support 
middle longitudinal reinforcement. The beam longitudinal 
reinforcement was anchored dcwnward in the beam-column 
connection of Specimen E-1 as shown in Fig. 2.6.a, following 
the standard Japanese construction practice. The beam and 
columns of Specimen E-3 were identical to those of Specimen 
Bol. 

Figures 6.1 to 6.3 show crack patterns of the three 
spceimens: 
(a) at a relative-story displacement angle (inter-story 
Gisplacement divided by inter-story height) of R = 1/200 


in the first cycle, (pb) eae n= 1/100 in the third cycle, 


ie 


and (c) after test. 

Flexural cracks were observed first in all Specimens. 

No flexural crack was observed in columns, but small vertical 
cracks were observed above the slab near the beam-column 
joint in Specimens E-l and 2. These cracks were not Significant. 

The behavior of the beams was dominated by flexure. 
Wider craks were observed on the slab surface of Specimen 
E-l than that of Specimen E-2. 

In Specimen E-2 (Fig. 6.2), narrow cracks appeared on 
the entire slab face in the first cycle. In Specimen E-1l 
(Fig. 6.1), wide cracks developed near the beam-column joint, 
and cracks spreaded toward loading point in later loading 
cycles. 

Diagonal shear cracks were observed in the web of beams 
in Specimens E-1l and E-2 only in one direction, starting from 
the slab side. These cracks were developed when the slab 
was in tension, causing the slab longitudinal reinforcement 
to yield. No diagonal crack was observed in the beam of 
Specimen E-3 because the flexural capacity of the beam was 
low. 

Cracks observed on the slab surface of Specimens E-1l 
and E-2 near the joint were perpendicular to the beam long- 
itudinal axis, whereas the cracks near the loading point 
were in a fanshape pattern with a pivot at the loading point. 

A wide flexural crack was observed in the beam lower 
portion along the column face, in Specimens E-1l and E-2 when 
the Specimens were displaced by R = 1/100. The second wide 


crack appeared approximately 10 cm from the first crack in 
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the two Specimens at R = 1/50. 

Spalling of the cover concrete occurred at the bottom. 
face of the beam adjacent to the beam-column joint at the 
final stage of loading in Specimen E-l. No reinforcement 


bucking was observed. 


6.2 Hysteresis Relationship 

Force-deflection relationship of the specimens is 
discussed in terms of story shear and inter-story displacement. 
The method of computing story displacement is described in 
Section 4.1. Figures 6.4, 5 and 6 show the relations for 
Specimens E-l, 2 and 3, respectively. The inter-story 
displacement amplitudes of corresponding loading cycles 
were not identical in Specimens E-l and E-2 because the 
loading control by the beam end displacement was not success- 
ful during the two tests. 

The amplitudes of story shear in the positive and negative 
directions were different in Specimens E-1l and E-2, directly 
reflecting the strength contribution from the slab. When 
the slab was in tension, the sign of story shear was defined 
positive. The maximum positive shear was more than 2.5 
times the maximum negative shear. The difference was also 
attributable to the more negative reinforcement in the beam. 

The load-deflection relationship of Specimen E-1 was 
almost identical to that of Specimen E-2 up to a displacement 
amplitude of approximately 40 mm (story displacement angle 
of R= 1/50). In the second cycle at R = 1/50, Specimen 


E-l showed a larger reduction in resistance at the peak 


displacement than Specimen E-2. The resistance of Specimen 
E-l at story displacement angle R = 1/50 was 90 per cent of 
the maximum resistance, whereas the resistance of Specimen 
E-2 was. 96 per cent. The resistance of Specimen E-1l deterio- 
rated further, and became lower than that of Specimen E-2 
beyond 40-mm displacement. Larger deterioration of resistance 
must be attributable to the smaller amount of web reinforce- 
ment in the beam and beam-column connection and the method 

of anchoring the beam logitudinal reinforcement in Specimen 
E-l. However, the deterioration appeared in a displacement 
range much larger than the expected maximum earthquake 
response displacement. 

Both Specimens E-1 and E-2 showed a slip-type hysteresis 
shape in the positive direction (slab in tension) at the story 
displacement angle greater than R = 1/100. Up to this dis- 
placement, the hysteresis loops were stable. 

~ Specimen E-3 showed a hysteresis loop shape different 
from Specimens E-1l and E-2; i.e., the hysteresis loop was 
fat and stable, cominated by flexural behavior. This 
Gesirable behavior was obtained because beam shear was 
limited by low flexural capacity of the beam without the 
slab contribution. 

When the first and repeated cycle hysteresis loops at 
the same displacement amplitude are compared in Specimens 
E-1 and E-2, the hystersis loop in the repeated loading cycle 
was much thinner than that in the first cycle especially in 
the posititve loading direction (Figs. 6.4 and 5). ‘This 


reduction in resistance must be attributable to the formation 
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of wide flexural cracks in the previous cycle. The difference 
in the hysteresis loops of the first and repeated cycles was 
much less in Specimen E-3. 

The contribution of column and joint-panel deformations 
to the story displacement was observed to be small. The 


story displacement was governed by the beam deformation. 


6.3 Strain» Measurement 

Strains in the reinforcement were measured at the critical 
sections of the beam, columns and slab. The location of the 
gauges on beam reinforcement are shown in Fig. 6.7. 
These strains are plotted against story shear, and shown in 
Figs. 6.8 to 10. The gauge numbers are cGefined in Fig. 6.7. 
The strain of beam longitudinal reinforcement mn the center 
of the beam-column connection was also measured to examine 
the bar slip. 

The beam longitudinal reinforcement at the center of the 
beam-column joint did not yield. This indicates that the 
bond failure did not propagate to the center of the beam-column 
JOLntSe 

The column reinforcement did not yield because the 

yield moment of the coiumn was higher than that of the beam. 

The distribution of strains in slab reinforcement is 
shown in Figs. 6.11 and Fig.6.12 for Specimens E-1l and E-2. 

The beam reinforcement was observed to yield first 
during the test, and then the slab reinforcement adjacent 
to the beam yielded in the first loading cycle (R = 1/200). 


The yielding of slab reinforcement propagated outward. The 


outermost slab reinforcement yielded in the loading cycles 


at story displacement angles of either R = 1/100 or LOU 


6.4 Measured and Calculated Strengths 
All Specimens were designed to have flexural failure in 
the beam. Calculated yield and maximum moments at the beam 
critical section are compared with the observed values in 
Pere ye 
Measured yield moments were defined as the moment 
observed when the yielding of tensile reinforcement (beam 
or slab) was first reported. Maximum moments were taken as 
the largest moment attained during all cycles of loading. 
(1) Approximate Yield Moment Equations Suggested in Commentary 
of AIJ Standard (2); 
for beam; 
M = eee ro, 1) 


y 
and for column; 


= Se Baga ott + O0.5ND(l1 - N/bDF .) (Gn 2) 
where 
a, = area of tension reinforcement including 
slab reinforcement; 
oye = “tensile yield stress of reinforcing steel: 
ad: beam effective depth; 
D : column overall depth; 
N : column axial load; 
b : column overall width; 
in > compressive strength of concrete. 


(2) Computation Based on Flexural Theory. 
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A linear strain variation across the section was assumed. 
The stress-strain relationship of the longitudinal steel and 
concrete were considered in the analysis, as shown in Fig. 2.8. 

Equation 6.1 gave the same beam positive yield moment 
for the three specimens, whereas the flexural theory (Fig. 2.8) 
Calculated the positive yield moment of Specimens E-1l and E-2 
approximately 12 per cent greater than that calculated by 
Eas G6alpivefileotangetic contribution of slab concrete and 
compressive reinforcement in beam and slab. The positive 
yield moments observed in Specimens E-1l and E-2 were 
approximately 30 per cent greater than that calculated by 
EGcmm6. ais: 

The calculated beam negative yield moments by the two 
equations are comparable, and gave good maui ea of, the 
observed yield moments of Specimens E-2 and E-3. The observed 
positive yield moment of Specimen E-1l was approximately 8 per 
cent lower than that of Specimen E-2. 

The flexural theory predicted well the observed positive 
and negative yield moments of Specimen E-3 without a slab. 

The ultimate moment was calculated by the flexural 
theory at an arbitrarily chosen curvature of 0.001 l1/cm. 

The calculated beam positive ultimate moments were almost 
identical to the observed maximum moment for the three 


specimens. The calculated beam negative moments were 6 to 


10 per cent lower than the observed values. 


6.5 Summary 


Second-floor exterior beam-column subassemblages were 


tested in this series to study the effects of (a) the amount 

of web reinforcement in beam, columns and beam-column connection, 
(b) the sub-ties for column middle reinforcement, (c) the 

method of anchoring beam longitudinal reinforcement in the 
beam-column connection, and (d) the slab, on the behavior 

within a practical deformation range. The effects of the 

amount of web reinforcement and the method of beam longitudinal 
reinforcement anchoring could not be separately identified 

in this series of tests. 

A nonlinear earthquake response analysis of the prototype 
structure predicted the maximum story deformation angle to 
be. approximately 1/100 radian even from the severe eartnquake 
motion. Within this deformation range, the amount of web 
reinforcement suaqgested by the Japanese Code (2) and the 
method of anchoring the beam longitudinal reinforcement 
commonly done in Japan was sufficient. Better hysteretic 
performance was observed from the specimen (I-2) designed 
using the U.S. practice beyond this deformation range. 

No crushing was observed in the critical region of 
column in any of the three specimen. Therefore, the sub-tie 
to.support middle reinforcing bars of column section was 
judged to be redundant. | 

The slab contributed significantly to pie Sie ereoes 
of the subassemblage, especially when the slab was stressed 
in tension. The increased resistance caused diagonal cracks 
in the beam, and wide flexural efack in the lower part of the 
beam along the column face, causing a slip-type hysteresis 


loop at large deformation amplitudes. 
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7 SRESUBTSIGOE sT-SBERLES “TEST 
7.1 General Behavior 
T-Series Specimens (T-l, 2) represented a top floor 


exterior column-beam subassemblage of the prototype building 


as shown in :Figs« 27 °(a)icand=(b).> Specimen T-l was) designed 
according to the Japanese design code (2). Specimen T-2 was 
designed according to *rheats, Swi design code, .(3)*). The:U ay. 


design code (3) requires more web reinforcement (Table 2a; 
sub-ties to support column middle reinforcing bars (Fig. 2.7.b), 
and beam longitudinal reinforcement to be anchored within 
the beam-column connection (Fig. 2.7.b). 

In this series, the beam and column were subjected to 
changing axial force in proportion to their Shear. force 
because of loading system shown in Fig. eke) The axial 
load in a beam is not unrealistic because the shear acting 
in the exterior column must be carried by the connecting 
beam and slab in the form of axial thrust. 

Figures -7.1 and 7.2 show cracking patterns of the two 
specimens in (a) the first cycle at story deformation angle 
R = 1/200, (b) the third cycle R = 1/100, and (c) the final 
stage of loading. Crack behavior of the two test specimens 
was as follows; (a) flexural cracks occurred first in the 
column and then in the beam, and then (b) inclined shear cracks | 


occurred in beams when the slab was in tension. 


When the slab was subjected to tensile stress (negative 
moment), the number of flexural cracks in the column was 


larger than that in the beam. Flexural cracks were observed 


throughout the height of the column. Wide column cracks 
appeared near the beam-column connection. Flexural cracks 
also occurred in the slab. Inclined shear cracks appeared 
in the beam web immeciately below the slab after the third 
cycle of loading’. 

A wide flexural crack was observed in the column 
immediately below the beam-column connection in Specimen T-2. 
Such a wide crack was not detected in Specimen T-l. This 
difference may be caused by the method of anchoring the beam 
reinforcement in the joint; i.e. the beam longitudinal bars 
were extended beyond the beam-column connection, and crossed 
the crack in Specimen T-1, reducing the crack width, whereas 
only the column reinforcement crossed the wide crack in Specimen 


T2etigs. /.leand.2).. 


When the beam bottom was in tension (positive moment), 
flexural cracks started from the bottom of the beam. The 
number of cracks appeared in the column was smaller than in 
the beam for this case. 

General crack patterns oF the twc specimens were similar 


at the corresponding loading stage. 


7.2 Hysteresis Relationship 
Figure 7.3 shows story shear - interstory displacement 


diagrams of the two specimens Figures 7.4 and 7.5 show beam 


shear - deflection relationship. The points at which the 


longitudinal reinforcement in the beam, calumn and s‘lab was 


first obserced to yield are indicated in the figures. 


In Specimen T-1, loading was controlled by the inter- 
story displacement, while loading of Specimen T-2 was controlled 
by the beam deflection. Therefore, the loading history of 
the two specimens was different. The hysteretic behavior of 
the two specimens could not be comparec directly. 

The deflection of the column contributed appreciably to 
the story displacement in the two specimens. This can be 
noted by the difference in the shapes cf the story shear- 
displacement relationship (Fig. 7.3.a. and b) and of the 
load-deflection relation of beam (Figs. 7.4 and.7.5). 

The lateral load Carring capacity of Specimen T-1l was 


higher than that of Specimen T-2. 


7.3 Strain Measurement 

Strains in the longitudinal reinforcement was measured 
at the critical section of the beam, column and slab. The 
yielding of beam bottom reinforcement was observed in the 
third cycle of loading ( R = 1/100) in Specimen T-1, and in 
the first cycle. ( R'=-1/200) in Specimen T-2. The yielding 
of beam top reinforcement and slab reinforcement was observed 
in the fourth cycle ( R= 1/50 ) in Specimen T-1, whereas 
the yielding of beam top reinforcement and slab reinforcement 
was not observed in Specimen T-2. Instead, the column long- 
itudinal reinforcement yielded in the second cycle ( R = 1/100). 
The middle longitudinal bar in the column also yielded in 
Specimen T-2. 

The story shear-beam longitudinal bar strain relations 


(Figs. 7.7 and 7.8) of the two specimens are similar at the 


=i SiGe 


corresponding critical sections. 
Thesstorys shean-cod umnslongitudinalwstrains relations 
are shown in Figs. 7.9 and 7.l0sed Thel distri bitaon=of strains 


an slab reinforcement is showmiim Pigs’. 7.412 ahdl7<12. 


7.4 Measured and Calculated Strengths 

Measured and calculated yield and maximum moments at 
the critical section of beam and column are given in Table 
5.1. Measured yield moments (positive and negative) were 
determined when the yielding of longitudinal reinforcement 
(in beam, or column, or slab) was first observed. Measured 
maximum moments (positive and negative) were the largest 
moment observed during the test. 

Yield and maximum moments were calculated in the two 
ways: 
(1) Approximate Yield Moment Equation in Commentary of AIJ 

Standard (2). 

Axial force in the beam and column varied proportional 
to applied load. Therefore, the equation for a column was 


used for both the beam and column. 


ae = 0.8a, ak oe O53 ND, 1( mes N/DDF ) @iick) 
where 
a, i area of tension reinforcement including slab 
reinforcement; 
oy : tensile yield stress of reinforcing steel; 
D : overall column (or beam) depth; 
N : gross column (or beam ) axial load;': 
b : overall column (or beam) width; 


Fy : compressive strength of concrete. 
The axial load and bending moment were related in a Specimen, 


Mo = 7620:55 N for beam 


Lis SEN for column. 

An iterative method was used to solve for the yield moment 
and corresponding axial load. 

(2) Computation Based on Flexural Theory 

A linear variation of strain across the section was 
assumed. The stress-strain relations of the steel and 
concrete were considered in the analysis as shown in Fig. 2.8. 
Axial loads corresponding to the calculated yield moment 
(Eq. 7.1) and the observed maximum moment were used in the 
analysis. 

As mentioned in the last section, the yielding was not 
observed in the beam top reinforcement nor in the slab rein- 
forcement. Instead, the column longitudinal reinforcement 
yielded when the load was applied downward. 

When the slab was in compression, the calculated and 
yield moments are comparable. A significant difference can 
be observed between the measured maximum moments of the two 
specimens in this direction. 

When the slab was in tension, a large difference was 
observed between the yield moments measured in the two tests. 
This was caused because Specimen T-l yielded in the beam, 
whereas Specimen T-2 yielded in the column. The approximate 
AIJ equation estimated approximately 15 per cent higher 
value in Specimen T-2. The flexural theory predicted the 


yield and maximum moments conservatively. 


It is considered that such a different hinge mechanism 
was caused by the different anchoring details for beam rein- 
forcement. In Specimen T-l, the beam reinforcing bars were 
anchored into the column top. Therefore, it was effective 
in resisting bending and reducing the crack width at the 


column top. 


8 SUMMARY 

A full-scale reinforced concrete structure is to be 
tested as @ part, of @the USSe—-Japan Cooperative Research 
Program Utilizing Large Scale Testing Facilities in the 
United States and Japan. In the course of designing the 
full-scale test structure, difficult choices were encountered 
due to difference in design codes and construction practices 
in the United States and Japan. 

One-half scale reinforced concrete test specimens were 
tested to study the behavior of beam-column subassemblages 
designed under the United States and Japanese construction 
practices. 

The objectives of the investigation is not to decide 
upon the general betterness of the design. codes used in the 
two countries. Eut the study was carried out simply to 
find the minimum sufficiency of design details for a reason 
of economy within the domain defined by the capacity of 
testing facilities of the full-scale structure and by the 
response of the. full-scale structure under a probable major 
earthquake motion. 

A nonlinear earthquake response analysis of the full- 
scale structure was carried cut to estimate the maximum 
deformation under a probable major earthquake motion. The 
maximum story deformation angle of the structure was 
calculated to be approximately 1/100 radian from an earthquake 
with the peak ground acceleration of 0.45g. This relatively 
small displacement amplitude may be attributable to the 


stiffness of a structural wall in the structure. 


ae ees 


The major differences in the United States and Japanese 
design codes, if applied to this particular structure, are 
(a) the amount and arrangement of lateral reinforcement, 

(b) the placement of sub-ties tc support the intermediate 
reinforcement in a column, and (c) the method of anchoring 
beam longitudinal reinforcement in the beam-column connection. 
The amount of lateral reinforcement in the specimens was 
determined by the minimum requirements in the two codes, 
because the expected shear forces in the members were small. 
In this case, the U.S. design code (3) required larger amount 
of web reinforcement in the beam, column and beam-column 
connection, and also required the sub-ties in the column. 

It is a Japanese construction practice in the exterior 
beam-column connection to bend down the beam bottom longitudinal 
reinforcement and to anchor the reinforcement in the lower 
column because the concrete is normally cast from a slab 
face to the next slab face above at a time. 

Three types of beam-column subassemblies were chosen 
to study the difference in the behavior of Specimens designed 
by the requirements of the two codes; 

(a) second-floor interior beam-column subassemblies, 

(b) second-floor exterior beam-column subassemblies, and 
(c) top-flcoor exterior beam-column subassemblies. 

The effect of slab on the behavior of the subassemblies was 
also studied. 

Within the deformation range predicted by a nonlinear 
earthquake response analysis, the larger amount of web 


reinforcement and the placement of sub-ties in a column 


required by the United States design code did not improve the 
performance of the specimen because the deformation was not 
large enough to cause wide shear cracks, nor to cause the 
crushing of concrete in beams and columns, nor to cause the 
buckling of column reinforcement. Therefore, the minimum 
recuirements suggested by the Japanese design code (2) were 
judged to be sufficient. 

Slightly better performance was observed from the specimens 
designed under the United States design code in a displace- 
ment range where the story deformation angle exceeded 1/50. 

The slab contributed significantly to the resistance of 
the subassemblage, especially when the slab was stressed in 
tension. The increased resistance caused diagonal shear 
cracks in the beam, and a wide flexural crack in the lower 
part of the beam along the column face, causing a slip-type 
hysteresis loop at a story deformation angle greater than 


1/100. 


« hd « 


(1) 


(2) 


(3). 


(4) 


(5) 
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Fig.2.5.(Continued): Details of I-Series Test Specimens 
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Fig.2.6. (Continued): Details of E-Series Test Specimens 
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(b) E and T Series Test 


Fig.3.1. (Continued): Loading Apparatus 
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FOREWORD 


BRI Research Papers are published in English by the Building 
Research Institute and distributed to major research organizations and 
libraries throughout the world. I consider that the introduction of BRI 
research results to building researchers and engineers throughout the 
world is one of our roles for international research cooperation. 
Therefore, it gives me great pleasure to announce that the Research 
Paper No.100 has just been published. 

BRI is the only national organization of the Japanese Government 
concentrated upon housing, planning and the building science along the 
following five goals: 

btjolhesprevention of disasters, 

2. The improvement of the living environment, 

3. The rational organization of building production and the 

development of new building techniques, 

4. The effective use of energy and resources, 

5. The promotion of international research cooperation. 

Since 1960 when the BRI Research Paper No.1 was published, the main 
research results achieved by the Building Research Institute have been 
introduced through these publications. 

Ouse research “approach is extending from "how to build good 
buildings" to "what are good buildings". In this point of view, this 
Research Paper No.100; "ASSESSMENT OF AN URBAN SAFETY EVALUATION 
FUNCTION" by Dr. Yoshitsugu Aoki is a new type of research which belongs 
to the latter approach. This study Surveys what safety level in an 
urban inhabitants could expect using a skillful mathematical technique 
called urban safety evaluation function. 

When we think of buildings or urban environment from the stand 
point of safety, we must consider a kind of balance, which is safe and 
not too safe. The excessive safe design of building takes too much cost 
and sometimes loses the architectural function. These excesses “could 
alternatively be put towards improving the other items for total our 


living environment. Thus, I consider, the safety evaluation function 


4 


proposed by Dr. Aoki in this short paper would be a very important 


criterion for evaluating a well-balanced approach for buildings. 


Katsurou Kamimura 
Director General 


Building Research Institute 


April 14 1983 
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ASSESSMENT OF AN URBAN SAFETY EVALUATION FUNCTION 


Yoshitsugu Aoki 


Synopsis 


The aim of this paper is to make clear the concept of safety 
evaluation of urban area and to identify an urban safety evaluation 
function. 

Firstly, the concepts on the urban safety level are analyzed. The 
definitions on a loss vector, a frequency function, a safety evaluation 
function, and a safety level are given. 

secondly, the formula on the relationship between the safety 
evaluation function and the preference of individuals are obtained 

Furthermore, the safety evaluation function are estimated through 
the survey of preference. By the criterion AIC, the most reasonable 


safety evaluation function is a first-order threshold type function. 


a a a 


Dr. of Engineering. Deputy Head, Planning and Information Department, 


Building Research Institute. 
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1. INTRODUCTION 


Satetyoe1s, (one ).ofs the semostis importante i factors.» in 2a’ omodern 
high-density urban environment, and one which might be improved by urban 
planning or government policy. However, the concept of urban safety is 


(1) leas 


not easily understood or measured with sufficient confidence. 
difficult to express numerically the total urban safety level. 

itewer COULA “Express urban, satety as a function “of physical 
characteristics, and improve those characteristics in the actual 
planning process, we can formulate the urban safety planning problem as 
an optimization problem. We should expect to obtain useful information 
from the result of this optimization procedure. 

Generally, the urban safety level is expressed as a function of 
life loss, property loss, and service loss in an urban area. Research 
on estimation of typical losses is frequently performed. With the 
results of loss assessment studies, we could obtain an actual urban 
safety function, for use in the decision and planning process of urban 
development and renewal. However, this would be possible only if we 
could define and estimate a function which evaluate the above losses 


numerically and synthetically. 
2. PROBLEM FORMULATION 
2.1. Definition of Urban Safety Level and Safety Evaluation Function 

In the usual case, safety level is expressed either as a loss or as 
@ypotential risk, of loss. That» is, there are two types of safety 
concepts in our perception, as follows; 


A) safety level = -(loss) 


or 


B) safety level -(possibility)x(loss). 


The distinction between these two types of concepts is important in 
safety analyses, but is not clearly distinguished in our perception. 
Furthermore, the above loss is a kind of vector which includes 
components for life, property and service loss. Then, in our 
perception, there is a mapping from a loss vector to a scalar value. 

In order to formulate the above concepts, some definitions are 


needed. 
(Def. 1) HOSS" VECtOre x. a [x 
(Def. 2) Frequency, function: fix), .feb + (0-1), 
(Def. 3) Safety evaluation function: g(x), g:L >R 
(Def. 4) Safety level: s = f ja GS ah Si Pere Bae 

L 


Note that the frequency function is defined for each urban area and 
each hazard type. Therefore the safety level is determined for each 
urban area. 

The safety level for each urban area could be calculated using the 
given frequency function of that area if we have available the safety 
evaluation function. Thus, our main objective is to obtain the safety 
evaluation function. 

The safety evaluation function might satisfy the following 
condition because our results ought to be effective at decisions in an 


actual planning; 


(CSL) The safety level of an urban area k is greater than that)or 
area k' if and only if area k is preferable more than area k'! 


from a viewpoint of safety. 


However, it is possible that preference for an urban area from a 
viewpoint of safety is different for each individual. That is, the 
Safety evaluation function varies with each individual. Then, we must 
define a (social) safety evaluation function g (unsubscripted), where 


each individual i has a safety evaluation function g.- 


es get 


cel) 


(Def.5) (Social) Safety evaluation function g(x) = E.(g; 0). 


The safety level of an urban area k calculated with safety 


evaluation function of individual 1: gs > is, preater than one 


of area k' if and only if area k is more preferable than area 


k' from a viewpoint of safety for individual i. 


2.2. Relation between Safety Evaluation Function and Preference 


Individual 


The safety evaluation function g(x) can be _ expressed 


following way according to Def.5; 


Cr 


g(x) 


I 
ihe) 
— 

tad 
— 
+ 
ca?) 
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an ect Rew eewmeaes CC DB ecw Ce Come KH eRe eKeeeeeweeeoee se ew © & 


of.4an 


in the 


Furthermore, the safety level of urban area k evaluated by 


individual i is expressed as; 


N 
| 


far) tse.) dx, 
J £ (x)+g(x) dx + J Ee) ce Fdx oe J yab eis a 


The variable e, expresses a kind of personal feeling for safety, 


fection Of any urban area. In this case, the safety level 


area k evaluated by individual i can be expressed as follows; 


where Si. denotes the safety level of area k defined by; 


Ss. = 3 FO) 80x) re ie ee Mel aise nas ape 
L 


a ins 


and it is probable that this variable is independent from the frequency 


of urban 


On the other hand, the preference for urban area k over area k' 
Pra from a viewpoint of individual safety, becomes a random variable 
under the condition that the number of individuals is sufficiently 


large: 9fihat* is; 


Note that the variable e. is a random variable in this case. 
In the case ‘such that we can assume the distribution of random 


variable e. to be normal, we obtain the following equation. 


sei ah ((s,-S,,)/9). DacleaRINe oles ds avenetee [cuete: ate tree ene Rees (7) 


where 0 is the standard deviation of random variable e., and @ denotes 
the standard normal distribution function. This equation and Eq. (5) 
show the relation between the social safety function and the preference 
of individuals. 

In order to define a numerical value for safety level, we can 


define an origin Xo and a safety level scale in the following way. 


Equation (7) above can be written in the following way with the inverse 


function of the standard normal distribution function. 


By the least square method, we can estimate the safety evaluation 
function that minimizes a following value, if we could obtain the data 


Pee by minimizing Q given as 


wey ye 


3. CASE STUDY 


In this section, we shall examine the safety evaluation function in 


the case that the loss vector has only one component; life loss. 
3.1. Survey of Safety Preference 


In order to obtain the data for the ratio Pree we have designed a 
preference survey. The data about life loss by accidents for 6 urban 
areas, as shown in Table 1, have been presented to 90 individuals who 
have lived in a downtown area of Tokyo city. We asked them for their 
Seaecy preterence above 6 areas: gAdsafety preference ordering for 6 
areas was obtained from each individual. From these answers, the ratio 
Pixs was calculated as Table 2. 

Dna the other hand,__the—frequency—-function—of— 1ifte—loss—for--each 


urban area can be constructed consistent with the data presented to 


individuals as shown in Table l. 
3.2. Estimation of Polynomial Type Safety Evaluation Function 


When we assume that the type of safety evaluation function is 


_ polynomial; 


m 
g(x) =} ae ee ts Q. AE ee Cae | | Serene (12) 
1= 


m ; 
i 
= Lay J x "£.O9 WD oak lak Re aie nid AO ier eins Soe eee (13) 


where the integral represents the i-th moment of loss. ( It remains for 
the analyst to choose the order of the function which, yields the best 
results. ) 

The condition of the least square estimation, i.e., the condition 
to minimize Q defined by Eq.(11) is; 

8Q _ 


poe Or d= 1m: ie.tO. E OMO EL eae See a Ce Re (14) 


1 rm 


Table 1. Data Presented in Safety Preference Survey 


Life Loss per Year 


Urban Areas 


Year 

A B C D E 
1970 10 10 9 8 5 
1971 12 13 10 1 7 
1972 8 6 7 5 3 
1973 13 16 10 12 8 
1974 ag 9 8 7 5 
1975 14 20 1 16 10 
1976 7 4 6 3 2 
1977 9 i 6 4 
1978 10 8 6 4 
1979 1 12 9 9 6 


~~ 


seal 
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Table 2. Ordering Relation Matrix for Safety Preference Survey 


The value of Pek! 


Urban Area Code ] 2 3 4 6 6 
and ee a a a en 

Name of Area A B C D E F 
len ee ee ee ee ae 

] A - 

2 B O22 - 

3 ‘k 0.94 0.99 - 

a D 0.58 O99 0,34 - 

5 E 0.85 0.99 0.69 0.99 - 

6 F 0.64 0,97 0.49 bathe, O.39 ~ 


a ae 


or 


m : ‘ 

= il a 1 

(L&R gabe art cat ad xe xara x) alert 

oe . ns 2 ‘ { : ‘ia, loot 

’ (15) 

x [ [x £, CO dx- fx" , ekdx Hig kl 90: 
L i 


This set of equations is a m-th order linear equation system for unknown 
parameter a. and it should be solved easily. 
The results using orders 1, 2, 3 and 4 are shown in Table 3 and 


Fig.1. These results will be discussed in a succeeding section. 
3.3. Optimal Order of the Polynomial Type Safety Evaluation Function 


It remains to determine what order is the best for the polynomial 
type safety evaluation function in the results above. Dr. Akaike (7? 
proposed a criterion for such a problem, named AIC or Akaike's 
Information Criterion. 

According to his theory, the best model would have the lowest value 


of AIC, defined by the following formula. 


AIC = -2 In( maximum likelihood of model for given data ) 


+ 2 ( number of parameters in model ae 


Injour case, the polynomial type safety evaluation function has m 


parameters and the likelihood of the model for the preference data P 


kk! 
can be written as follows; 
n i { oT (PL) OSs Sir) ie 
Pa Lo Dep ee 
kk! 25% 2p 
in which n is the number of summations in Eq.(15). In our safety 


preference survey, the number of urban areas evaluated by individuals is 


6. Therefore, 


= oes 


Table 3. Estimated Polynomial Type Safety Evaluation Function 


Cr eee 


Model Code Order of Model Estimated Model 


i 1 g(x) = -0.03150x 


0.05799x - 0.00571x° 


ine) 
ine) 
wo 
— 
x 
— 
iT] 


0.01252x + 0.00003x* - 0.00020x2 


ies) 
(oS) 
ce) 
— 
x< 
—_ 
iT 


4 : g(x) = -0.07173x + 0.01787x* - 0.00161x° = 0.00004x" 


ee 


Table 4. Values of Akaike's Infomation Criterion for 
Polynomial Type Models 


Fa ee ee eee 


Model Code AIC for Each Model 
1 0.9] 
2 -16.44 
8 -15.66 
4 -13.99 


eee 


-—0G > 


Nn Staats -S- BaF See =P ee ee BS ee A VAAN a (18) 


A 


The maximum likelihood T must satisfy the condition; 
ple pve oy or sel 1 eal at ban 8 na ee ce rie Te a eC gran ee ae ae ee (19) 
ele) fe) 


From this condition and Eq.(17), 


and 


In fF = =n-iIn[v2r]o-yn-info} = A. 
20 
Therefore, Akaike's information criterion AIC for our polynomial 


type model of safety evaluation function is calculated as follows; 
Ca SS PRR oe 28 is a a ae re (22) 


where C is constant. 

The value of AIC for the above estimated polynomial type models of 
safety evaluation function are shown in Table 4. According to AIC, the 
second order model is the best among the polynomial type models. 

However, this safety evaluation function is not monotonous as seen 
in Fig.l. That is, this function increases in a range in which the loss 
is small. It implies the strange result that more loss is preferable in 
that part of the range. To avoid this strange conclusion, another type 


of model for the safety evaluation function is needed. 
3.4. Threshold Type Safety Evaluation Function 


Instead of the polynomial type of safety evaluation function, 
consider the following threshold type model of safety evaluation 


function; 


=. PRE 


m ; 
g(x) = ) a. (h(x,6))", OO ee es (23) 


where 


h(x,.0) = x - 6 if x > OF 


In this case, the safety level can be written as; 

. i 
s, = ) | asfta(x): (G0) 1a | a eee ne (25) 
k oq ey k 


and the condition to minimize Q is; 


i m : ‘ 
Pr =D ag Ph” (6) fy Oo) dx- fh (x, 8) fy , (20 dx) ] 


x [{h* (x,6) £, (x)dx-fht (x, 9) £ (XX es (26) 


We can solve the above equation system numerically and estimate the 
threshold type safety evaluation function for each order ime=» bj 2,03 
and 4 as reported in Table 5 and Bigaw2e 

The threshold type models have a new parameter 0 in addition to 


parameters a.. Akaike's information criterion is; 
AIC) = 0 sin (Q) sc 2s 2a Coe So hte ae) (27) 


The value of AIC for each threshold type model of safety evaluation 
function is calculated by the above equation and the results are shown 
in Table 6. 

According to AIC, we can conclude that the best model among the 


threshold type models is; 


@ (xX) -=)=0,10987 (30, 9.6) esas sae: eae | el pam (28) 


where 


=) al Jee 


Table 5. Estimated Threshold Type Safety Evaluation Function 


Model Code Order of Model Estimated Model 
5 1 g(x) = -0.1099h(x,9.6) 
6 2 g(x) = -0.01245h(x,7.1) - 0.00693h°(x,7.1) 
7 BN g(x) = 0.00057h(x,6.8) - 0.00951h*(x,6.8) 


+ 0.00016h?(x,6.8) 


8 n g(x) = -0.06049h(x,4.6) + 0.02483h2(x,4.6) 
- 0.00388h2(x,4.6) + 0.00014h7(x,4.6) 


where h(x,0) = x - 6 for x> 0, h(x,6) =0 for x <0. 


Table 6. Values of Akaike's Infomation Criterion for 
Threshold Type Models 


Model Code AIC for Each Model 
5 -16.83 
6 -15.83 
i -13.88 
8 -12.48 


F.g.2. Threshold type safety evaluation function 


x- 9.6 if x2 9.6 
0 Sp ee ORO eee, ee pele 2 A (29) 


h(x,9.6) 


This threshold type model of the safety evaluation function is 
better than any polynomial type model from a viewpoint of Akaike's 
theory. Moreover, this function is monotonously non-increasing. This 
property of the safety evaluation function does not contradict our 


feeling that more loss is more dangerous. 
4. CONCLUSION 


We obtained an equation (7) that shows a relation between _the 
safety preference for an urban area and the safety level defined by 
(per -4)yand Eq. (5). .iIn the case where the loss vector has only one 
component: human loss, the safety evaluation function is obtained as a 
threshold type function presented by Eq.(28) and Eq.(29) from the case 
study safety preference survey. By this result, it is possible. to 
consider that there is a threshold such that any losses perceived as 
negligible. 

The safety level of any urban, Sho can be calculated with the 
safety evaluation function. For example, if we use the safety 
evaluation function which is estimated in the case study, we can show 
the safety level S, as follows when we can obtain the frequency function 


Or loss x, f(x); for the urban area k; 


S- 8 1098 MIRO O86 dx eens eae es teeta (30) 
k 
X29.6 


These results show the possibility of obtaining the safety 
evaluation function for the case where the loss vector has 
multi-components, permitting us to assess the safety level of an urban 
area synthetically. Further research on this possibility is needed in 


order to develop safety planning methods for modern urban areas. 
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FOREWORD 


Progress of research on the behavior of building structures is 
remarkable; member and subassemblage tests have provided detailed 
information on the behavior of individual building components. 
Analytical study, such as computer simulation, has also enabled us to 
examine the complicated behavior of building structures. While these 
findings have been incorporated into our design practices, one problem 
remains as to how the results should be interpreted in view of the real 
behavior of structures. To assess the suitability of our engineering 
judgements, there is a increasing demand to test building structures on 
a large scale. | 

Building Research Institute set up in 1978 a new test laboratory 
for large scale testing of building structures, and subsequently has 
developed new systems that facilitate such large scale testing. This 
paper is the first issue of a series of reports in which we plan to 
introduce the details of the test facilities and systems in this 
laboratory. This first paper depicts an overview of the test 
laboratory, loading and measuring techniques, and computer controlled 
test system. 


We wish readers' critical comments on the content of this paper. 


Katsuro Kamimura 


Director General 
Building Research Institute 
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AT BRI LARGE SCALE STRUCTURE TEST LABORATORY 


by 
Shin Okamoto 
Takashi Kaminosono“ 
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Masayoshi Nakashima‘ 
and 
Hiroto Kato® 


Building Research Institute, Ministry of Construction 
Tsukuba, Ibaraki JAPAN 


ABSTRACT 

Study on the behavior of building structures has greatly been 
progressed. Tests of structural members and assemblies and scale models 
Of buildings have been conducted extensively all over the world. 
Numerical techniques for analysis of structures have also been developed 
a great deal thanks to innovation of digital computers. While many 
findings obtained from these studies have been reflected to various 
design practices, one problem remains as to how effectively they have 
been interpreted. Many engineering judgements are needed to incorporate 
member and assembly test data into structural design. Assumptions and 
simplifications are inevitably employed in numerical analysis. In order 
to evaluate the validity of these judgements and assumptions and to 
examine ultimately the overall behavior of building structures, one 
needs to execute tests of building structures on a large scale. 

This paper introduces the structural test laboratory at the 
Building Research Institute, Ministry of Construction, Tsukuba, Japan. 
The laboratory has various facilities specially designed for large scale 
testing. They include the large test floor, loading system, data 
Processing system and computer supervised testing system. 

1. Director, Production Department 


Research Engineer, Production Department 
Assistant Research Engineer, Production Department 


The laboratory has an open space of 30 meters by 60 meters in plan 
and 40 meters in height. In the middle of the laboratory, a 25 meter 
high wall is built to provide reaction to horizontally installed 
actuators with which earthquake forces are to be simulated. The 
Capacity of the wall is 40 MN and 720 MN x m in allowable shear and 
bending moment. Either of two test floors, one on each side of the 


é floor 


wall, can accommodate a building specimen of as much as 300 m 
area and 25 meter height. 

Servo-controlled actuators are prepared for load application, 
having capacity of + 1000 KN in force and + 1000 mm jin stroke. The 
needed power is provided by hydraulic pressure, and the actuators are 
controlled by-electronic servo-controllers. 

Various types of displacement transducers (DLT's) are provided for 
the measurement of large displacements. Specially designed digital type 
DLT's are unique. They are magnetic response type and can measure 
displacements up to + 1000 mm with accuracy of 0.01 mm. Signals of 
these DLT's are completely noise free. 

For systematic loading, data processing, and safety check, two 
computers are introduced. One serves for servo-control. During the 
test, the computer keeps supervising the motion of actuators for 
accurate command control and for safety check. The other computer, on 
line with data acquisition units, is used for data processing.» Jhis 
computer is linked to the  servo-control computer. Commands’ for 
servo-controllers can then be created in this data processing computer 
and sent to the servo-control computer. 

These two computers, servo-controllers, actuators, measuring 
instruments, and data acquisition units constitute one closed loop test 
System. This closed loop system, involving minimized manual Operation, 
enables us to perform a new test method developed for simulation of 
building response under earthquake motions, referred to as_ the 
Ppseudo-dynamic (PSD) test method. This paper also introduces large 
scale tests conducted at the laboratory by means of this new method. 
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I INTRODUCTION 


Research on building structures has been progressed for decades, 
and many findings obtained from works associated with the behavior of 
various types of structures have effectively been incorporated into 
design practices. Two approaches are common in those research works. 
Experimental work is the first approach, where scaled models of 
structural members, subassemblages, and frames are tested with various 
loading and support conditions, and these results are correlated with 
the behavior of real building structures. The other is the analytical 
approach. This approach ranges from purely mathematical formulation of 
physical phenomena (belonging to the applied mechanics) to numerical 
analyses using large size digital computers. 


In these days, experiments have been facilitated a great deal 
thanks to the innovation of technologies of various kinds. Load 
applying actuators have been upgraded in both capacity and accuracy. 
Dynamic servo-controlled actuators are now available in many structural 
laboratories. Very fine measuring instruments make it possible to 
monitor displacements, rotations, and strains of test structures with 
great accuracy. The actuator control, and storing and reducing data 
have also been advanced by means of digital computers. The progress of 
analytical study is also remarkable. Specifically, work that utilizes 
large size digital computers has proven its great capability. Iterative 
solution techniques enable to solve various non-linear problems. The 
finite element method is one of those powerful techniques. 


With all of these sophistications in experiment and analysis, there 
is one thing still lacking before making final judgement on the 
evaluation of the behavior of real building structures. That is the 
test of building structures as building structures. Component tests 
provide useful information on the individual characteristics of these 
members, but by these tests, little can be known with respect to the 
Overall behavior of building structures in which many members are 
combined. Structural integrity, in other words, the effect of member 
connections on the overall behavior, is after all difficult to predict 


by the member tests. The Subassemblage test is a useful approach to 


investigate closely the behavior of a building structure as a unit. 
There is still one problem; when selecting an subassemblage, we must 
assign some artificial boundary conditions in the subassemblage. In a 
real structure, the part that is extracted as a subassemblage is 
Supported by other members and subassemblages. In many occasions, it is 
difficult to perfectly simulate those boundary conditions. The ultimate 
validity of the adopted (or assumed) boundary conditions can only be 
checked by comparing with the behavior of the real structure. 
Structural tests of scaled models are often employed in order to 
evaluate the overall behavior of building structures. Tests of scaled 
models always accompany one problem as to how individual constituents of 
real structures are to be scaled down. For example, in reinforced 
concrete structural models, the size of aggregates and reinforcing bars 
cannot be scaled down correctly once the scale ratio is lowered to a 
certain level. 


Numerical analysis such as the finite element analysis is a very 
powerful tool to study the overall behavior of building Structures, and, 
as a matter of fact, many practical design provisions, Suggestions, and 
recommendations have been made on the basis of these analytical results. 
Nevertheless, we must be reminded that assumptions of various levels - 
stress-strain level, moment-curvature level, end force-end displacement 
level, or story level - are inevitably involved in these analyses. 
These assumptions can first be verified when the analytical results are 
compared with the behavior of real structures. 


All of those statements suggest a need of executing tests of full 
(or at least large) scale building structures. It does not indicate, 
however, that we need to test large scale structures all the time. Such 
tests are very costly, and therefore need only to be done on a timely 
basis when final confirmation is required. 


The Building Research Institute of Ministry of Construction 
constructed in 1978 a laboratory for tests of large scale building 
Structures, and subsequently has developed various systems that 
facilitate such testing. In addition to conventional test equipment, 
special facilities are indispensable for large scale tests; first. 


need a large open space to accommodate a large scale building specimen. 


Actuators that apply forces and deflections which are commensurate with 
the capacity of the test building are required. Not only accuracy but 
large capacity in stroke is also required for some measuring 
instruments. Since we deal with a large scale specimen, many 
displacement, rotation, and strain signals must effectively be read, 
stored, and reduced into meaningful forms. 


This paper reports on this large scale structural laboratory and 
the outline of the techniques developed for large scale testing. A 
number of large scale tests performed in this laboratory are also 
described. 


II LABORATORY FOR LARGE SCALE TESTING 


In order to execute large scale tests, first we need a good 
environment for the testing. The test laboratory shelters a Space which 
is about 30 meters by 60 meters in plan and 40 meters in height (Fig. 
1). Figures 2 and 3 show the overall dimension of the test laboratory. 
In most tests, specimens must be rigidly supported on the test floor. 
The test floor, therefore, is designed so that it is Strong as well as 
Stiff enough to sustain securely the forces applied through the 
Specimens. The test floor, made of prestressed concrete, has a 
dimension of 44.6 meters by 20 meters in plan and 6.6 meters in 
thickness and is able to sustain 1.0 MN/m? of bearing force and 720 MN x 
m of bending moment (Fig. 4). 


To simulate wind or earthquake forces applied to a Specimen, it is 
often needed to install load applying actuators in a_ horizontal 
direction. For this purpose, a wall, referred to as reaction wall, is 
arranged in the middle of the test floor (Frqwes)ievathe: wall "4's? 20 
meters by 25 meters in surface dimension and 6.6 meter thick. The wall 
is also made of prestressed concrete, having an allowable capacity of 
40 MN and 720 MN x m in shear and bending moment respectively. Having 
the wall placed in the middle of the floor, the test floor is 
essentially divided into two areas (Floors A and B in Pid tags mi nese 
floors are respectively 20 meters by 24.6 meters, and 15.4 meters by 20 
meters in plan dimension. Therefore, the larger test floor (Floor A) 


can easily accommodate a building that is as much as 300 mé in floor 
area and 25 meters (approximately eight story equivalent) in height. 
Since the test floor is divided into two areas, it is possible to 
execute many tests simultaneously. 


Test specimens, load applying actuators, and other test setups must 
be securely tightened to the test floor as well as the reaction wall. 
For this tightening, the floor and wall have holes that are 100 mm in 
diameter and spaced at every 500 mm pitch in both directions (Fig. 6). 
Some holes penetrate through the thickness of the floor or wall, whereas 
the other holes are threaded and stop at one meter in depth. By using 
specially designed prestressing steel rods, specimens, actuators, and 
setups can be clamped to the floor or wall. 


To make facilitated connecting cables extending from various 
measuring instruments to data acquisition units in the measuring room 
(Fig. 3), permanent cables are embedded throughout in the floor and 
wall. These cables extend from the data acquisition units and 
terminates at cable relaying boxes, called “Junction Boxes", that are 
installed at strategic locations in the floor and wall. Cables 
extending from measuring instruments then can readily be connected to 
junction boxes in the vicinity of the test area, and the Signals are 
sent forward to the data acquisition units. 


III LOADING SYSTEM 


For large scale testing, load applying actuators must have Capacity 
sufficient to provide large force to a test specimen and at the same 
time to follow the specimen's deflection. The most suitable power 
element to generate such large force is hydraulic. All of the actuators 
for the large scale testing then are controlled by an electro-hydraulic 
servo-control system. 


Two types of actuators are available in the laboratory. The first 
type of actuators, designated as JB type, has a capacity of + 1.0 MN in 
maximum applying force and +1000 mm or +500 mm in full stroke (Fig. 7). 
In each of these JB type actuators, one load cell and one inductive type 
displacement transducer (DLT) are attached. Signals from the load cell 


or DLT are first sent to one of the junction boxes and further to the 
data acquisition unit, or to the servo-control unit if those signals are 
to be used for test control. The maximum piston velocity of these 
actuators are 2 mm/sec., and, therefore, they are primarily used for 
quasi-static loading tests. 


Another type of actuator, designated as JA type, is able to run to 
+ 300 mm with as much as 350 mm/sec. in piston velocity and primarily 
used for vibration and fatigue tests (Fig. 8). The maximum force that 
could be applied by this actuator ranges from +500 kN to +750 — kN in 
accordance with the piston velocity of from 350 mm/sec. down to the 
quasi-static level. One load cell, two DLT's, and one accelerometer are 
installed in this actuator. One of the DLT's covers the full range 
(+300 mm), while the other measures within a range of +60 mm but with 
higher accuracy. As in JB type actuators, signals of the load cell, 
DLT's, and accelerometer in this actuator are sent to the data 
acquisition unit and to the servo-control unit. A function generator is 
prepared for this actuator. The generator can provide sinusoidal, 
Square, ramp, sawtooth functions with a frequency range of between 0.01 
to 10,000 Hz. 


The actuators have a swivel base and swivel head to provide freedom 
at both ends of the actuators. The type, number, and capacity of these 
actuators are tabulated in Table 1. 


Hydraulic power units provide the needed power for the actuators. 
Six pairs of oi] reservoir and hydraulic pump are installed at the 
basement of the laboratory (Fig. 9). Each of the pumps is capable of 
delivering 0.3 m>/min. of oi] at a pressure of 21 MN/m<. 


The motion of an actuator piston is controlled by the hydraulic 
Servo-control system. The system consists of a feedback selector panel, 
input module, servo controller, and servo-valve. This system together 
with the actuator forms a closed-loop. Input, given as either 
displacement or load, is applied to the input module manually or from 
various external sources such as the function generator, magnetic type 
reader, or computer that are on-line with the servo-control System. 
(The computer on-line testing system will be discussed in detail in a 


later section.) An input signal, either displacement or load increment, 


“+ 


is converted to a voltage change in the servo-controller (Fig. 10) (the 
full range corresponds to +10 volt). This electrical command signal is 
then converted by means of the servo-valve to a regulated flow of high 
pressure hydraulic fluid. 


IV MEASURING TECHNIQUES 


In a quasi-static structural test, displacement and rotation 
measurement is of the most importance to evaluate the structural 
behavior of the test specimen. Various types of displacement 
transducers (DLT's) are prepared in the test laboratory. A number of 
DLISs. that can read small to medium displacements (a range up to about 
+100 mm) are available. They are either strain gage type or inductive 
type. In addition to those DLT's, DLT's that can measure large 
displacements are inevitable for large scale testing. For this purpose, 
two types of large range DLT's are provided in the laboratory. The 
first type of DLT's is inductive type and has a capacity of either +200 
mm, + 300 mm, or +500 mm in full stroke (Fig.- 11). The second is 
digital type with a full stroke of either +500 mm or +1000 mm (Fig. 12). 
In each of these digital type DLT's, plus and minus magnetic elements 
are coated alternately at a pitch of 0.01 mm along its length. A cursor 
of the DLT shifts in accordance with the movement of a measuring point, 
and the displacement is measured by counting the number of magnetic 
elements that are passed by the cursor. These DLT's are digital type 
and free from any noise signals, which are often generated in analog 
type DLT's that need amplifiers. The detailed description of the DLT's, 
including the available number, stroke capacity, and sensitivity, are 
Shown in Table 2. 


DLT and strain gage signals are first transmitted to scanning boxes 
set on the test floor. Cables extending from those scanning boxes are 
to be connected to junction boxes (described in the previous section). 
The signals, scanned with 0.1 sec. per channel, are then sent forward to 
the static data acquisition unit in the measuring room (Fig. 13). At 
the present time, a total of 700 channels can be connected to those 
Scanning boxes (seven boxes with each box having 100 channels). 


In the data acquisition unit, analog signals are amplified and 
converted to digital values by an A/D converter. Those digitized data 
are stored in a cartridge disk of a data processing computer which is 
connected to the unit. Data reduction can be done with the computer, or 
alternatively, those data can be sent to a larger computer in the 
computer center for more refined data reduction. 


All of inductive type DLT's and many of strain gage type DLT's can 
also be used for measurement of dynamic tests. Signals of DLT's and 
accelerometers recorded in a dynamic test are assembled to bridge boxes 
on the test floor and sent forward to the interface of a dynamic 
measuring unit in the measuring room. Then the signals are scanned by a 
multiplexer in this unit with a speed of 0.000050 sec. per channel and 
digitized by an A/D converter. Those digitized data are stored in a 
cartridge disk as for the measurement in static tests. Analog signals 
can also be recorded directly in magnetic tapes before passing the A/D 
converter. Table 3 provides details on the capacity of the data 
acquisition units. 


Since the test control room is away from the test floors, and 
furthermore it is very difficult to keep tracking all local behaviors of 
a large scale test structure, many portable TV cameras are provided in 
the laboratory for the purpose of monitoring those local behaviors. 
These cameras can be mounted at strategic positions so that critical 
locations of the structure are to be continuously monitored. Pictures 
sent from those cameras can be displayed in color displays in the 
measuring room (Fig. 14). Those pictures can also be recorded in video 
tape recorders if it is found necessary. 


V COMPUTER ON-LINE TESTING SYSTEM 


When a test specimen is large, the level of applying force becomes 
proportionally high. Many load, displacement, and strain data also need 
to be acquired during the test. More importantly, it should be noted 
that malfunction in testing may cause a disaster since the Specimen no 
longer is a miniature. Under these conditions in large scale testing, 


an unified control of loading and data processing as well as continuous 


safety check are very critical. To achieve effective execution of large 
scale tests, the overall testing has been designed to be computer 
controlled, named the computer on-line testing. This computer on-line 
testing system (Fig. 15) comprises one servo-controller control 
computer, servo-controllers, servo-control led actuators, data 
acquisition units, and one data processing computer. All of them, 
linked together, constitutes one closed loop as shown in Fig. Il. The 
servo-controllers, actuators, and data acquisition units have already 
been described in the previous sections. 


The servo-controller control computer, as the name explains what it 
is, controls’ the servo-controllers of actuators. The loop of test 
control - control by actuator loads, by actuator displacements, or by 
external measuring instruments - can be selected by this computer. The 
piston velocity of actuators can also be set up from this computer. If 
two or more actuators are used in one test, the piston velocity of each 
actuator can be adjusted so that all actuators would reach their 
specified command values simultaneously. Accuracy of the attained 
command values (command displacements, loads, or external signals) 
relative to the specified command values is continuously checked by this 
computer. First, allowable error bounds are set up. The computer then 
keeps monitoring the control command values as to whether the measured 
values are within the specified error bounds. When all values fall 
within these bounds, the actuators stop their motion, and the computer 
Starts a stability check. The computer scans the command values at a 
specified time interval (of an order of 0.00005 sec.) and repeats the 
scanning for a certain period (say, for one second). If those values 
constantly are within the error bounds, the computer presumes that the 
Specified values are all achieved and sends a message to the data 
acquisition units, for starting measuring. If the values exceed the 
bounds during the scanning, the stability is assumed not to be achieved, 
and the computer gives a pause signal to the servo-controllers. Several 
options are prepared in this situation; the stability check may be 
repeated again, the error bounds may be enlarged, or the test may be 
continued as if the stability check were cleared. 


The computer keeps monitoring the signals of actuator loads and 


displacements in order to ensure the safety of the test. When one of 
them exceeds its predetermined (upper or lower) limit, a command is 
given to stop the loading. The computer also monitors the hydraulic 
pressure of oil, the currencies and voltages of the servo-control 
interfaces, servo-valves and others. Upon detecting an abnormal signal, 
the computer gives a pause signal to the servo-controllers. 


The data processing computer are on-line with on one side the data 
acquisition units and on the other the servo-controller. control 
computer. This data processing computer is composed of its main body, 
One system typewriter, one card reader, one line printer, one X-Y 
plotter, and one graphic display. The computer receives from the data 
acquisition units key data, which are selected by a fortran computer 
program made by a test operator. Data reduction is performed according 
to the fortran program, and the reduced data can be output on the output 
devices. The displayed pictures can be copied in hard sheets if 
requested. Sometimes, the test (loading) program may not’ be 
predetermined; in other words, command signals are to be determined on 
the basis of data acquired in previous steps of the test. This computer 
is capable of executing such post-determined program test. The computer 
receives data needed for computing the next command signals, executes a 
fortran program that is to compute the signals, and send the Signals to 
the servo-controller control computer. A flow chart of this computer 
Supervised closed loop testing system is given in Fig. 16. 


VI PSEUDO-DYNAMIC TESTING METHOD 


The Pseudo-Dynamic (PSD) testing method is a technique to simulate 
the earthquake behavior of structures and was originally proposed by 
Tanaka et al. (Ref. 1). The basic concept of this method is to solve 
equations of motions of a structure by means of the direct integration 
technique. This method, however, is not a pure numerical analysis but a 
combined simulation technique of experiment and numerical analysis. To 
explain this method, let us consider an example shown in Fig. 17. 
First, the test structure is set on a test floor, and actuators are 
installed at appropriate positions (Fig. 17 (a)). The test structure 


then is assumed to be a discrete mass system (Fig. 17 (b)). The 
equations of motion of the system can be written in a matrix form as: 
CS eee ee re nee eee Sante 
Where {Xn}, {Xn}, and {Xn} respectively indicate the acceleration, 
velocity, and displacement vector of the discrete system at the n-th 
step. [M] and [C] are the mass and viscous damping matrix. {fn} denotes 
the restoring force vector of the system at the n-th step. 


Various integration techniques have been developed to solve the 
equations of motion for linear and nonlinear structures. The techniques 
are divided primarily into two groups. The first group, referred to as 
the implicit solution technique, includes the techniques such as the 
Houbolt method, the Newmark method, and the Wilson method. The second 
group, referred to as the explicit solution technique, includes the 
central difference method, the two-cycle interation with the trapezoidal 
rule, and the Runge-Kutta method. Among those techniques, the central 
difference technique is used to solve the equations in the PSD testing 
method. In the central difference technique,- the acceleration and 
velocity vectors at the n-th step are expressed in terms of the 
displacements at the (n-1)-th, n-th and (n+1)-th steps as: 

{Xn} __{xnt1} - 2{xn} + {Xn-1} | wo ech pee ee sale ie | ee oe Eq. 2 
[Ni 
ete {xn+1} - {Xn-1} 
Zee 
By substituting the above two equations into the equations of motion 
(Eq. 1), and solving for the displacement vector at the (n+1)-th step, 
we can obtain the following relationship: 


(xne1} = (Om) +S (cI) *L20M}om} + (5 Cc] - (MI) m2} 


wieeletigc cs 
+ (At)*{fn} - [M] (At)? oXn] 


That is, the displacement vector at the (n+l)-th step can be derived 
from the displacement vectors at the (n-1)-th and n-th steps, the ground 
acceleration at the n-th step, and the restoring force vector at the 
n-th step. The PSD test procedure of the structure shown in Fig. 17 can 
be described as follows: 


Suppose the test specimen is in the n-th step, the n-th step ground 
acceleration oXn, the displacements at the (n-1)-th and n-th steps 
{Xn-1} and {Xn}, and the restoring forces (measured from actuator 
load cells) {fn}, are all known. Substituting these values into 
the equations of motion (Eq. 4), we can solve the displacements at 
the (n+]l)-th step {Xnt+1}. Next, the specimen is loaded 
quasi-statically to the calculated displacements by the actuators. 
When the displacements are attained, the restoring forces {fntl}, 
are measured by the actuator load cells. The test is advanced by 
repeating this procedure. A flow diagram of the test is also shown 
Tang. 618 


Some good reasons to employ the central difference technique can be 
drawn in view of the above formulation and procedure. *First,” in this 
technique, no iterative process is required even in the inelastic range. 
It should be noted that, in an implicit technique, iteration is 
inevitable in the inelastic range because the equations of motion are to 
be formulated for the (n+1)-th step rather than for the n-th step. Such 
iteration is not suitable for the PSD method. Second, as seen from Eq. 
4, the computation does not involve any matrix calculation if the mass 
and viscous damping matrices are diagonal, which we generally adopt in 
PSD tests. As a result, the time and Storage space required for 
computation is reduced significantly. According to previous studies 
(Refs. 2 and 3), the central difference technique is stable and 
reasonably accurate as long as the time interval for the calculation is 
smal]. The stability condition is given as (Ref. 3): 


At < Tin aide: GAs tSIGE Sk ARIS Ae cas BIRT a es 


Where At is the time interval, and Tmin is the smallest natural period 
of the mass system representing the test structure. Further, use of 
One-sixth of the above time interval is recommended to guarantee 
sufficient accuracy. 


The most prominent advantage of this method over numerical analysis 
is that no assumption is required for the stiffness (restoring force) 
Characteristics of the structure. The PSD method, therefore, is very 


powerful when analyzing dynamic behavior of inelastic systems. A second 


* 


benefit of this method is that it is possible to keep tracking the local 
behavior and damage propagation of the test structure since the test 
structure is loaded quasi-statically. This statement can readily be 
understood if the test is compared with a shaking table test, in which 
the testing time is very short (of an order of 10 to 30 seconds). In 
the PSD test, selection of the viscous damping matrix [C], however, is 
debatable as in numerical analysis. Idealization needs to be made such 
that [C] be proportional to the initial stiffness matrix, to the mass 
matrix or to the combination of the two. Such idealization can be a 
disadvantage if compared with shaking table tests. It may be cited, 
however, that miniature model structures used in a majority of shaking 
table tests are also poor in simulating the damping characteristics of 
their prototypes. 


To execute a PSD test, the restoring forces of the test structure 
must be measured, and are to be sent to a computer which calculates the 
equations of motion. Then, the calculated displacement increments for 
the next step must be applied to the test structure. It is clear that 
the computer supervised closed loop test system described in the 
previous section can control the execution of the PSD test 
systematically. Furthermore, it should be emphasized that the 
earthquake behavior of a structure which is too large to be tested on a 
Shaking table can directly be simulated with this PSD method, for the 
laboratory of the Building Research Institute can serve for a large 
scale test structure. To the authors' best knowledge, the world largest 
Shaking table, a counterpart of this laboratory with the PSD test 
method, is 15 meters by 15 meters in table dimension and can carry at 
most 10 MN of weight. From these numbers, we can notice a sizable 
difference in the allowable maximum scale of test structures between the 
Shaking table test and the PSD test. 


VII LARGE SCALE TESTS CONDUCTED AT LABORATORY 
Table 4 lists large scale structural tests conducted at the 


laboratory from 1978 through 1982. As seen from the third rightest 
column of the table (Objectives), a majority of the tests was made to 


ees el 


investigate the earthquake behavior of various types of building 
structures. 


VIII CONCLUDING REMARKS 


This paper has introduced various facilities specially designed for 
the large scale testing at the Building Research Institute of Ministry 
of Construction. We demonstrated distinguished features of those 
facilities particularly those closely associated with the performance of 
large scale tests. For effective and systematic control of loading, 
data processing, and safety monitoring during the test, a computer 
Supervised closed loop system has been devised and explained in detail 
in this paper. We believe that this large scale testing system has high 
potential for research of various aspects of structures. Some of the 
research projects being planned are as follows: 

i Study on the behavior of high-rise buildings: There have been very 
few experimental works to have investigated the wind and earthquake 
resistance mechanism of high-rise buildings. Since the laboratory 
has a space that can accommodate a specimen of as high vas, 25 
meters, a test specimen can rise to 25 story high if one third of 
scale ratio is adopted. 

2. Study on the behavior of three-dimensional buildings: itecd 
building that has a slender rectangular plan Shape; that is, has 
many spans in one horizontal direction relative to the other 
direction (like an apartment building), the distribution of 
earthquake force to each force resisting component is very complex 
particularly if the force is applied beyond the elastic limit of 
the building. With one third of scale ratio, a building model that 
has aS many as ten spans can be tested in this laboratory. 

3. Study of soil structure interaction: Behavior of foundations, 
footings and piles is still very difficult to evaluate because of 
the complexity of soil properties and the interaction between the 
Super structure and the soil. Considering that soils are difficult 
to scale down properly, a full scale test can be performed on the 


laboratory test floor of 400 me, 


NOTATION 


The following symbols are used in this paper. 


Rey : Viscous damping matrix 

{fini : Restoring force vector at the n-th step 

[M] : Mass matrix 

{Xn} : Acceleration vector at the n-th step 

{xn} : Velocity vector at the n-th step 

{xn} : Relative displacement vector at the n-th step 
oXn : Ground acceleration at the n-th step 

Tmin : Smallest natural period of discrete mass system 
At : Time interval for direct integration 


Table 1 Capacity of Actuators 


ACTUATORS 


Number of |Maximum Maximum| Maximum 
Actuators Ram 


Velocity 
(mm/sec. ) 


+1000 


+500 Re 350-0 +750 
(dynamic) 

+750 

(static) 


Table 2 Capacity of Displacement Transducers 


DISPLACEMENT TRANSDUCERS (DLT) 


Type Number of| Full Sensitivity 
DLT ea 
(mm) 


Strain Gage +300 0.1 mm 


wp - 
(2) 42) 


+100 0.01 mm 


ine) 
(=) 


ou ae 0.002 mm 


ine) 
Sy 


Inductive 


2 
fa 0.5 percent of 
A full range 
2 
Digital 3 +1000 0.01 mm 
5 
3 0.001 mm 
2 
Ms i esa 
] 


Table 3 Capacity of Data Acquisition Units 


DATA ACQUISITION SYSTEM 


1. Static Strain Data Acquisition System 


Number of channels: 


Number of scanning boxes: 


Number of strain measuring units: 2 


Scanning time: 0.1 sec./channel 


Measuring range: +30,000 x ite u 


2. Dynamic Strain Data Acquisition System 


a. for digitized data 
Number of channels: 30 
Number of bridge boxes: 3 
Number of signal conditioners: 30 
Number of low path filters: 30 
Scanning time: 50 u sec./channel 
Measuring range: +30,000 x 107° u 


Maximum frequency measurement: 5,000 Hz 
b. for analog data 


Number of data recorders: 14 


2861 “AON 0} anbLuydsa|] Sal peo) auweu4 padeug 

al Z86L ABW] Sd JO UoLzeN|eAZ/ayxenbyzUeJ [aais Auois Gg | 2l 
oWweU4 do owed 4 
2961 “Bny 07} -9U0) passaizSedg|peo] *Luo} 91949U0) passau 
vl Z86L Ltudy|Jo woLraeyag OLuistas| 913e9S -3Sadd Ku0ys 2 


sa0eug-» 
YZLM oWe44 [9294S 
40 abelquasseqns 


au 4 

[9981S UL Sadeug-> 
40 AOLARYAG ILWSLAS 
[ [PMAeOUS 


YZLM BULpLing Oy 
$0 AOLARYSG ILWSLIS 


peo] * 140} 
343045 


Z86L “AON 0} 
2861 Lludy 


BuLp|ing 


c86L “uel 03 
dy A409s L 


O86. °3das 


aie 086, aunt 03 anbLuyda| 3Sal peo] 
LL 0861 °99G| GSd JO UoLzENLeAR |ayenbyqe4y aueuy Jy f401S pb | 8 
ot O86L “bny 07 anbiuyde| Sa] peo 
Ol Og6L ALNC| dSd JO YOLZeN|eAR|ayxenbyzueJ awerjs Jy A409S Z l 
eee O86L “qed 0} awes4 OY peo] awed 4 
6/61 390/40 woLAeyag OLuIsLas |ayxenbyzue4; AeURLd JY AW0IS Z@ | 9 
LLM LLM 
are O86L YOuRW 02} -4eaYs paldnoy JY} peo] *L4OH -ueays paldnod dy 
6/61 °99G|/4O UOLAeYsG ILWSLIS 9L3e45S 40 abel quasseqns G 


BULPLENg [aued 


O86L YOURW OF)  pazedluqesaid JY}PeO] *L4OH Butpling Laued 
l 6/61 aune| Jo uoLaeyag Jtwstas| 919295 payeotuqejaud OY | 
6/61 °30eS 03 anbiuyda| 4Sal aweA} paodeug ree 
9 6/61 Liddy]. dSd $0 uoLzen|eAR|ayenbyqued | €/L [aaas AuoIs 2 | € 
BuLpling 
6/61 ‘da4 03] payedtuqesaud Jy;peo] “L40H {LLNS (A40Is ¢) BuLpLing 
5 9/6, “Sny|4o woLAeyag JLwstas| 917895 | payediugesaud OY | 2 


wazshS ssnul de 


peo] “eA 
40 abel quasseqns 


3L324S 


aweu4 a0eds 40 
Ay.oede) 93} ewl3iN 


6Z6L “Uet OF 
8/6. ‘3das 


Rou URI d 


# POLuad Burpeoy jolzey| 2YHLeH 
UawLdads 


"Jay 189] SaAL IeLqO 40 adh] jayeos| uawLdeds 


aunzonuys 40 adh] 


Ksozpeuoge] 7e pajyonpuoy) szsaj ales abue] yp a1 qel 


Fig. 1 Large Scale Test Laboratory of Building Reseach 
Institute, Ministory of Construction 
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Fig. 2 Elevation of Laboratory 
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Fig. 3 Plan of Laboratory 


i 


Fig. 4 Floor and Wall Load Carrying Capacity 


Fig. 5 Reaction Wall (25m x 20m x 6.6m) 
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Fig. 6 Dimension of Test Floor and Reaction Wall 
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Fig. 7 JB Type Actuator (for Static Loading) 


Fig. 8 JA Type Actuator (for Dynamic Loading) 


Fig. 9 Hydraulic Pump Unit at the Basement of Laboratory 
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Fig. 10 Servo - Controllers in Control Room 


Fig. 11 Analog Type Displacement Transducers 
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Fig. 12 Digital Type Displacement Transducer 
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Fig. 13 Data Processing System 


Control board 


Test Monitoring Unit (Displays and Video Recorders) 


wazshs Bburysa}, doo P8SOL) pastauadns uazndwoy Giecbiy 


gees 
LINN NOILISINODY AlddNS Y3MO0d JITNVYCAH 
viVd JILVLS 


. TANVd 10YLNO9 
OULNOD OAYSS \ <7 , 
ra OT TNVYCAH ie 


TOYLNOD OAYSS 


a: IAA 
. cee gil = Kee 


NOTLONAF 


LIND NOTLISINOOY ae 4 
ae Sie ) 
VIVG DIWWNAC , L : eee ce 
: ON S<Seer 91 INVYCAH-0Y 19.373 


X 


O@ ININNVIS 


STANVd 


f- SAWIdSIC 3 
XO" SOACIA “SYOLTNOL 
NOTLONAE 


VYAWWD AL 
| 


fie] 
yorvnioy [C— 


YSONCSNVY L 
LNAWFIVIdS Id 


Buljsay auly-ug uwayndwo) jo weuberq mol4 91 ‘B14 


YaLNdWOD = TONLNOD — 
YITIONLNOI-OANSS : II - Add 
YdLNdwod 

QNISSISOUd ViWG = I - Ndd 


X0d 
NOT LIANNOD 


YANNWIS 


“dW “dWV 
s,11d s, 11d 
JOTWNY WLI 


LINN 
| ee NOTLISINOIY 
NOILOIANNOD OAYAS oe 


LINN 
YaMOd 
OT INVYCAH 


J10SNO9 
TOYLNOSD 
YOLVALSV 


J1OSNOI 
TOULNOD 
OT INWYGAH 


AW1dSTG 
Ya LOVYWHO 


LINA 
JIVANT 
YdLNdWOd 


II - dd ie) ide 


YaLTYMadAL 
WALSAS 


AW1dSIG 
JIHdWud 


ASIC 
JIG YLYya 


YALNTYd 
ANT] 
udavad 
auv9 


YALTUMadAL 
WALSAS 


Sale Paes 


Wa 4SAS 
sseW quapeAtnb 


Bulysal 


(q) 


otueutg - opnasd JO udLze}UaSauday /| “Bi 


aingonsaSesole (ey 


f 1 
——— 
TW be —— | 
a 


UsWLOedS 407 PNY 


|o)| 


"and responce analysis 


i Gomme System for data acquisition 


Input of earthquake ground 
motion x 


Calculation of displacments_ for 
next step {Xnerf 


(MX } +(C) Pen} +fFa}= LO) [Keonf 
I 


me. 
| Transfer of displacements 


Reactions Wee and 
Xnsi} 


displacements 


Actuator Control System 


computer system 


2. servo controller 


; Xnei 
\ Fr+1 
NB eh 
\ { | 
: i 
1X +1 
+ 


fe 
tr 

| 

| 

| 


data acquisition system 


Fig. 18 Flow Diagram of Pseudo - Dynamic Test Operation 
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FOREWORD 


This paper describes the result of the research on the calculation 
procedure for room air temperatures including effect of natural 
ventilation. This research had been made as part of a development project 
conducted by the Ministry of Construction for 1977-1981, entitled, 
"Development of Total Energy Conservation System for Dwellings". 

In Japan, energy conservation for dwellings is one of the major 
national problems. The method of estimating heating and cooling loads in a 
dwelling house can be regarded as one of useful fundamental tools for this 
problem. 

fs this paper, taking “account of the structure and the use of 
dwelling houses in Japan, the procedure for calculating room air 
temperatures is formulated where natural ventilation takes place. Validity 
of this procedure has been checked by comparing the calculated 
temperatures with the observed ones. 

The computer program of heat loads made according to _ this 
calculation procedure was named BRIMAP and has been used to develop a 
design system and a simplified evaluation system for energy conservation 
in dwellings. I would like to express my hope for the development of 
research on this field and acknowledge with appreciation those who have 


supported on this study. 


October, 1983 
K. Kamimura, Director General 


Building Research Institute 
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DEVELOPMENT OF PROCEDURE FOR CALCULATING ROOM AIR TEMPERATURES 
INCLUDING EFFECT OF NATURAL VENTILATION * 


by 
Yuzo SAKAMOTO ** 
BUILDING RESEARCH INSTITUTE 
MINISTRY OF CONSTRUCTION 


Tsukuba, Ibaraki, Japan 


ABSTRACT 


Natural ventilation has an important effect on thermal behavior of a 
dwelling house in Japan. A calculation procedure was developed to predict 
the room air temperature and the heat load of Japanese housing. It 
includes the process to calculate natural ventilation precisely and 
quickly. The room air tempertures of two experiment houses were predicted 
by this procedure. The prediction was compared with the observation to 
valuate the accuracy and the validity of this procedure. The agreement 


between both is good and satisfactory. 


*This paper was reproduced from the manuscript prepared for the Fourth 
International Symposium on the Use of Computers for Environmental 


Engineering Related to Buildings, 1983, Tokyo. 


* : Nata ; ee ee: 
Chief Reseach Member, Building Services Division 
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1. INTRODUCTION 


Generally speaking, the dwelling house of Japan is air-conditioned 
partially and intermittently. It is very rare that rooms’ are 
air-conditioned all day long. Air-conditioning is usually performed only 
in the room used. Furthermore, natural ventilation (including infiltration 
in this paper) has an important effect on thermal behavior of Japanese 
housing. The typical example is seen for a wooden house built by the 
conventional construction method of Japan. Namely, because its walls, 
floors, windows, and other components have many cracks, the heat loss by 
infiltration is considerable. In addition, windows and doors are opened 
during cool summer nights, even if an air condition system is installed. 
This paper describes a calculation procedure developed to predict room air 
temperatures or room heat loads in the dwelling house of Japan, which has 
the features mentioned above. The actual examples of prediction by this 
procedure are shown and compared with the observation to check accuracy 
and validity. 

Air temperatures or heat loads of multi-room houses can be 
calculated by solving the system of equations of room heat balance. It is 
possible to formulate the multi-room problem using the response factor 
method!) This formulation and its solution have been presented by Eguchi2), 
The calculation procedure developed in this work is based on the same 
formulation as his, but includes the precise calculation process for 
natural ventilation and for thermal interaction with the ground on new 
assumptions and approximations. 

This calculation procedure emphasizes the importance of natural 
ventilation in the calculation of thermal behavior in dwelling houses. 
Natural ventilation has tended to be assumed as a steady quantity in these 
problems. However, in this procedure, it is treated as unsteady and the 
precise calculation is performed for its effect. Volume flow rate of air 
by natural ventilation is calculated by solving the system of equations of 
the ventilation network constituted on the basis of the crack method. An 


approximate method is adopted to solve this system of equations quickly. 


2. MATHEMATICAL PROCEDURE 


2-1 PROCEDURE FOR CALCULATING ROOM AIR TEMPERATURES OR HEAT LOADS 


(1) Room Heat Balance Equation 


Various heat transfer processes take plase inside a room. If the heat 


exchanges by radiation among the surfaces 


neglected, 


enclosing the room can be 


the heat gain or loss of room air can be divided into five 


types as shown in Fig.1. When the number of the present room is "i", these 


types are described as follows: 


255° 


7. 
13 


net heat gain by transmission through components located between i 
snd j, such as wall and a floor 

net heat gain by ventilation and infiltration between i and j 

net heat loss by transmission to the component fixed on the ground 
net heat loss by storage of heavy objects, such as furniture and a 
column 

heat gain by convection from appliances, human bodies, and so 


fort: 


ae "S#  room-i room-j 


solar radiation ( 


Fig.l Five Types of Heat Gain or Loss in a Room 


A heat balance equation for air of i is given by, 


M 
n n-1 n 
CR etait oniite(ges +onty) eae! seg? tet (1) 
ask ce ‘ay ij ey ui i i 

j=1 

Where 

o. = air temperature of i 

C5 = heat storage capacity of air volume of i 


M = number of rooms including outdoor spaces, and the superscript n 
indicates time nAt ( At is time increment) . 


The net heat gain by ventilation and infiltration is represented, 


n 1G go gee a non 
ne = Ge - era (2) 
where 


Vas = volume flow rate from j into i by ventilation and infiltration 
cy = heat storage capacity of air per unit volume. 


According to the response factor method, 


n co m n-m ig layg m n-m igh igl n 
= » ee (On ee Bae OO + Je ie +a. ] 
Qi, TA ay [ {Ys 5) ; J. /a) isk : i 7) ; (3) 
k me 
Set et (Ley egy ay. St} (4) 
ue abs aie alg aly al gh 
k m=O 
where 


Aiak = area of a component k between i and j 


Y and Z.., = thermal response factors of this conponent 
ay ijk 


AI, = surface area of a heavy object k inside i 
i: 


YI... and ZI. = thermal response factors of this object 

J. = mean radiation absorbed on the surfaces of components and objects 
¥ inside i 

a® = overall surface conductance. 


If the heat storage capacity of the component is negligible, 


O O 
= cies = kG 
Yi5k 444k ijk 
m m 
= = f > 1 
oe “nee Op ace We ia = 


where Ki ak is the total heat transmission coefficient of the component. 
Notice here that J. includes solar radiation and radiant energy from 
appliances and others, but that long wave radiation emitted from the 
surfaces enclosing i and the objects inside i is not included. The solar 
radiation incident on the arbitrary surface of the outside component can 
be calculated from normal direct radiation and sky radiation, both of 
which are regarded as input meteorological data. Direct solar gain ina 
room is calculated from the solar radiation calculated above by taking 
account of the shading effects of window glass and shades. If the room i 
is an outdoor space, Js also includes the contribution of nocturnal 
radiation, which is treated as input data. The nocturnal radiation can be 
evaluated from cloud amount and atmospheric vapor pressure according to 
Yamamoto's empirical formula>), if the observed one is not available. Thus 
the same radiation data is input in this procedure as in the program of 
HASP 4) developed by the Society of Heating, Air-Conditioning and Sanitary 


Engineers of Japan. 


(2) Thermal Interaction with the Ground 

The net heat loss by transmission to the component fixed on the 
ground can be assumed as follows. 

If daily fluctuations are neglected, natural soil temperature can be 


approximated by, 


0G(z,t) = As exp (-Bz)* cos (2mt/T, - 62 - >) + OM (5) 
where 

8G = natural soil temperature 

8A = annual amplitude of @G at z=0 

8M = annual mean of 6G 

Z = depth from the ground surface 

t = time 

B = Yu /(aT)) 

To = period of 6G = a year 

a = thermometric conductivity of soil 

> = phase difference. 


If it is assumed that the temperature under the building is equal to 66 


below z=d, G. is given as: 


1 ENS ee n om, _n-m 
) G, = 0,{2G, (8, - 6s ) hia wan: gas : ZG.0E, } (6) 


m=1 


where 


O 
Ge = AAG. / (ZG. +) 


8S” = 6G(0,nAt) 


q. = A,¥t/a,T *OA+cos(2mAt/T. - 6 + 1/4) 
ae 1 ue @) O 
feaeere es — Ge’ / (aac, } 
ak aL 1 


AG, = surface area of the component fixed on the ground inside i 


2G. = thermal response factor of this component 
(but surface boundary layer is excluded and component’s thicknss 
is da) 

he = mean thermal conductivity of this component 

a; = mean thermometric conductivity of this component. 


(3) System of Equations of Room Heat Balance and its Solution 
Substitution of equation (2), (3), (4), and (6) in equation (1) and 


collection of the unknown air temperatures on the left-hand side give: 


4 rn n n 
ede eae By 
where 
: -(ZIA Y cay.) fore. = 4 
ij p ijk i5k ‘an 
- : z De ueeeuyoey er uC ASEATs {210 NT) 
eee yf ott atiy, Re pecs Ft es prada a ik 
n ‘ 0 m n-m n-m O n 
n. = Lf ZA. af z Ysa, (Pu + J, /a) + Yi 5340570 
Techs) er eee. ec en) 
a | aL ale a5 ak 
m= 
© m m n-m qi) 
- ZAI, 1 i? (ZI. = YT 1) (6. +3, /0) 
k m= 
O @) n n “ fay fal n-l n 
CZ et Toe , ee Oe ve 0 Coe” (KE FR 
zk ae rid. ul yen en ee nS 


N = number of indoor rooms (indoor rooms are numbered from-l to N). 


Because Bi is square matrix of order N, equation (7) forms a 
system of linear equations in which the unknowns are the indoor air 
temperatures and can be solved by the successive over-relaxation method >) 
The solution of equation (7) is natural temperature without air 
conditioning. On the other hand, if some rooms are air-conditioned, their 
sensible heat loads can be calculated as the residuals defined in the 
successive over-relaxation method while regarding their air temperatures 


as known. 


2-2 PROCEDURE FOR CALCULATING NATURAL VENTILATION 


(1) Development of a new Approximate Procedure 

Natural ventilation (including infiltration) arises as a result of 
pressure differences generated across openings and cracks by the action of 
external wind and by buoyancy due to air temperature differences. However, 
the relationship between air flow rate by ventilation and pressure 
difference is nonlinear for an opening or a crack. Therefore, it will take 
much computer time to calculate the amount of natural ventilation in each 
room hour by hour, if the system of equations of ventilation remains 
nonlinear. In this paper, a new approximate computer procedure for solving 
the system of equations of ventilation is developed and applied to the 
calculation of natural ventilation. It is expected that this procedure has 
an appropriate accuracy and an acceptable computer time for the practical 


computation of air-conditioning loads. 


(2) Net Volume Flow of Air 

In the notation of this section, attention must be paid to the 
omission of the superscript indicating the time point. Applying the crack 
method, the net volume flow rate Vi 5k of air through an opening or a crack 


"k" located between i and j is approximately represented by, 


Vijk ~ F554 H 55,14 


ee ak 
5k ahs (8) 


where the expressions of Cast Oak! and "i4k are& given in *Table 1}The 


expression of a crack or an opening is introduced by the following 


approximation: 
Tee vols |e, O-. 
Ros | ay ah (Kap. - gzdo,.) x |AP,. - gzho,. a oe 
= o,,,{AP.., - qzA 
isk ag OT i Te) 


Ngee ce, 
x = Siar ie De, 1/2 
JAP, Go (Zia + Li ay / ) | ijk 


Though equation (8) is nonlinear in the relationship between Mae and 
Oa 5k! it can be approximated by applying,several linear relationships: 


Vijk ech ie 


ijx?igk + Pag) (9) 
mes nhs : 
where ijk and ijk have several values different in the range of aTake 
and for the value of "iik respectively. Here the team "range" means a 
range within which one of several linear relationships expressed by 
equation (9) is applied, and the number of ranges is equal to that of the 
relationships. Five linear relationships were applied to equation (8) for 
a value of deste in the actual computer program. 
(3) System of Equations of Ventilation Network and its Solution 
The balance equations of room air volume is described for the room 


i" as follows: 


M 
LAST Sed, Ded O (10) 


Wei ee 
Substitution of equation (9) in equation (10) and reference to the 
expression for  i5k in Table 1 lead to the system of equations of the 


natural ventilation network in which pressures are unknown, 


N 
nad a ae ea: | (iy 
j=1 i (| a 
where 
i5 = Bay F phON aplia Mee J 
M 
hag = yee 
j=l * 
M M 
i, =" & (to". ba ghiek Lo m ) + Dy 
; ijk ik yo a eG cis 
jal eee td 1j, 1jk ijk ijk SEN s 
eS P43 AGK 
6 = 
Lak ijk ee 745k i anes 


Bhavel iy. 


Table 1 &E i 
xpressions of Sankt Hee ijk 


horizontal 
crack 


vertical 


crack 


Definition of Terms 


a. = infiltration constant per unit) lenth of 
crack 


= gravitational constant 


eh, = length of®crack or opening (any opening is 
assumed to be rectangular) 


Nese = flow exponent of crack or opening 
1) 
Pe eS state pressure on the ground level 
1 
.., = area of opening 
ijk 
“eeeeeneront of crack or opening 
ijk : 
tins F discharge coefficient of opening 
13 
APS as Pim 2 
Ly 5 au 
Pas = air density difference between i and j 
oe AS lea n-l 
S 8. au a3* 
Pp (8, F Aches. 
Po = mean air density (constant) 


* This assumption is not elaborate, but is inevitable 
to make calculation procedure as simple as possible. 


The pressures of outdoor spaces (i=N+1~M) are taken as the boundary 


conditions and are given by, 


2 
= 2 
Phe CD, p Au 


where 
CD, = wind pressure coefficient 
U = wind velocity at the roof of building. 


(12) 


The wind pressure coefficients are evaluated from an empirical function 


which depends on an angle between wind direction and normal of outside 


surface. This function, which is shown in Fig.2, was formulated on the 


basis of Sekine's experimental data®) 


wind pressure coefficient 


outside component 


outdoor 
space j 


wind 


Tewrou 


Fig.2 Wind Pressure Coefficient Adopted in Calculation 


of Natural Ventilaion 


- 10 - 


Equation (11) can be solved by the Choleski method”) because the 
matrix ary is symmetric. The solved pressures are substituted into the 
expression of “i45k" which results in new values of ijk: The ranges 
re-defined by the new values of “i5k have to be compared individually with 
the former values. Disagreement between the former ranges and re-defined 
ones requires the re-solving equation (11), whose coefficient matrix and 
constant vector are re-calculated. Thus this process must be iterated 
until the agreements between both ranges are completely reached. In this 


way the final solution of pressures can be attained. 


(4) Calculation of Volume Flow Rate of Air 


The volume flow rate from j to i is given by, 


ae a Bake (13) 
k 
where Bask is volume flow rate which goes across k into i. Applying the 


final pressure solution, Diag can be calculated in the following four 


cases restricted by two factors: 


ht ela Sere 


1) [Ao ; ;<0 and AP; 3205 5x! or [Ao ; ,20 and AP. > ] 


j=" ijk 
) 


: Pape hb 
eee | tisk desk ijk 


2) [Ap 5 ,<0 and ¥ 55x64? 55<¢ ] 


1 ae big 
Oi 4k n +1 
Ss eee ee Shee 2 = + 

ae. = , (AP... ct eae we SB 

+ 

ijk eee BACLT ees ere 5 ijk 
6) ares aA 

3) [Ap ; ;<0 and AP 5 554 i5x! or [Ao ; 52 and ee 


ak 


4) [do ; .>0 and $5 4p <AP i SVs 44) 
Pigk 


u A = 
1jk ee ae neh Ee, 8 
J Sone 1Ghp Lee 


ee Site ee L 
Haare - O44K Ny tere 


sats Sole 


As soon as 435 is calculated, the volume flow rate from i to j is 


obtained by, 


MEE ae her (14) 


Thus the volume flow rates calculated by this procedure are applied 


to the room heat balance equations. 


= l2+- ’ 


4 


3. VALUATION OF DEVELOPED PROCEDURE 


(1) Purpose and Process 

It is necessary to valuate the accuracy and the validity of this 
calculation procedure by an appropriate process, because it includes many 
physical assumptions and many mathematical approximations. In this paper 
the valuation was attemped by comparing the prediction by this procedure 
with the observation in actual dwelling houses. This valuation process 
seems to be most reliable to valuate general accuracy and validity easily. 
In fact, room air temperature was adopted as a comparable quantity for the 
valuation. It is the main reason for this adoption that room air 
temperatures can be measured more easily and more accurately than room 


heat loads. 


The calculation procedure was coded with Fortran language, and the 
abbreviated name of BRIMAP* was applied to this computer program. The 
observation on room air temperature was made in two different types of 
dwelling house. Both are experiment houses and were built using the 
conventional construction method of Japan. One is a house without thermal 
insulation,which is called by the name of House-A. The other is a well 


insulated house named House-B. However, both have the same features except 


machine 
room 


8,190 


6,370 


Japanese 
room 


living-dining 
room 


5190 balcony 


first-floor plan second-floor plan 


Fig.3 Plan of House-A and House-B 


ie 


* Building Research Institute/Multi-Room Air-Conditioning Load Program 
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for the insulation of the outside wall, the floor, and the ceiling and the 
material of the window frame. In addition, both were built adjacently on 
the same site of BRI at Tsukuba, Japan. The plan of both houses is shown 


in Fig.3, and the differences between both are listed in Table 2. The 


Table 2 Difference of Component Specification between House-A 


and House-B 


(unit of thickness is mm) 


component House-A House-B 


outside 
wall 
2F 

ceiling 


1F-floor 


window 


NNANANARRARARART 
MOXOS 


mortal 30 
wooden board 8 
air space 


plaster board 12 


wooden board 8 
air space 


glass wool 50 
plaster board 12 


room room 
wooden board 5 glass wool 100 
insulation wooden board 5 


insulation 


poana 12 hoara fa 


SOOO OOS 


Pern 4ypee 


MINN 


KA 


SS 
A 


Sk 


208, 
aN 


carpet 4 carpet 4 


wooden board 12 
glass wool 50 


single pane with] single pane with 
wooden frame aluminium frame 


wooden board 12 


S ha-+ 


observations were taken under the two kinds of living conditions 


simulating an actual living pattern. One was to heat some rooms 


intermittently with fan-heaters. The other condition was to open and close 


doors, windows, and curtains. Both conditions were common to both houses. 


The time schedules for both conditions are illustrated in Fig.4, which is 
accompanied with the meteorological condition observed in the vicinity of 


both houses. The meteorological condition consists of atmospheric 


temperature, direct solar radiation, sky solar radiation, and 


wind vector. These conditions of living and weather were input for the 


prediction. The minimum unit of time was half an hour in both observation 


and prediction. 


heating up living-dining room (3kw) —=— 
heating up Japanese room (2kw) 

heating up bedroom-A (3kw) 

‘heating up bedroom-C (2kw) — 
opening all curtains 
opening the door between l.d. and Jpn ecu SED 
opening doors and windows for cleaning = 


o4---- ee oe -- 


0 


ie 6 12 time 18 6 12 time 18 
< » 
eee Be 5 eee a ee r 
atmospheric s Feb. 18 ‘81 , Feb. 19 '81 
© temperature = Pa 
U 6 
p direct solar 2 = 
radiation _ s = a 52 ee eee 
@ 
% sky solar p 4 
radiation i 
en 
nocturnal 
e radiation 5 2 
as] wind vector 
0 Sm/s wind vector @: 
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SS 


Fig.4 Living and Meteorological Conditions in the Observation 
and the Prediction 


eT 


(2) Results 

Fig.5 represents the comparisons between prediction and observation 
in some room air temperatures of both houses. The observed room air 
temperatures in this figure are the temperatures which were defined by 
averaging the profile formed of the temperatures observed with seven 
thermocouples lined vertically. 

An agreement between the prediction and the observation is fairly 
good and seems to be satisfactory for the practical calculation of air 
conditioning loads. The difference between House-A and House-B is well 
simulated. This difference is brought about by the insulation with glass 
wool. The sudden descents of temperature at 9:30 a.m. every day, which 
were caused by opening most windows and most doors for cleaning, are also 
simulated suitably. However, considerable disagreements between the 
prediction and the observation are found here and there, where the room 
was heated. These disagreements are considered to be caused by ignoring 
the heat exchange by long wave radiation among the room surfaces and the 
vertical distribution of room air temperature in the calculation 


procedure. 
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Fig.5 Comparison of Room Air Temperatures between the 


Prediction and the Observation 


-17 - 


4. CONCLUSIONS 


This paper presented a calculation procedure for the multi-room 
problem including the effect of natural ventilation. The calculation 
procedure for natural ventilation is a new approximate procedure developed 
according to the crack method. The calculation procedure for multi-room 
temperaures or heat loads is based on the response factor method. 

In order to valuate the accuracy and the validity of this procedure, 
the room air temperatures predicted by it were compared with those 
observed. Generally speaking, an agreement between the prediction and the 
observation was good and acceptable for practical calculation. However, 
there were some disagreements in the period of heating the rooms. These 
disagreements suggest that it is necessary for more accurate simulation to 
develop the calculation procedure including the effect of the heat 
exchange by long wave radiation and the vertical temperature distribution 
in a room. Nevertheless, this calculation procedure is considerably valid 
and available for seasonal heat loads. The computer program BRIMAP has 
been utilized for designing and evaluating the energy conservation in 


houses of Japan. 
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